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Abstract 


Observer-based sliding mode control is investigated for application to aircraft 
reconfigurable flight control. A comprehensive overview of recon flgurable flight 
control is given, including a review of the current state-of-the-art within the 
subdisciplines of fault detection, parameter identification, adaptive control schemes, and 
dynamic control allocation. Of the adaptive control methods reviewed, sliding mode 
control (SMC) appears very promising due its property of invariance to matched 
uncertainty. An overview of sliding mode control is given and its remarkable properties 
are demonstrated by example. Sliding mode methods, however, are difficult to 
implement because unmodeled parasitic dynamics cause immediate and severe 
instability. This presents a challenge for all practical applications with limited 
bandwidth actuators. One method to deal with parasitic dynamics is the use of an 
asymptotic observer in the feedback path Observer-based SMC is investigated, and a 
method for selecting observer gains is offered An additional method for shaping the 
feedback loop using a filter is also developed. It is shown that this SMC prefilter is 
equivalent to a form of model reference hedging A complete design procedure is given 
which takes advantage of the sliding mode boundary layer to recast the SMC as a linear 
control law. Frequency domain loop shaping is then used to design the sliding manifold. 
Finally, three aircraft applications are demonstrated. An F-18/HARV is used to 
demonstrate a S1SO pitch rate tracking controller. It is also used to demonstrate a 
MIMO lateral-directional roll rate tracking controller. The last application is a full linear 
six degree-of-freedom advanced tailless fighter model. The observer-based SMC is seen 
to provide excellent tracking with superior robustness to parameter changes and actuator 
failures. 



Chapter 1 


Overview 

1.1 Introduction 

Reconfigurable flight control is an area of research that has seen exponential 
growth within the control community in the last two decades The early works in 
reconfigurable flight control— just a few research papers— began to appear in the 1960’s 
and 1970’s The 1 980’s saw a growing interest with a few dozen papers, and the 1990’s 
saw an explosion of effort in reconfigurable flight control There are hundreds of 
reconfigurable-related papers from the 1990’s, many of them from the latter part of the 
decade. Considerable progress has been made, but the problem is a difficult one and 
continues to see increasing interest on the part of researchers worldwide. 

It is difficult to find a precise definition in the literature of reconfigurable flight 
control. There is an entire body of work in the area of adaptive control in which the 
term reconfigurable never appears. The basic purpose of an adaptive control law is to 
control a system with unknown and/or time-varying parameters However, this is 
exactly what a reconfigurable control law is intended to do. In fact, all the 
reconfigurable control approaches use some type of adaptive control scheme or adaptive 
control redistribution Strictly speaking, a reconfigurable control law is adaptive 
Likewise, an adaptive control law is one that reconfigures itself On the surface, 
reconfigurable and adaptive seem to be synonymous terms. However, there does appear 
to be a difference between the terms when one considers what each one implies 
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Reconfiguration in the literature always implies a response to some type of system 
failure or damage. This failure may be the loss of a control effector (i.e. a control 
surface locked in a hard-over or zero position, a broken linkage allowing the control 
surface to float, degraded controller effectiveness, or degraded actuator rate response) 
The failure may be damage to an aerodynamic surface resulting in large changes to the 
stability and control derivatives. The failure could also be related to the thrust or thrust 
vectoring or could be simply a failed sensor. Whatever the failure type, reconfigurable 
flight control implies that a sudden, large, unknown failure has occurred. Adaptive 
control does not always connote this kind of failure. Changes to system parameters are 
generally implied (or explicitly stated) to be smooth and non-catastrophic — usually due 
to mechanical wear, unmodeled dynamics, and structured uncertainty due to 
environmental variations throughout the operating envelope. With this in mind, an exact 
definition of reconfigurable flight control is offered Reconfigurable control is an 
automatic control system which is able to compensate for sudden, potentially large, 
unknown failure events in real-time using on-line adaptive control laws and/or adaptive 
redistribution of control effort with the objective of guaranteeing system stability and 
achieving some level of required performance and handling qualities. 

The motivation for developing reconfigurable flight controls is clear. Failures 
during flight are inevitable — especially in combat aircraft If the flight control system is 
capable of stabilizing the aircraft and providing acceptable control, it may be possible to 
return the damaged aircraft to base, salvage the airframe, and save lives “Lessons 
learned from the Vietnam Era show 20% of aircraft losses were due to flight control 
damage. Loss of hydraulics, actuator damage, and surface damage are responsible for 
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most of the flight control system losses.”' In addition, high performance fighter aircraft 
tend to be so inherently unstable that a pilot is incapable of controlling it without the 
stabilizing influence of the control system. This, combined with the increasingly 
complex workload on the pilot, dictates the need for a highly fault-tolerant control 
system Commercial transport aircraft can also profit from the potential safety benefits 
of reconfigurable flight controls. “In the majority of cases surveyed, major flight control 
system failures .. have resulted in crashes, with a total of over 1200 fatalities ” 2 

The approach for designing a reconfigurable flight control system has undergone a 
fundamental shift in philosophy in the last five to eight years Initially, the concept was 
to design a control law robust enough to guarantee system stability in the face of the full 
range of potential failures. Then, upon the event of a failure, the system would identify 
the fault, isolate the failure mode, estimate new system parameters, and select new 
control law gains and/or control input distribution to regain some level of performance. 
The schedule of gains for each condition was calculated off-line. Today, most 
reconfigurable control schemes being proposed are fully adaptive — continuously 
reconfiguring themselves using dynamic on-line algorithms which do not require any 
a priori schedule of gains. Some perform continuous estimation of system parameters. 
Others do not even require this, deriving adaptive laws based on performance-related 
measures only. In the words of one set of authors, the control system is “self- 
designing ” 3 This shift in philosophy is leading to classes of controllers that are highly 
robust to both “normal” parameter variations as well as damage/failures. They are more 
flexible, cost efficient, and general in nature. It is almost certain that future advanced 
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aircraft flight control systems will incorporate some form of on-line reconfigurable 
control law 

1.2 Thesis Overview 

The purpose of this research is to demonstrate a design strategy for a multi-input, 
multi-output (MIMO) flight control system which is: 

• able to provide satisfactory performance for a coupled-axis, statically 
unstable, high performance aircraft with a highly redundant control effector 
suite 

• robust to model uncertainty, smooth parameter variations, and sensor and 
process noise 

• able to adapt to sudden, unknown, potentially large system parameter 
changes due to a failure event 

• able to provide satisfactory handling qualities and present no tendencies for 
pilot-induced oscillation (PIO) before and after a failure event 

• based on current approaches utilizing sliding modes. 

First, an overall view of reconfigurable flight control is given. Chapter 1 contains a 
comprehensive look at the current state of all key aspects of reconfigurable control. 
Based on this review, the sliding mode control method is chosen as a key focus area 
Chapter 2 provides a brief introduction to sliding mode control. Chapter 3 investigates 
potential solutions to key implementation issues associated with sliding mode 
controllers. Based on this work, a tutorial design procedure along with several 
application examples is given in Chapter 4 Conclusions and recommendations for 
future research are found in Chapter 5, 
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1.3 Overview of Reconfigurable Flight Control 

There are four main aspects of flight control reconfiguration that must be 

addressed for a complete treatment of the subject. 

• Failure detection 

• System parameter estimation 

• Flight control law reconfiguration 

• Control allocation 

These four subdisciplines have traditionally been treated as separate problems because 
each has fundamentally different objectives. However, recent works 4,5 are addressing 
the interaction between them and demonstrating the need to consider the system as a 
whole. For example, it is shown later that the control allocation process can create a 
condition in which system parameter identification becomes impossible without 
considerable preconditioning. In the sections that follow, each of these components of 
reconfigurable flight controls are introduced individually. In this respect, the following 
introduction follows the traditional approach. However, the problematic interactions 
between the components (where they are known to exist) are also introduced. 

1.3.1 Failure Detection 

In the time-sequence of events that occur after a failure, detection is logically the 
task that must be accomplished first. This task has traditionally been the job of the pilot 
It is up to the pilot to recognize that something is wrong and to initiate efforts to 
compensate for the failure. While this method may not be suitable for high performance 
applications, it still finds its place in practical applications. In a proof-of-concept flight 
test program, 2 an MD-11 is flown using a propulsion-only emergency control system. 
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To activate the system, the pilot essentially flips a switch. This requires the pilot to be 
the one to detect and isolate the failure. Although this approach works well in this 
application, automatic failure detection is often required. 

While little emphasis is placed on fault detection in the current literature, the early 
efforts in reconfigurable controls include fault detection as a primary concern A key 
component in the Self-Repairing Flight Control System (SRFCS) 1 study, conducted by 
the Air Force Wright Aeronautical Laboratory (AFWAL), McDonnel Douglas, General 
Electric, and Alphatech 6 in the early 1980’s, is their so-called Failure Detection and 
Isolation (FDI) procedure. This includes a local FDI algorithm to detect actuator failures 
and a global FDI algorithm to detect surface failures. The global FDI compares the 
outputs of a nominal system model and the actual aircraft This error signal is passed 
through a set of hypothesis-testing filters. The output of each filter represents the 
statistical likelihood that the failure hypothesized by the filter has occurred. This 
approach was flight tested on a NASA F-15 and was shown to have “high potential for 
the concepts evaluated.” 7 In the same time period, the Control Reconfigurable Combat 
Aircraft (CRCA) program, conducted by Grumman, Lear Astronautics, Charles River 

o 

Analytics, and AFWAL, was using a similar approach Again banks of Kalman filters 
operate in parallel and test a predetermined set of failure hypotheses. The CRCA 
program concluded with piloted simulations and flight tests in which the method 
performed well, but extensive adjustments were needed throughout the testing. This 
multiple model approach has the advantage of being able to provide very fast and 
efficient detection of failures among a set of preplanned conditions Unfortunately, as 

the number of hypothesized failures grows, it becomes increasingly difficult to classify 
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all the possible combinations of failure conditions This leads to either an unreasonably 
large set of classifications or an unacceptably large number of gaps in potential failure 
modes. Also, since failure detection relies on models of the nominal, unfailed system, 
any differences between the actual system and the nominal system model can result in a 
false detection of a failure False alarms are a major difficulty with this method. 
Despite these difficulties, the basic idea of this method is further developed in recent 
works and is termed the Multiple Model Adaptive Estimation (MMAE) method y 12 
MMAE, which encompasses a entire control scheme has been successfully flight tested 
on the VISTA/F-16 

Another approach appearing in the current literature for failure detection involves 
the use of the sliding mode control technique. In particular, a sliding mode observer is 
designed with either an additive perturbation parameter 13,14 or a multiplicative 
parameter. 15 This parameter is then adaptively estimated. Using thresholds on this 
parameter, actuator failures are determined. 

Probably the largest reason for the lack of extensive work in fault detection is the 
fact that many of the adaptive algorithms being proposed do not need it The algorithms 
either continually estimate system parameters (and thus automatically “see” a failure) or 
do not require the parameters explicitly at all. Most current works do not even make 
mention of fault detection. Those that do require knowledge of a failure event typically 
assume that this knowledge is provided to them “somehow.” 

1.3.2 Parameter Estimation 

System parameters are those elements of a dynamic system model which define 
its dynamic response In linear state space where x = state variables, u = input vector, 
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and y = output variables, the system parameters are the components of the A, B, C, and 
D matrices, as shown in equation (1.1) 

i(i) = A x(t) + B«(t) 

( 11 ) 

j>(t) = Cx(t) + D«(t) v ’ 

For aircraft, these components are the stability and control derivatives. System 
parameters might also be variables in transfer functions (e g., time constants, damping 
ratios, natural undamped frequencies, etc.) as in equation ( 1.2 ). 

0 (s + 1/T 6 ,Xs + I/T 6 ) 

— - ( s ) = is- 0 g — 2 2 2 — 

e (s + 2q S p© S pS +o> S p)(s + 2qph©ph s ( 1-2) 

Of course, the actual plant parameters are never known exactly, even with careful 
modeling. Unmodeled high order dynamics, linearization approximations, rigid body 
approximations, operating envelope effects, axis coupling, control effector limits, 
hysteresis effects, and stochastic randomness may all be present in the actual plant. 
However, the model used for designing the control system can not encompass all these 
effects perfectly This is a common problem for all control designers and is addressed 
by careful attention to modeling the dominant plant dynamics and concerted efforts to 
ensure controller robustness. A key point in a traditional design problem is that the work 
of identifying the system parameters is all done off-line. The control algorithm may be 
required to switch through a schedule of gains, but all these gains are calculated a priori 
based on off-line estimates of the plant parameters. 

In a reconfigurable control design, the controller must be able to handle large, 
sudden, unknown changes to system parameters These changes greatly exacerbate the 


design problem Some adaptive design methods require knowledge of the plant; but 
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after a system “failure,” the new effective plant is unknown. This leads to the 
requirement for on-line system parameter identification. There are many parameter 
estimation techniques in use; however, in reconfigurable control applications, the 
requirement of the algorithm to nin in real time narrows the choices of algorithms to 
only the most computationally efficient. Some on-line parameter estimation techniques 
include: Lyapunov design approach, gradient methods (with instantaneous cost 

functions and integral cost functions), pure least squares, pure least squares with 
covariance resetting, modified least squares with forgetting factor, gradient methods 
with projection, least squares with projection, and hybrids. 16 The most common ways of 
performing on-line system parameter estimation in the current literature include: time 
domain least squares parameter estimation, frequency domain parameter estimation, 
reduced-order transfer function estimation, and neural nets. Each of these methods are 
introduced in the following sections. 

All parameter estimation techniques are faced with two difficult problems' noise 
and data information content. 17 Noise poses a large obstacle in accurate parameter 
estimation because it is difficult to distinguish between noise and a sudden change in a 
system parameter. In the time domain, large data records are required and special 
techniques, like “forgetting factors,” are used Data information content is an especially 
prevalent problem in aircraft applications. Aircraft frequently spend large amounts of 
time in a steady state condition with relatively constant state and control variables This 
lack of persistent excitation leads to a lack of information content in the signals used by 
the parameter estimation algorithm and can cause the calculations to become ill- 
conditioned In the most extreme case, if the input signal u = 0 Vt > 0, it is clear that no 
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useful information will be available for the identification algorithm. In fact, it is well 
known that it is impossible to establish parameter error convergence to zero without 
imposing some kind of conditions on the input to the system. 16 Another effect that leads 
to a lack of information content is data colinearity due to the control system or control 
allocation routine. 18 ' 21 Unfortunately, when the states and controls are nearly 
proportional to one another, or if the control deflections are linearly correlated, it is 
impossible to identify individual stability and control derivatives from the measured data 
alone. 17 ' 18 These two data content problems— persistent excitation and data colinearity— 
are major issues in parameter identification. 

1. 3.2.1 Time Domain Least Squares Parameter Estimation 

Clearly the most commonly used method for system identification is some 

variation of a least squares approach in the time domain. There are many references 
which utilize this approach ’ ’ ' The basics are as follows: 

Consider k= 1 . . .n discrete measurements of a scalar signal y(k), and a regressor 
vector w(k), where 

y(k) = w T (k)0 (13) 

and <?is the vector of unknown parameters to be determined. For a linear system in state 
space, w consists of stacked state variables. At time instant n, a least squares criterion is 
defined: 

n 2 ? 

J[6>(n)]= £ |y(k)-wT(k)0(n)j A, n- k+a(#(n)-6>(n-l)j (1.4) 

When a = 0 and X = 1, this is the standard least-squares parameter identification scheme. 


The constant X is the so-called exponential forgetting factor and is often included to 
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minimize the cumulative effects of noise. When the constant a ^ 0, variations of the 
parameter vector B from one time step to the next are penalized. This is an example of 
one of several common modifications which can be made to the least squares approach. 
This optimization problem has a well known standard solution which can be given in 
“batch” form or in regressive form. This formulation is able to follow time-varying 
parameters and is a popular starting point for most reconfigurable control efforts. The 
main problem that plagues this approach is the need for persistently exciting input 
signals. One solution to this is active noise injection, called dithering. This, however, 
can have undesirable effects upon handling qualities. 

I.3.2.2 Frequency Domain Parameter Estimation 

Another approach is to frame the problem as a least squares error minimization in 

the frequency domain. 17 This algorithm is typically able to converge to the correct 
parameter values in less than one period of the dominant dynamic mode. The basics of 
the technique are as follows. 

Take the Fourier transform of the linear, time invariant state space equations to obtain 

jco x(gj) = A x(<a) + Bm (©) 
y(oo) = Cx (co) + Du (©) 

where A e 9l nxn , B € 9T xm . Next, form the cost function for the r 1 * state equation of 
the vector equation ( 1.5). This is given in ( 1.6 ). 

f ' 

Jr = I J®k x r (k)-A r x(k)-Bu(k)| 2 ( 1.6) 

k=l 

where A r and B r are the r 111 rows of matrices A and B; x r (k) is the r th element of the 
vector x for frequency ©r; and x(k) and u(k) are the Fourier transforms of the state 
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and control vectors for frequency Ok These transforms are calculated recursively using 
the discrete Fourier transform in the following manner, 

X|(<o)=X,Al 

where At is the sampling interval, and iAt represents the current time at sample i. 


Define the vector of unknown model parameters in the r‘ h row as 

Mk, a ,,2 - A,, B,., B, , ••• B,J T (1.8) 


Also, define 



jffl, x r (l) 


'x T (l) 

5 T (1)' 

Y = 

j© 2 x r (2) 

: 1 

x = 

? t (2) 

*J T (2) 


IX 

i 


_x T (f) 

u T (f). 


Then, the problem can be cast as a standard least squares regression problem with 
complex data, 


Y = X0 + g (MO) 

where e is the complex equation error in the frequency domain. The least squares cost 
function in equation (1.6) now takes the form 

J r = i(Y-X0)*(Y-X0) (1.11) 

Finally, the parameter vector estimate which minimizes this cost function is given by 

0, = [Re(X*X)f ' Re(X*Y) ( 1.12) 
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Note that 9 r e ^ (n+m)xf ; and equation ( 1.12 ) must be solved for each row of the state 
vector each time the routine is called. Fortunately, the matrix X X is the same for each 
row, so the inverse [Re(X*X)] ' e on |y nee( j s 5 e calculated once per pass. 

The existence of this inverse is not guaranteed, and in fact, Re(X X ) is nearly singular 
for the first few time steps due to a lack of information content in the data during the 
initial seconds. 

There are several advantages which make this method attractive: 

• Using a carefully chosen limited frequency band for the Fourier transforms 
restricts the analysis to the band of frequencies where the system dynamics 
occur. This automatically filters out all unwanted high frequency data due 
to noise and higher order system dynamics. Also, excluding zero 
frequency eliminates trim values and measurement biases, so it is not 
necessary to estimate bias parameters. For aircraft (where the rigid body 
dynamics typically lie in the band of approximately 0.01-1.5 Hz), the 
frequency band which seems to work well is 0. 1-1 5 Hz. 

• The algorithm requires no initial values for the parameters and does not 
need to be regularized with a priori values or constraints. 

• The algorithm is robust to noisy state measurements and data dropouts 

• The algorithm does not require persistent input excitation; and it is not 
upset by intermittent input excitation. 

• The algorithm has fixed memory requirements regardless of data record 
length. 
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• The method has low computational requirements due to the recursive 
Fourier transforms and the simplicity of the algorithm 

After trying to duplicate this work, 17 some additional comments about this method 
are in order. Some of these findings present obstacles which must be overcome if this 
method is to be used in a reconfigurable control setting 

• As with any parameter estimation algorithm, the stability of convergence of 
the parameter values is highly dependent on the control input excitation 
signal. This signal must be frequency-rich in order to achieve acceptable 
results. Inputs such as a step, a sinusoidal doublet, or a series of sine 
waves with several distinct frequencies, produce reasonable results as long 
as they excite the dominant system modes. However, more frequency-rich 
inputs produce better results. Also, the excitation must be of sufficient 
duration to achieve convergence Experiments with short, periodic pulses 
yield very poor results. 

• The convergence of the parameter values are also dependent on the range 
of frequencies chosen for the Fourier transforms. Both the range of 
frequencies and the fineness of the frequency spacing affect the solution 
convergence. Since both of these are arbitrary, it is not clear how to select 
the optimal frequency band and spacing. Experiments show that the 
frequency range used for an F-16 needs to be completely different than the 
range used for a DC-8 Unfortunately, simply choosing a wider frequency 
band is not necessarily the best answer. A wider frequency range increases 
the computational and memory requirements (for a given bin spacing) and 
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reduces the accuracy of the algorithm Low frequency bins tend to “pick 
up” bias signals High frequency bins tend to pick up noise and high 

frequency modes Ideally, the frequency range chosen should be just large 
enough to capture the dominant system dynarhics being modeled. This 
could be an interesting task in a reconfigurable control setting where the 
system dynamics could potentially have modes considerably different than 
that of the healthy vehicle 

• The algorithm assumes constant system parameters. If a parameter 
changes suddenly (or gradually) during the estimation process, the 
algorithm converges to incorrect values This is because the algorithm has 
“infinite memory," by nature of the recursive Fourier transform being used. 
This problem could be overcome by restarting the recursive Fourier series 
at given time intervals and assuming the parameters remain constant 
throughout this time interval If the Fourier transforms are not reset, the 
parameter estimates will be some average of the changed and unchanged 
dynamics, weighted by the information content in the data. 

• The algorithm assumes zero initial conditions for the states Nonzero 
initial states cause the solution to converge to completely incorrect values. 
When using relatively low frequency bins (as in the case of aircraft 
parameter estimation), any constant bias (zero frequency component) in the 
signal “spills over” into the adjacent frequency bins, thus polluting the 
frequency domain data used for the parameter estimation Non-zero initial 
states appear in the signal as just such a bias. One way to overcome this 
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difficulty is to subtract a low-pass filtered version of the original signal, 
leaving only the frequency components to be used for the parameter 
estimation. The low pass filter introduces a small time lag, but that is not 
an issue because the value being computed is low frequency in nature. 

Overall, this method appears attractive — especially for initial flight test work 39 It 
is interesting that the infinite memory of the Fourier transforms is both a key strength 
and a key downfall of the method. For constant parameters, the infinite memory means 
persistent excitation is not required. Unfortunately, it also means it can not track time- 
varying parameters in its current form. The initial efforts to modify this method to 
accommodate varying parameters are somewhat promising; however, there are 
unresolved implementation issues 

1. 3.2.3 Reduced-Order Transfer Function Estimation 

A method which is attractive in a reconfigurable control setting is reduced-order 

transfer function estimation. Rather than trying to estimate all the parameters for the full 
system model, a lower order model is assumed and its parameters are estimated. This 
concept can be incorporated in both the time domain and frequency domain least-squares 
approaches 25 Also, a more simplistic ad hoc approach can be applied For example 
refer to Figure 1-1 If the system is pulsed with a step input, by measuring outputs like 
rise time, period, maximum amplitude, response slope, and delay time, it is possible to 
estimate the parameters in a standard second order transfer function with delay, as in 
(1.13). 
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Figure 1-1: Reduced Model Order, ad hoc Approach 


P(s) = 


Ko (s + z) e" TjS 
z s 2 +2^co n s + (a 2 


(1 13) 


This approach is used successfully in several references 


40-43 


1.3.2.4 Neural Nets 

Neural nets are receiving a great deal of attention in the current literature and are 
being used in a hugely diverse range of application types. Even though neural nets are 
finding their way into a multitude of applications, there seems to be some resistance to 
their use in systems employing a human operator This is probably due to at least two 
factors. First, due to their highly nonlinear structure, stability convergence is often 
difficult to prove for certain kinds of networks. It is also sometimes difficult to extract 
the “knowledge base” contained in the net after training and to predict results for cases 
outside the training set. Second, even though many networks can be proven to model a 
linear system to any arbitrary degree of accuracy, 44 there is little incentive to favor a 
neural net over traditional linear control methods for which there is a wealth of theory 
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and experience However, in cases where the system is complex and nonlinear, neural 
nets are finding their place. They appear in the literature dealing with reconfigurable 
controls due to the inherent nonlinear nature of the problem Functional link neural nets 
combined with an adaptive critic algorithm are utilized by Cox c/a /. 45 ’ 46 in a system 
identification/adaptive algorithm scheme. Although not exactly a system parameter 
identification application, neural nets are also being used successfully to remove the 
inversion error in dynamic inversion algorithms. 

I.3.2.5 System Identification for Correlated Effectors 

Whether a time domain approach or a frequency domain approach is used, 

correlation of the control effectors is a problem All approaches utilize some type of 
regressor matrix which must be inverted When the control effectors are correlated, 
there will be collinear columns in the regressor matrix — thus making accurate parameter 
identification impossible . 18 There are at least two ways the effectors can become 
correlated 

First is feedback control. Feedback control correlates the effectors’ displacement 
with the aerodynamic angles. One way to deal with this problem is to break the 
parameter identification into two steps First perform a singular value decomposition on 
the regressor matrix, remove all small singular values, and estimate the reduced set of 
parameters. Next, reconstruct the full parameter estimate with some a prion 
knowledge . 18 - 21 ’ 34 ’ 49 

A second way control effectors can become correlated is the use of control 
allocation. Redundant control suites require some means of distributing control 
demands to the multiple effectors Unfortunately, almost all common control allocation 
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schemes result in highly correlated effectors. There are several ways to deal with this 
problem. First, use the extra degrees of freedom to add decorrelating excitation For 
example, when using a psuedo-inverse allocation approach, randomly vary the weighting 
matrix in the optimization cost function 18,19 (see Eqn ( 1 .22 ), pg 34). Another approach 
is to use a priori models in the control allocation and create a reduced set of pseudo- 
effectors. 19 The control effectiveness matrix for these pseudo-effectors is then identified 
rather than for the full effector suite Another alternative is to provide the control 
allocation algorithm with enough information to compensate for the modeling errors 
rather than trying to provide the control law with correctly identified parameters. 19 

Accurate system identification is a complex problem Even under the best 
conditions, identification is difficult and takes time. Then, with the addition of issues 
like persistent excitation and correlated control effectors, an already difficult problem 
becomes much worse. System identification, if required by the control law, is a key 
bottleneck in the reconfigurable control problem. 
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1.3.3 Reconfiguration Control Algorithms 

Adaptive control strategies can be loosely grouped into two categories, indirect 

and direct adaptive control In indirect adaptive control, the process/plant model is 
constructed (usually by means of an on-line observer/parameter estimator), and an 
appropriate control law is then calculated based on this model This is also called 
explicit adaptive control 16 because the design depends on an explicit plant model. Direct 
adaptive control, on the other hand, synthesizes the controller utilizing performance 
criteria only without the explicit construction of the plant model 44 This is sometimes 
called implicit adaptive control 16 because the design is based on an estimation of an 
implicit plant model These definitions sound straight-forward; however, the distinction 
between direct and indirect methods can become easily confused In fact, a direct 
adaptive scheme can be made to appear identical to an indirect adaptive scheme “by 
including a block for calculations with an identity transformation between updated 
parameters and control parameters. In general, for a given plant model, the distinction 
between the direct and indirect approach [only] becomes clear if we go into the details of 
design and analysis.” 16 There are several works which compare different adaptive 
schemes! Bodson etal 22 ' i0 compare three model reference adaptive schemes 1) 
Indirect method, 2) Direct method based on output error; and 3) Direct method based 
on input error They provide a nice comparison, but draw no strong conclusions In a 
later paper, Bodson et al 23 compare four adaptive schemes in light of actuator 
saturation. Steinberg 31 compares seven different non-linear adaptive control laws. 
Again this work provides a good top level overview of the methods, but makes no strong 
conclusions. 
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1. 3.3.1 Indirect Adaptive 

While the indirect controller may appear to be more complex, it does allow some 
separation between the system parameter estimation and the controller synthesis (the 
certainty-equivalence principle). 44 The overall concept of an indirect control scheme is 



Figure 1-2: Indirect Adaptive Controller 


Two indirect methods which receive the most attention in the literature are 
receding horizon optimal control (RHO) and multiple model estimation. Eberhardt et al. 
use a RHO controller on the Innovative Control Effector (ICE) aircraft combined with a 
least squares parameter ID algorithm utilizing a lower order equivalent system flying 
qualities model. 19,25 Pachter et al. use a RHO controller with a one-step-ahead actuator 
rate constraint enforcement. 29 RHO control with a modified least squares parameter ID 
is also demonstrated on an F-16/MATV. 3 In Ward et al . 33 , a RHO controller with least 
squares parameter ID is used in an inner loop, and a polynomial neural net is used for an 
outer loop shipboard landing task of an unmanned airvehicle Maybeck et al. use the 
method of Multiple Model Adaptive Estimation (MMAE). 10,12 Adaptation occurs via a 
control redistribution, and the parameter estimates are provided by the MMAE 
algorithms. Napolitano et a/. 51 also use multiple Kalman filters to estimate the model of 
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the damaged aircraft and use these results in an adaptive control law. Boskovic et at. 52 
recast the aircraft model in terms of a damage parameter. They then use multiple 
observers to estimate this damage parameter. The control law involves multiple 
reference models and output feedback Other examples of indirect control methods are 
model reference adaptive controllers, 53 ' 54 and a simple adaptive PLD controller with 
modified least square parameter ID. Another interesting approach is to form a multiple 
objective cost function which includes optimal control and parameter identification in a 
single cost function. 55 Since parameter identification and control performance are 
competing objectives, this approach allows a direct tradeoff between the two. Rather 
than performing continuous parameter identification, one method proposes taking the 
pilot out of the control loop, inputting control steps/doublets, and estimating lower order 
plant models. 41,42 This model is then used in an adaptive multi-loop scheme with inner 
loop linear dynamic inversion and outer loop quantitative feedback control. 


1. 3.3.2 Direct Adaptive 

Direct adaptive methods almost always include some form of model reference 
following as illustrated in Figure 1-3. 



error 


Figure 1-3: Model Reference Adaptive Control System 
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There are, in general, four elements in the MRAC system. 44 

• The plant, which may be nonlinear, time-varying, with unknown parameters. 

• The reference model, which is usually a lower order linear dynamic model 
which generates a desired closed-loop system output response 

• A controller with time-varying components. 

• Some type of adaptive algorithm which adjusts the controller based on the 
error between the model reference outputs and the actual plant outputs. 

One of the direct methods that has dominated the literature is the method of 
dynamic inversion (Dl). As the name implies, the controller attempts to invert the plant 
in order to cancel its dynamics and then replace these dynamics with those of a reference 
model There are several successful examples of this approach. Dynamic inversion with 
a neural net to regulate the inversion error is demonstrated on the Tailless Advanced 
Fighter Aircraft (TAFA) for the RESTORE program. 24 This approach was successfully 

o * 

tested with full piloted simulations. Others have also demonstrated DI with neural nets 
to remove the inversion error. 34 ’ 47,48 Bacon et al 5G employ non-linear dynamic inversion 
using acceleration and position feedback. Another example of a direct adaptive method 
is an approach called backstepping. Several works utilize this method 57,58 Ferrara 
et al 59 use a classical backstepping approach for the first (n-1) steps, and then use a 
second order sliding mode controller to find the control for the n th step. Other 
approaches include decentralized adaptive neuro-fuzzy design, 60 an adaptive 
proportional plus integral control for the AFTI/F-16 with gain adjustment based on 
errors between the model reference and actual output, 61 and a structured model reference 
adaptive technique in which the kinematic differential equations are assumed to be 
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known exactly— which reduces the amount of uncertainty for which adaptation is 
required. 62 Finally, a method which is receiving notable attention is the method of 
sliding mode control (SMC). Shtessel et al. demonstrate a multi-loop application of 
SMC on an F-16 63 ' 65 and the ICE aircraft. 66 This method provides some of the best 
tracking results of any of the other works, and it does so without any parameter 
identification. Actuator limits are handled by adaptively varying the sliding surface 
boundary layer. Sliding mode control is well known for being very robust (in fact, 
invariant) to certain kinds of uncertainty. This makes it a very attractive choice for a 
reconfigurable control problem. 

I.3.3.3 Others 

There are other schemes which address adaptation to large system failures but do 
not fit well in either a direct or indirect adaptive category. For example, Burken 2 
proposes an emergency augmentation system for an MD- 1 1 transport. A control law 
which utilizes asymmetric thrust capability is designed off-line If the aircraft 

experiences control actuator failure, the pilot turns on the emergency control 
augmentation system. In another work , Burken et al. examine two methods as applied 
to the X-33. In the first method, system failures are modeled as disturbances, and an 
LQR controller is designed to handle these disturbances. No reconfiguration takes 
place — it is simply a robust design. In the second method, reconfiguration is 
accomplished by a quadratic programming control reallocation and a nominal control 

/TO 

law based on the healthy system. Lyshevski presents a robust control law based on 
dynamic programming and Lyapunov to handle system failures Again no active 
reconfiguration takes place. 
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1.3.4 Control Allocation 

A conventional aircraft creates moments about its three axes via three control 
effectors: elevator, ailerons, and rudder. While the ailerons and rudder do not create 
strictly decoupled moments, the three conventional control effectors can sometimes be 
loosely considered decoupled for design purposes. Control allocation for this kind of 
configuration is not an issue. A pitch moment is created exclusively by the elevator; a 
roll moment by the ailerons; and a yaw moment by the rudder. If the coupling of the 
ailerons and rudder is not negligible, the control allocation problem becomes one of 
scheduling the ailerons and rudder such that a decoupled moment is generated (e.i. roll 
moment without yaw). This is traditionally accomplished mechanically by the use of 
aileron-rudder interconnects. 69 

With the advent of sophisticated, highly unconventional aircraft with a large suite 
of coupled control effectors, the problem of control allocation has become a significant 
design issue. The aircraft shown in Figure 1-4 is a conceptual design for which 
considerable work is currently being done in the literature. 24,30,34 ' 35 ’ 47 ' 52,70 Note the 
multiple, coupled control effectors available to generate moments. 



forebody blowing 


axi-symmetric pitch/yaw 
thrust vectoring 

trailing edge flaperons 

— ailerons 


passive porosity strip 
(directional/roll control) 


Figure 1-4: Boeing Tailless Advanced Fighter Aircraft 
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In order to create a desired moment about a given axis with an unconventional design, 
the moment demand must be divided up among several control effectors according to 
some control allocation scheme. Of course, issues like relative control power, preferred- 
positions, and rate and position saturation must all be considered. Control allocation 
often becomes an issue in reconfigurable control applications because it is often a 
control effector which is assumed to be failed or degraded. This, then, can necessitate an 
on-line re-allocation of the controls Also, highly redundant control effector suites are 
typically assumed in reconfigurable applications because some level of control 
redundancy is required in order control the aircraft after a failure 

In general, the problem to be solved in control allocation is the following: 

B u = m d ( 1 14 ) 

where B is the 9f nxm control power derivative matrix, m& is the desired moment vector, 

and u is the control vector for which a solution is sought 4 Since control allocation is 

used for systems with multiple control effectors, it is assumed that the dimension of u is 

greater than the dimension of Wd. Note, the “moments” in the moment vector need not 

all be actual moments in this context. It is understood, without loss of generality, that 

“moments” are taken to mean “the desired control effect.” In the unconstrained case, 

Equation ( 1.14) has, in general, an infinite number of solutions. However, when 

control constraints (position and rate limits) are considered, there may exist a unique 

solution, no solution, or an infinite number of solutions. The determination of the 

control vector u, given the desired moment, is the linear control allocation problem 

Control allocation is, by itself, a large and challenging field of interest and is 

receiving considerable attention in the literature Several of the most common 
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approaches are ganging, direct allocation, linear and quadratic programming, pseudo- 
inverse, cascaded generalized inverse, and daisy chaining. Each of these, as well as 
some of their variants, are introduced below. 

1.3.4.1 Ganging 

The simplest control allocation is a static scheme. The controls are “ganged” 
together by some predetermined static distribution matrix. Due to its simplicity, this 
scheme is often used when the main focus of the research is the demonstration of some 
control law The weights in the allocation can be somewhat arbitrary, although certain 
actuators can be favored for some given reason This method does not optimize any 
dynamic cost function and can not span the entire attainable moment space. However, 
its simplicity makes it an attractive approach in many instances. There are a number of 
examples of recent works which use ganged controls. 27,2!i ’ 40 ’ 41 ’ 47 ’ 48 ’ 61 ’ 71 ’ 72 

1.3.4.2 Direct Allocation (Attainable Moment Subset) 

The Direct Allocation method is based on the Attainable Moment Subset (AMS)— 

the set of all “moments” which can be produced by a set of control effectors constrained 
to move within their given limits, u, < u; < Ui 69 Note that this is a position limit 

and does not account for rate limits. The AMS is found by mapping the w-dimensional 
control space into the desired moment space. For the 3-moment problem, Durham 69 ’ 73 
gives a geometric interpretation of this subspace as a 3-dimensional closed surface 
polytope with vertices, edges, faces, and facets. There are 2 m2 m\\2\(m-2)\] facets in the 
control space Since most of these facets map to the interior of the AMS, only m(m-l) 
facets lie on the exterior boundary. Each facet is a parallelogram with four corners 
called vertices and four sides called edges. There are m(m- 1)+2 vertices in the AMS 74 
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First, the AMS must be calculated for the given system Then, the solution on the 
boundary of the AMS is found in the direction of the desired moment. This represents 

the maximum attainable moment in the given direction, and the associated control vector 

* 

is given the designation u d If the desired moment lies within the AMS, the desired 
moment is attainable, and the control vector ua is found by scaling the boundary solution 
down while preserving the desired direction, according to 


u d = u* d 



( 1.15) 


a 

where “a” is found by solving equation (1.16 ) 73 If the desired moment lies outside the 
AMS, the desired moment is unattainable. Note that since this method defines the entire 
attainable moment space, if the moment lies outside the AMS, it is unattainable by any 
allocation method within the given control effector constraints. Solutions that do not 
violate effector constraints are called admissible, solutions that do violate effector 
constraints are called inadmissible. If the solution turns out to be inadmissible, the 
solution is taken on the boundary of the AMS in the desired direction. This method 
turns out to be computationally intensive due to the difficulty of finding the AMS 
boundary solution Because of this, there are a few different approaches to the actual 
application of this method. Three are introduced here 

1.3.4.2.1 Sequential Facet Search 

The most difficult part of the Direct Allocation method is identifying the point of 
intersection of a line in the direction of the desired moment with the bounding surface of 
the AMS. The brute force way to find this intersection is the sequential facet search. A 
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facet is defined by three vectors, m \ > m i — j ^ and At the intersection of a facet with a 

line in the direction of the desired moment, nii , the simple vector equation is: 

* * * 
a w c j = m l + b/??j_ j +c Wj_j< 

or 




* * 
m x : 


1 



_c 

a 


( 1-16) 


During the process of identifying all the facets on the AMS, the solutions of 

I * * * j-1 

m \ : w,_ j : Wj-kJ are calculated and stored for each facet. Then, the values of a, b, 

and c can be evaluated for a given m<s. A candidate facet is then selected and tested to 
see if it intersects the line in the desired direction If it does, it will satisfy the 

73 

conditions 

a>0, 0 < b < 1, 0<c<l (117) 

If the candidate facet is not the correct one, the next one is checked, and so on. 
Unfortunately, this search can be unreasonably large. For example, with 20 controls, 
there are almost 50 million facets in the subset of constrained controls. 75 Of course, the 
search for the correct facet is narrowed considerably if done in a selective manner 
Durham offers a more efficient algorithm to perform the facet search It is 
considerably better than sequential searching from a random location, but it is still 
computationally cost prohibitive for real time on-line implementation 
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1.3.4.2.2 Edge-Bisecting Facet Search 

Recently, Durham presented another more efficient method for locating the 
desired facet 76 In this method, the AMS subset is first transformed such that the 
direction of the desired moment is aligned with the x-axis The edge of the 
2-dimensional figure that crosses the x-axis is identified and its z-component is 
calculated. The polytope is then rotated about the x-axis and another edge is identified. 
This is repeated until the z-component of the identified edge changes sign. The last two 
edges identified (one “behind” the x-axis and one “in front” of it) are candidates for 
defining the desired facet The facet is tested according to the conditions in ( 1.16) and 
(117). If the facet is the correct one, the problem is solved. If not, the direction of 
rotation is reversed and a smaller rotation angle is used. This process is repeated until 
the correct facet is identified. This method of bisecting edge searching results in 
computational requirements on par with other control allocation methods currently in 
use. The number of computations increase roughly linearly with the number of controls 
(as in other methods), whereas the sequential facet search computations increase roughly 
quadratically. 76 

1.3.4.2.3 Spherical Coordinate Transformations 

Peterson and Bodson present two options for finding the correct facet based on a 

spherical coordinate transformation 74 The first option requires more on-line 
calculations but requires less memory. The second option eliminates virtually all on-line 
computations but requires a significant amount of memory. The spherical coordinate 
transformation effectively turns the AMS into a 2-dimensional system, where special 
techniques can be used to accelerate the search. Each option is comprised of off-line and 
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on-line computations. In Option 1, ranges are computed for the coordinates of the facets 
(in spherical coordinates) off-line. These ranges define search boxes in which the facets 
are located. These boxes are used on-line to quickly assess whether the desired moment 
is likely to lie within a given facet. If the test is successful, the facet is found quickly 
and the solution is obtained Option 2 makes use of a lookup table. The azimuth and 
elevation angles are quantized and facet numbers are associated with a given pair of 
spherical coordinates. The creation of the table is performed off-line The size of the 
table can easily reach 1*10 6 elements, depending on the fineness of the quantization. 
The on-line computations to find the correct facet then consist of converting the desired 
moment into spherical coordinates followed by a simple table lookup. This guarantees a 
known fixed amount of time to locate the desired facet. 74 Both methods offer 
considerable computational improvements over the sequential search of facets based on 
Durham’s 3D tests. An interesting note, however, is that although the spherical 
coordinate table lookup method provides very fast on-line performance, it is completely 
unsuitable for a reconfigurable control setting If a control effector fails or is degraded, 
the entire lookup table would need to be recalculated Since this is a large off-line task, 
it is not a reasonable candidate for used in an adaptive environment. 

1.3.4.3 Linear Programming 

Page and Steinberg 4 provide an excellent summary of the Linear Programming 
approach. This summary is duplicated below. The Linear Programming method has two 
steps. First, the following linear program is solved, 
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T 

min l ~W^u s 
u 


subject to 


u s 


0 

- u 


“ w max 

u 

> 

w min 

- B U 4- u s 


~m d 

Bw *f u s 


m d 


(1 18) 


where u s is a vector of slack variables, and W A is a positive weighting vector If J > 0, 
the desired moment is unattainable and the solution (denoted u) is that which minimizes 
the weighted 1-norm distance between B u and m A . If the solution to equation (118) 
yields 1=0, then the desired moment is attainable and a second linear program is solved: 

min J = ivj u s 

U 



(1.19) 


13// = m ^ 


where u s is again a vector of slack variables, is a positive weighting vector, and w pre f 
is a vector of control preferences. In this case, the weighted 1-norm distance between 
the control vector solution (w*) and the control preference vector (w pr ef) is minimized. 

4 18 77 79 

There are a number of examples of works which use this method. ’ ’ 
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1.3.4.4 Quadratic Programming 

The Quadratic Programming method is similar to the Linear Programming 
method, except it has a quadratic cost function. Again, it is a two step process. 4 First, 
the equation 

I T 

min J = — (B u - w c j) Wj (Bw - m^) 

u 2 

( I 20) 

subject to u mm < u < « max 


is solved, where W d is a positive definite symmetric weighting matrix. If J > 0, the 
desired moment is unattainable and the solution (denoted u ) is that which minimizes the 
weighted 2-norm distance between Bw* and m<\ If the solution of equation ( 1 .20 ) yields 
J = 0, the desired moment is attainable and a second quadratic program is solved. 


1 7 ' 

nun ] — — (u — Wp re f ) (M ~ ^pref ) 

u 2 

subject to u mm <u< u max 

Bu = 


( 1 . 21 ) 


where W u is a positive definite symmetric weighting matrix. In this case, the weighted 
2-norm distance between the control vector solution («*) and the control preference 
vector (wpref) is minimized. There are a number of works which use this 

method. 4 - 34 ' 66 - 67 ’ 70 


1.3.4.5 Pseudo-Inverse 

There are several variants of the Pseudo-Inverse. In each, the basic idea is to 
minimize the same quadratic cost function that the Quadratic Programming solved 4 
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Note that this is posed with no limits on u. With - oo < u < <x> , the solution to ( 1 .22 ) is 

"=Vef +'V l 7 1 B T (BW,7 l B T )- , (m d -B Wpref ) ( 1 23 ) 

This solution is a biased weighted pseudo-inverse, for which it gets its name. For the 
special case of identity weighting ( W u = I ), this is referred to as the Pseudo-Inverse 

method. With constraints on the control effectors, equation (1.23) may produce an 
inadmissible solution. There are two common ways of dealing with the saturated 
actuators. First, only the individual commands that have violated the given constraints 
are scaled down so each is at its corresponding limit. Note that this individual clipping 
does not preserve the direction of desired moment. The second method works to 
preserve the direction of the resulting moment by scaling all commands by a single 
factor until no constraints are violated by any effector. This method is called the 
Direction Preserving Weighted Pseudo-Inverse method. There are a number of works 
which use some variant of the pseudo-inverse allocation. 4,5 ’ 19,21 ' 24 ’ 25 ’ 35 ' 76 ’ 80 ' 83 

I.3.4.6 Cascaded Generalized Inverse 

This approach is similar to the Weighted Pseudo-Inverse except in the way it 

handles violations of effector constraints If an actuator command reaches its limit, it is 

set to its limit, its effect is subtracted from the desired moment, and it is removed from 

the problem— which is then solved again. 4 This cascading process it repeated until 1) no 

new controls saturate (the desired moment is attainable); 2) all remaining controls 
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saturate (the desired moment is unattainable); or 3) fewer controls remain than the 
number of desired moments. When (3) occurs, the left pseudo-inverse solution is used 

u = (B T B) -1 B T m d (124) 

There are several examples of works which use this method. 4 ’ 76 ’ 84 

1.3.4.7 Daisy Chaining 

In the Daisy Chaining method, the available controls are separated into two or 
more groups. In order to meet the demand of the desired moment, one group is activated 
while the others remain constant. If any elements of the operating group reach their 
limits, that group of controls are held at their last position and the next group of controls 
are activated. 69 At first glance, this approach sounds similar to the Cascaded General 
Inverse, but it is different in that entire groups of controls may remain completely 
unused when daisy chained As an example, consider two groups of controls, three 
aerodynamic controls, u\\ and three thrust vectoring controls, uj . The control 
effectiveness matrix is then partitioned into two 3x3 matrices, Bi and B 2 . 

B«=[B] B 2 ) 1 = B]«] + B 2 W 2 (1.25) 

_ u 2_ 

Since in this application, B[ and B 2 are square and assumed to be invertible by design, 
their inverses are unique. For a given desired moment, the aerodynamic controls will be 

used first until the point of saturation. That is, while w i min - u \ ^ u l max , 

u\ = B\^ uj = 0 (1.26) 
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If any of the aerodynamic controls reach their limit, the group of aerodynamic controls 
are held at that position ( w l(sat) ), and the thrust vectoring controls are brought on-line. 

“1 = u l(sat) “2 = B 2 ! ( m d - B 1 «l(sat>) ( 1 27 ) 

A notable problem arises in the Daisy Chaining method when control deflection rates are 
considered “Cooperative control efforts are those in which all available controls are 
simultaneously varied to meet a time-varying demand. The total rate of change of the 
moment produced is a linear combination of the individual control rates. For a given 
rate of change of the required moment, in magnitude and/or direction, a cooperative 
effort among all available controls will require lower individual control rates than will a 
noncooperative effort. Daisy chaining is a noncooperative allocation scheme and 
potentially will command unattainable deflection rates that would not be commanded by 
cooperative control allocation methods .” 69 There are several examples of works which 
use this method . 09,85 ' 87 


1. 3.4.8 Discrete Time Methods 

As introduced above, all the methods account for position limit constraints on the 
control effectors. They do not, however, account for rate limits. Rate limits are 
typically incorporated within a discrete time framework by considering how far an 
effector can travel within a single time step . 4 Therefore, the control effectors are 
constrained by: 


where 


Au. 


< Au < Au mar 

max 
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Au m,„ = max 

Au m ax =ITlin 

and At = - t k ; Au = u k +i - u k ; u mm and u max are effector position limits; and 

u max an£ l ^max are ff e negative and positive effector rate limits. Each of the control 

allocation methods discussed above can be modified to include these discrete time 
effector rate constraints. 

1.3.4.9 Control Allocation Performance Measures 

Much of the early comparisons between the different control allocation methods 

involve an evaluation of the ability of the allocation method to reach 100% of the 
volume of the Attainable Moment Subset. Bordignon and Durham 87 offer methods to 
compute the volume of the AMS reachable by different allocation schemes They show 
that, while the Direct Allocation method can access 100% of the AMS, the pseudo- 
inverse methods can only access 13-42% of the AMS, and daisy chaining only reaches 
about 22%. They conclude that the Direct Allocation method must be superior since it 
can achieve moments the other methods can not. This seems to be a reasonable 
conclusion; however, recent studies show otherwise. Page and Steinberg 5 ran several 
different control law designs with four different control allocation schemes in order to 
investigate potential adverse interactions between the control laws and the control 
allocation algorithms. They conclude that “the choice of control allocation technique 
can have a dramatic impact on system performance.” They show by example that 
unexpected interactions with the Direct Allocation method (which produces admissible 
solutions for 100% of the AMS) occur in some cases which drive the effectors to 


K~ At > ( U mm- U k) 
U Lx At > ( U max- U k) 
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saturation and greatly degrade the system performance. For these same cases, the 
Weighted Pseudo-Inverse technique (which produces admissible solutions for only 10% 
of the attainable moment subset), did not cause saturation, and performed very well 
One conclusion of their work is that it is important to integrate the control allocation 
effort with the adaptive control laws because unexpected interactions can occur A 
follow-on study compares the open-loop and closed-loop performance of sixteen 
different control allocation methods using a single dynamic inversion control law 4 It 
shows that because the allocation approach is generally designed to meet an 
instantaneous moment demand and not to optimize any closed-loop properties over time, 
the open-loop measures of control allocation performance do not necessarily translate to 
closed-loop performance. The Direct Allocation method again did not perform as well 
as might be expected and the Discrete Time Weighted Pseudo-Inverse method provided 
the best closed-loop tracking performance for the application investigated 
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1.4 Literature Summary 

Section 1.3 has attempted to summarize the varied approaches to reconfigurable 
control system design that are described in the literature. It is obvious that the 
requirements of failure detection, failure isolation, system identification, and control law 
reconfiguration present significant challenges to the control engineer and to the practical 
implementation of a reconfigurable flight control system. It is particularly evident that 
most, if not all, of the approaches reviewed require varying amounts of time to reach a 
“reconfigured” state. With the advent of modern combat aircraft with highly unstable 
unaugmented dynamics, this reconfiguration time can become a critical issue in the 
viability of the design approach Of all the works reviewed, those which utilize sliding 
mode control appear to have the most promising results. In these works, the tracking 
performance after failure is impressive; there is no reconfiguration time at all, and no 
parameter identification is required. For these reasons, sliding mode methods are 
selected for investigation for their applicability to the design and implementation of a 
practical reconfigurable control system. 
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Chapter 2 

Sliding Mode Control 

2.1 Introduction to Sliding Mode Control 

Of all the reconfigurable control schemes appearing in the current literature, those 

based on the concept of sliding modes appear to be the most promising. While results 
from many of the other reconfigurable control methods are very good (the aircraft 
remains stable, and tracking performance, while degraded, ranges from acceptable to 
very good), the recent works utilizing sliding mode control (SMC) show no performance 
degradation at all 63,64,66 This is because these controllers are invariant to matched 
uncertainty/disturbances. If the controller is invariant to certain system parameter 
changes, there is no need to perform system failure detection or parameter identification. 
Since parameter identification is the main bottleneck in control reconfiguration, 
employing a controller that does not require it provides a large advantage Not only this, 
but its variable structure allows SMC to adapt to parameter disturbances 
“instantaneously.” An ideal sliding mode controller with no parasitic dynamics or 
actuator limits easily handles noise, parameter changes, and unmodeled nonlinearities 
with absolutely no degradation in tracking performance. However, since real systems 
always have parasitic dynamics and actuator limits, the actual implementation of an 
SMC design becomes challenging. The general concepts of sliding mode control, 
including basic theory, properties, design techniques, implementation difficulties, and 
several examples are offered this chapter. 


Chapter 2: Sliding Mode Control 


2. 1 Introduction to Sliding Mode Control 


2.1.1 Sliding Mode Control History 

The basic concepts of sliding mode control first appeared in the Russian literature 
in the early 1930’s. Kulebakin (1932) used essentially a sliding mode controller in the 
context of voltage control for a DC generator on an aircraft He called it “vibration 
control” of the voltage signal. Nikolski (1934) proposed sliding mode relays for 
controlling the course of a ship He actually used many of the terms currently in use 
today, including phase plane, switching line, and sliding mode. 88 Emelyanov and 
Barbashin continued the pioneering work in Russia in the early 1960’s, ft wasn’t until 
the mid I970’s that the ideas of sliding modes appeared outside Russia when a text by 
Itkis 89 and a survey paper by Utkin 90 were published in English. 91 Vladim Utkin is one 
of the key figures who championed the concepts of sliding modes through the 1970’s 
and 1980’s, and continues to publish numerous works in the area. The I980’s saw a 
large increase in interest in SMC including several often-cited survey works 92-94 
Several newer survey works 95 ' 97 are also regularly cited Utkin has several texts on the 
subject. ’ Slotine includes a chapter in his text on non-linear control; and Edwards & 
Spurgeon have a recent text 91 devoted entirely to SMC Applications of SMC in the 
literature are too numerous to list By 1993, general application areas included: robotic 
control, motor control, aircraft and spacecraft control, flexible structure control, load 
frequency control of power systems, servomechanisms, pulse-width modulation control, 
guidance, process control, phase-locked loop control, power converters, and remote 
vehicle control. 95 Sliding mode controllers are even being used for controlling the 
convergence rates for neural net learning algorithms. 100 Other interesting applications 
which have appeared in very recent works include, direct robust exact differentiation, 101 
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large-scale, time-delay systems; 102 missile autopilot; 103 multiple unmanned air vehicles 
in close-formation flight; 104,105 and of course, reconfigurable flight control. 31,63 ' 66 ’ 106 

2.1.2 Concept of the Sliding Mode 

Sliding mode controllers are a subset of a class of controllers known as Variable 
Structure Controls (VSC)— although the terms are often used interchangeably. VSC 
systems are characterized by a control structure that changes according to some 
predefined rule which is a function of the states of the system. In order to illustrate this 
concept, consider the double integrator in equation (2.1) with a feedback control law in 
equation ( 2.2 ). See Edwards & Spurgeon 91 for much of the following development 

y (0 = u(t) (2 1) 

u(t) = -ky(t) (2.2) 

This results in a pure undamped harmonic motion as shown in Figure 2-1 with 
y(0) = 0, y(0) = 1 , and k = 4. The phase plane portrait of this system is shown in 
Figure 2-2 Obviously, this control law would not be appropriate for this system since 
the state variables do not move toward the origin for any value of gain chosen The gain 
only changes the orientation and eccentricity of the ellipse in the phase plane. 



Time (sec) 

Figure 2-1: Time History of Oscillator Figure 2-2: Phase Plane Plot of 

Oscillator 
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Consider, instead, a Variable Structure Control law as shown in (23). 

u(t) = 


[ - k, y(t) if yy < 0 
[-k 2 y(t) otherwise 


(2.3) 


This causes the control law to change depending on the quadrant of the phase plane in 
which the state variables are moving. When in quadrants 1 & 3, the gain is kg when in 
quadrants 2 & 4, the gain is k 2 . The resulting motion, for y(0) = 0, y(0) = 1 , ki = 0.5, 
k 2 = 4, is shown in Figure 2-3 and Figure 2-4. 




Time (sec) 


Figure 2-3: Time History of VSC 
System 


Figure 2-4: Phase Plane Plot of VSC 
System 


This simple example shows that by introducing a rule for switching between two control 
structures, a stable closed-loop system can be obtained. This occurs even though neither 
of the two control structures individually are able to provide convergence of the state 
variables to zero. 

Next, instead of using the phase plane quadrants for the switching rule, consider 
using a switching function: 

a(y, y) = m y + y (2.4) 

Where m is a positive scalar. 
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Then, the control law is defined in terms of the switching function as follows 


f-1 ifo(y,y)> 

I 1 ifo(y,y)<0 


The resulting motion (with m = 1) for this system is shown in Figure 2-5 and Figure 2-6 



01 23456789 10 


Time (sec) 


Figure 2-5: Time History for SMC System Figure 2-6: Phase Plane 

Plot for SMC System 


This is called a Sliding Mode controller because, after reaching the switching function, 
the state trajectory follows (or “slides along”) the switching line to the equilibrium point. 
This is the goal of any sliding mode controller: drive the states to the switching surface 
in a finite time and constrain the states to remain on the surface a = 0 for all subsequent 
time. As illustrated in Figure 2-7, the time period where the states are moving toward 
the sliding surface is known as the reaching phase, and fhe phase where the states follow 
the surface a = 0 is called the ideal sliding mode 
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Figure 2-7: Phases of Sliding Mode Control 

“A 

This is also clearly illustrated in the time history of the switching function from the 
previous example. In this case, the system reaches the sliding mode in about t = 0 8 sec. 



Time (s) 

Figure 2-8: Time History of Sliding Surface 

In order to maintain the motion on the sliding mode, the control effort requires 
infinite frequency switching because the control law is undefined on the surface a = 0. 
If this infinite frequency switching were possible, the motion would be constrained to 
the sliding mode, and the dynamics of the closed-loop system would be that of the 
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switching surface a = 0. For this example, setting a = 0 and rearranging gives the 
motion for the closed-loop system after reaching the sliding mode: 

y(t) = -my(t) (2.6) 

This is a first-order decay which depends only on the design parameter m. To further 
illustrate the fact that the dynamics depend only on m , consider adding a non-linear term 
to the double integrator system 

y(t) = -a sin(y) + u(t) (2.7) 

This new system looks like a normalized frictionless pendulum. Using the exact same 
sliding mode controller as before and including the non-linear plant term, the resulting 
phase plane plot looks like. 



Figure 2-9: SMC Phase Plane Plot for System with Non-linear Term 

The previous solution of the double integrator (a=0) is shown for comparison purposes. 
Note that there is a difference between the two cases during the reaching phase; 
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however, once on the sliding mode, both systems exhibit identical closed-loop response 
The non-linear term (a sin(y)), which can be interpreted as a structured modeling error or 
disturbance, has been completely rejected. Although not proven here, there are four 
basic notions about SMC that can be observed 95 

• The line (or hypersurface) that describes a = 0 defines the transient 
response of the system during the sliding mode. 

• During the sliding mode, the trajectory dynamics are of a lower order than 
the original model. 

• While on the sliding mode, the system dynamics are solely governed by the 
parameters that describe the line a = 0. 

• The trajectory of the sliding mode is one that is not inherent in either of the 
two control structures alone 

A formal statement of the sliding mode control problem is now developed. 

2.1.3 Sliding Mode Control Problem Statement 

Consider the uncertain system with m inputs and n states given by: 

i(t) = A(x,t) + B(x,t)w(t) + /(t,x,u) (28) 

where A e lH nXn and B e 9f n * m ; B is full rank, and 1 < m < n The function 

/ : 9? xfR" x9i m i— > 9?" represents the parameter uncertainty or nonlinearities present in 
the system and is assumed to be unknown but bounded by some known functions of the 
state. 
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The objective is to define: 

• m switching functions, represented in vector form as o(x ) with the desired 
state trajectories 

• a variable structure control 

«0,t) = />sgn(a) (2.9) 

such that any state outside the switching surface is driven to the surface in finite 
time and remains on this surface for all subsequent time 

2.1.4 Equivalent Control Concept 

Since the control action in the ideal sliding mode is discontinuous, the resulting 
differential equation can not be analyzed with traditional methods. The approach which 
is almost always cited in the literature is the classic method proposed by Filippov 107 for 
differential equations with discontinuous right hand sides. While Filippov’s constaiction 
provides a theoretical basis for solving the discontinuous problem, a more useful 
approach is Utkin’s concept of the equivalent control 90 Loosely speaking, the 
equivalent control is the continuous control action needed to maintain the ideal sliding 
motion One of its key strengths is that it allows a continuous analysis and provides a 
means to prove certain properties of the sliding mode. 

Assume the switching surface is reached at some time t s and an ideal sliding 
motion is established. By definition, this means o(t) = 0 for all t > / y , which implies that 
c(/) = 0 for all 1 > t s . Since, as will be shown later, the order of the switching surface is 

designed to be exactly one degree less than the relative order of the plant, the input u 
will appear in the derivative of the switching function Setting this equal to zero and 
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solving for u results in the so-called equivalent control. As a simple example, recall the 
problem of the normalized pendulum: 
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y(t) = -asin(y) + u(t) 

The switching line is defined as: 

a(t) = my(t) + y(0 

Taking its derivative and substituting equation ( 2.10 ) gives: 

6(t) = m y(t) + y(t) 

= m y(t) - a sin( y) + u(t) 

Finally, setting d = 0 and solving for u gives u eq : 

%(0 = a sin(y)-m y(t) (t>t s ) 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


In general, assuming a linear, time invariant state space model and a switching surface 

o(x) = Sr (214) 

where S e 9f nxn [ s f u |] ran |( anc ] chosen by design to meet dynamic closed-loop 
requirements, it is easy to show that the equivalent control is uniquely defined by: 91 


u eq (t) = -(SB)-*S A x(t) 


(2.15) 


If the control action is purely discontinuous as in equation ( 2.9 ), Ue*, can be thought of 
as the averaged control signal applied. It can also be shown to be equivalent to the low 
frequency component of the discontinuous signal by passing that signal through a low 
pass filter 91 The equivalent control, however, is not actually applied in practice — at 
least, not alone If the plant model is exactly correct (i.e. there are no unmodeled 
dynamics or disturbances), u eq will maintain the sliding mode assuming the switching 
manifold has been reached (by definition). This will be demonstrated in an example 
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later. However, the discontinuous control function is required to maintain the sliding 
mode in the presence of model uncertainties and disturbances It is common to see the 
equivalent control included in the control law, as in (2.16). 

«(x,t) = M eq + psgn(<j) (2. 16) 

This has the advantage of reducing the activity of the discontinuous portion of the 
control. Even so, it is not usually implemented because it adds complexity to the 
controller. It is common to see it included in initial problem developments and then set 
equal to zero for the implementation. One of its main strengths is its use in proofs 
concerning sliding mode properties and stability. 

2.1.5 Properties of the Sliding Mode 

Unfortunately, very little can be found in the literature about the properties of 
transient response while in the reaching phase. Hung et al 95 make only a brief comment 
about it in their survey paper. Utkin 98 and Edwards & Spurgeon 91 also discuss it briefly 
but do not offer a full treatment. Choi, et a/. 108 offer some of the best treatment of the 
reaching phase, but is too involved to include here Since the reaching phase is typically 
very short, and since the states are moving toward a stable manifold, the reaching phase 
is usually neglected and its properties are of little interest. 

For the following analysis, assume a linear, time invariant state space plant and 
assume the switching surface is defined by equation ( 2. 14 ) with the square matrix SB 
nonsingular. The ideal sliding motion is found by substituting the equivalent control, 
equation ( 2.15 ), into the state space equatioft, which yields a free motion independent 
of the control action: 

x(t) = (l n - B (S B) -i S ) A x(t) Vt>l s and Sx(t s ) = 0 (2.17) 
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SMC Property 1: The sliding motion given in equation ( 2.17 ) is of reduced order and 
the eigenvectors associated with any nonzero eigenvalues of the system matrix 

A eq = (l n - B (SB) -1 S ) A (2.18) 

belong to the null space of the matrix S. 

Proof: See Edwards & Spurgeon. 91 

It follows that A eq can have, at most, n-m nonzero eigenvalues. 

Next, the invariance property is developed. Let the disturbance function / (from 
equation ( 2 8 )) be J{ t,x,u) = D ^(t,x) where the matrix D e is known and the 
function £, : 9f+ x 9?" — > is unknown. This function can be interpreted as an 
exogenous disturbance acting on the system or uncertainty in the system matrices A and 
B. The uncertain linear system can then by expressed as. 

i(t) = Ax(t) + Bw(t)+ D£,(t,x) (2 19) 

Any uncertainty in form given above with /?(D) c: /?(B) is called matched uncertainty. 
Any uncertainty which does not lie with the range space of B is called unmatched 
uncertainty 91 

Using the term in equation ( 2. 17 ), define for convenience: 

P s = (l n ~ B (SB ) 1 s) (2.20) 

Edwards & Spurgeon 91 call this a projection operator and show that is satisfies two 
important properties: 

SP s = 0 and P s B = 0 (221) 

With these definitions in place, the invariance property can be stated 
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SMC Property 2: The ideal sliding motion is completely insensitive to the uncertain 
function £,(t,x) in equation ( 2. 1 9 ) if R(D) cr R(B). 

Proof (from 91 ): The equivalent control for (219) assuming( 214) is found as before 
and is given by 

«eq(t) = -(SB) _1 (S A x(t) + S D <f(t)) for t > t s ( 2.22 ) 

Substituting this into (2.19) gives 

x(t) = P s Ax(t) + P s D<f(t, x) Vt > t s and Sx(t s ) = 0 (2 23 ) 

Now, if /?(D) c ^(B), then D can be decomposed into D = BR, where R e is a 
matrix of elementary operations. As a result, it follows that P S D = P S (BR) = (P S B)R = 0 
by the property given in equation ( 2.21 ). The sliding motion equation then reduces to 
x(t) = P s A x(t) V t > t s and S x(t s ) = 0 ( 2 24 ) 

which is independent of the exogenous signal. 

This is a very important result and is the primary reason sliding mode control is 
attractive for a reconfigurable control setting. Note, it says nothing about invariance to 
unmatched uncertainty In fact, SMC is not invariant to unmatched uncertainty; 
however, measures can be taken to ensure a certain level of robustness. 

2,1.6 The Reachability Problem 

So far, nothing has been said about guaranteeing that the system will reach the 
sliding mode or that it will remain on the sliding mode once it gets there. In fact, this is 
actually the primary concern and is the area that receives the most attention when 
developing the control law. Existence of the sliding mode requires stability of the state 
trajectory to the sliding surface a(x) = 0, at least in some neighborhood surrounding the 
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surface. This neighborhood is known as the region of attraction. 93 If the sliding mode 
is globally reachable, the domain of attraction is the entire state space. The most 
common approach for proving stability, at least in the SISO case, is a Lyapunov 
analysis. 

Geometrically, in order for the sliding surface to be attractive, the trajectories of 
a(t) must always be directed toward it Stated more formally, 

lim a < 0 and lim a > 0 in some domain Oc 93" ( 2 25 ) 

This can be stated in a more succinct manner as 

6a <0 (2.26) 

This is called the reachability condition 91 As illustrated in Figure 2-10, if state 
variables are within the domain O, they will be attracted to the switching surface and 
remain on the surface if the reachability condition is met 



Figure 2-10: Phase Plane Plot of Domain of Attraction 


Noting that (in the SISO case) 


1 d 2 
a 

2 dt 


= 6 o , it follows that the positive definite 


function 
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V(a) = ia 2 ( 2 27) 

is a good candidate Lyapunov function. If it can be shown that V(o) < 0 , the system is 
guaranteed to reach the sliding surface asymptotically 91 

A stronger condition which guarantees an ideal sliding motion in a finite time is 
the so-called ^-reachability condition given by 

a a < -r||a] ( 2 28) 


where ri is a small positive constant 


91 


Noting that ~~~° 2 =oa, and integrating 


equation ( 2.28 ) from 0 to 4, it follows that 

|a(t s )|-|a(0)|<^t s (2.29) 

This, then can be solved for the time (t s ) required to reach the sliding surface (a = 0). 


t < 1 Q(0) 1 

^ X — 


( 2 30) 


Typically, what is done in a design analysis is the candidate Lyapunov function 
V(ct) = — a 2 is assumed. The control law has parameters (for example, p in equation 


( 2 9 )) which need to be determined System uncertainties are parameterized and 
replaced into the state equation, which then appear in the Lyapunov function. Then the 
control parameter(s) are found in terms of the given uncertainty bounds such that the 
sliding mode is globally attractive. 


2.1.7 Design Approaches 

There are many design approaches in the literature — too many, in fact, to present 
here. In theory, there are an infinite variety of control strategies to achieve the sliding 
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mode There are diagonalization methods, hierarchical methods, relays with constant 
gains, relays with state dependent gains, linear feedback with switched gains, linear 
continuous feedback, and univector non-linearity with scale factor. 93 All approaches 
consist of two basic steps: 1) design of the sliding surface(s), and 2) design of the 

control law to achieve the sliding mode. Two common approaches which are the most 
intuitive and which will be employed in this work are given below 

One of the most commonly cited approaches is some variation of the so-called 
regular form approach. Edwards & Spurgeon 91 present this most succinctly, so their 
development is essentially duplicated here 

Consider the nominal linear model of an uncertain system, given by 

x(t) = Ax(t) + Bw(t) (2.31 ) 

where rank(B) = m and (A,B) is a controllable pair. Define an associated switching 
function 


o(t) = S x(t) ( 2.32 ) 

This system can be transformed into regular form via a change of coordinates defined by 
an orthogonal matrix T r such that: 


z(t) = T r x(t) ( 2 33 ) 

where T r is found by a QR decomposition of the input distribution matrix, that is, 


Then, defining 


and 


T r B = 



TAT. = 


21 



ST r T =[S, Sj 


(2.34) 

( 2 35 ) 
(2.36) 
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The system can be expressed in the well-known regular form: 


zi(t) = An zi(t) + A12Z2O) 

z 2 (t) = A 2 i zi(t) + A 22 zi(t) + B 2 w(t) 


(237 ) 


and 


o(t) = S,z,(t) + S 2 Z 2 (t) ( 2 38 ) 

Note that in this form, it is easy to identify matched and unmatched uncertainty. 
Anything appearing in the input line is “matched.” Therefore, the closed-loop response 
is insensitive to variations in A 2 i, A 22 , and B 2 . 

During the sliding motion, the switching function must be identically zero, so 

S, zi(t) + S 2 z 2 (t) = 0 (2 39) 

It can be shown that S 2 is nonsingular, so z 2 can be solved for on the sliding mode. 


= 


-S 2 -'S 1 z 1 (t) 
: -M Zj(t) 


(2 40) 


where 1VI e 9?" 


x(n-m) 


is defined as 


M =S 2 ~ l S, 


(2.41) 


The sliding mode is then governed by equations ( 2.42 ) and ( 2.43 ) 


^l(0 = A ll- z l( t ) + A L2 z 2( t ) (2.42) 

Z 2 (t) = -Mzi(t) (2.43) 

This is an ( n - m) ih order system in which z 2 acts in the role of a linear full-state 
feedback control signal. Closing this loop gives the free motion of the system: 

i 'l(0 = ( A ll + A 12 M ) z l( t ) (2.44) 
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Notice, once again, that the terms which contain the matched uncertainty do not 
appear — thus demonstrating again SMC’s invariance to matched uncertainty. Notice, 
too, that the problem of selecting the sliding surface has turned into a standard linear 
full -state feedback problem. In order to ensure the reaching condition, it must be 
ensured that Z| is asymptotically stable. Further, since all the unmatched uncertainty 
appears here, it is desired that this system be made as robust as possible. Edwards & 
Spurgeon detail two standard methods: robust eigenstructure assignment and quadratic 
minimization (the standard LQR problem). 


Another method for SMC design is based on the ideas of feedback linearization. 
Fernandez & Hedrick 94 and Slotine" approach the problem this way, as does Yuri 
Shtessel, a very prolific SMC proponent in the recent literature. 63 ' 66 ’ 106,109 ' 120 A major 
assumption with this approach is that the system must be square — equal number of 
inputs and outputs, and it must be feedback linearizable. If the system is square and 
feedback linearizable, it is possible to decouple the outputs with the given inputs. This 
turns the design into m simple SISO-like designs. 

Consider a non-linear square MIMO system 


x = f{x) + G(x) u 

y = Kx) 


( 2 . 45 ) 


where x e 9l n , y e 9t m , u e 9i m . Assume the functions J(x), h(x) and columns g,(x) 
Vi = i,m of the matrix G(x) e are smooth vector fields. Further, assume the 

system is completely linearizable in a reasonable domain xef The control system will 

be designed to track a real-time reference profile, TrO) . 
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This system can be transformed to a normal form 


99 


E(x) = 




L, r > h,(x) 

y 2 (r ’> 

= 

L f f2 h 2 (x) 


1 

L f r - h m (x) 


+ E(x) u. 


L g) (Lf _ 

'hi) 

L g2 (Lh-’ h ,) •• 



LgiOf 

'h 2 ) 

L ga (L r r , h 2 ) ; 

■ hJ L 'r' h 2) 

|e(x)| ^ o vx € r 

L g,( Lr r 

'hm) 

L g2 (L r f-'hm) 

• L g m ( Lr r' h m)_ 



(2.46) 


Where L r ; h , and L g (L'f *h ,) Vi = l,m are corresponding Lie derivatives 120 

Next design m independent sliding surfaces (note, these have orders exactly one 
less than the relative order for the corresponding state variable): 

a, = e^ _l) + c, f _ 2 e< r, ' 2) + . . . + c, ^ + c, 0 e, Vi = 1, m ( 2.47 ) 

d j e 

where e, = y r , (t) — y, (t), e[ j) = — L . The coefficients c Vi = l,m and Vj = 0, r, -2, 

dt J 

are design parameters which are chosen to achieve the desired eigenvalue placement of 
the decoupled differential equations of the output variables. It is also common to include 

an integral term, cje^x, in the sliding equations to account for potential steady state 

error which can occur when utilizing a sliding mode boundary layer (to be discussed 
later) 

The control law which can be used is 

u , = u «,. +P 1 sgn(a i ) ( 248 ) 

In order to prove system stability, assume the candidate Lyapunov function 
V,= — a 2 , take the derivative of the sliding functions, 0 ;, and solve for p, which 


provides global attractiveness to the sliding surface in finite time. 
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2.1.8 Implementation Issues 

While sliding mode controllers are very attractive from a robustness standpoint, 
there are some implementation issues which must be addressed. In fact, much of what 
appears in the literature for sliding mode control is related to the issue of 
chatter 59,109 ’ 112 ’ 113116 ’ 119 ' 128 Without a solution to the chatter problem, SMC would be 
relegated to only a very few applications where the control action can attain the high 
frequency switching required, such as control of electric motors and power converters. 

2.1. 8.1 Chatter 

Chatter is defined by many authors as the high frequency action of the actuators. 
Young & Utkin 97 take a more purist stand on the definition of chatter. The high 
frequency action of the control — by their definition — is not chatter, but rather the 
expected result of the SMC control law (recall the control output is undefined on the 
surface a = 0). An infinite frequency switching is actually required while on the sliding 
mode. This results in a smooth ideal sliding along the sliding manifold with smooth 
state trajectories. However, due to parasitic dynamics in the real system, such as 
actuator dynamics and time delays, the switching can not take place at the required 
frequency, and the state trajectories chatter along the sliding manifold After working 
with some example models, it is found that not only do actuator dynamics result in 
chatter, but they can drive the system unstable. At first, this may seem to contradict the 
invariance property (since the actuator dynamics are part of the control action). 
However, recall that the invariance property is only valid while on the sliding mode. 
Theoretically, this requires an infinite frequency switching. Actuator dynamics prevent 
the controller from maintaining the sliding mode. There are several methods proposed 
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in the literature to deal with the chatter problem. Two of the most common approaches 
are briefly introduced. 

The simplest and clearly most popular approach is the so-called boundary layer. 
The basic idea is to replace the discontinuous signum function with an arbitrarily close 
approximation The signum function looks like 


v = sgn(a) = 


if a > 0 
if a < 0 


Several such approximations include. 


v = sat 




V 8 J 


if a > e 
otherwise 


0 


saturation element 


v = < 


5 q-i 


a 


a 


0 


if a > e 


if 0 < a < 8 


a = 0 


power law 


v = 


cr -f-£ 



differentiable approximation 
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Some research has been done to determine the performance of these different 
approximations to the signum function. Results show that there are no significant 
differences between them in the application to SMC. 91 The result of using one of these 
continuous approximations is the states become attracted to a small boundary layer 
surrounding the switching surface, rather that to the switching surface itself. 
Convergence proofs show stability to the boundary layer and the states remain arbitrarily 
close to the switching surface. Figure 2-11 demonstrates the effect of a saturation 
element boundary layer on the state trajectories. 



Ideal SMC Boundary Layer SMC 

Figure 2-11: State Trajectories for Ideal Sliding Mode and Sliding Mode with 

Boundary Layer 


Of course, the ideal sliding mode is lost and the resulting motion is often referred to in 
the literature as pseudo-slidmg. While invariance is also lost, the system still retains 
much of its robustness. The boundary layer can be made arbitrarily small to approach 
ideal sliding Typically in the design process, the boundary layer thickness is tuned to 
achieve an acceptable balance between maximum performance/robustness and reduced 

chatter. This approach has the advantage of being very simple to implement and is the 
Chapter 2: Sliding Mode Control 2.1 Introduction to Sliding Mode Control 
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choice of many authors. 31,63,66 ’ I03 ‘ 105,121 ' 129 ' 131 One disadvantage is the resulting control 
becomes essentially a high gain controller — something that is typically avoided in 
control design. 

Another method to reduce chatter in SMC systems is an approach called second- 
order sliding modes, or 2-SMC. 59 ’ 101,113 ’ 122 ’ 127 ’ 128 The basic idea of this approach is to 
use the first derivative of the control as a pseudo-control during the design Therefore 
all the high frequency switching will occur in u , and when this signal is integrated, the 
actual control will be continuous. There are also higher order sliding mode controllers 
which have been proposed 128 This technique works well, but a key disadvantage is that 
derivatives of the states appear in the switching function — which means these 
derivatives must be measured or calculated 

2.1. 8.2 Actuator Limits 

Little attention is given in the literature concerning the effects of actuator limits in 
SMC. Unfortunately, actuator limits can have large, undesirable effects on the system. 
The type of actuator limits that are typically associated with degraded control 
performance are rate and position saturation. If a rate limit is encountered, it appears as 
a degradation of actuator bandwidth. The result is added phase lag — a condition which 
almost always results in chattering or loss of the sliding mode. If a position limit is 
encountered, global stability can not be guaranteed. Probably one reason position limits 
are not discussed is the very nature of the SMC control law. Consider, again, equation 
( 2.9 ) : w(x,t) = p sgn(o). The maximum position this control will see is ±p. If p is less 
than the maximum allowable control effort, there will not be any position saturation. 
Since p is a design parameter, it should be chosen to be less than the maximum 
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available. However, since p is typically a function of the maximum expected 
uncertainty or disturbance, it is possible to desire more control authority than what is 
available. This results in a reduced region of attraction, and if the states are outside 
of this region, stability can be lost. This is demonstrated in one of the examples in 
Chapter 3 Madani-Esfahani et al . 132 investigate the problem of estimating regions of 
asymptotic stability (RAS) in variable structure systems with hard bounds on the control 
action. The paper is highly theoretical and shows that, using their approach, the RAS 
cannot be found analytically using a single Lyapunov function. 

The problem of finite bandwidth actuators with rate and position limits is a major 
implementation issue and comprises the major of the effort in this research. Chapter 3 is 
develops the problem further and outlines some potential solutions. 
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2.2 Application Example: Inverted Pendulum on a Translating Cart 

In order to demonstrate the basic concepts of sliding mode control, the classic 

problem of an inverted pendulum on a translating cart is considered: 

System parameters. 



Cart mass 

= M 

Pendulum mass 

= m 

Pendulum length 

= L 

State variables 

Cart position 

= X 

Pendulum angle 

= 0 

Control inputs: 

Horizontal Force 

= u 

Pendulum Torque 

= T 


If the horizontal friction force and the frictional torque on the pendulum are included, the 
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non-linear equations of motion are given by. 

3 0)0-m L0 2 sin 0 = u 

(2.49) 

Linearizing about the equilibrium point, 0=0, the linearized state space equations are: 


(M + m) x + F x x + (m L cos 0)0 -m L0 2 sin 0 = u 
J 0 + F e 0 — m Lg sin 0 + (m Lcos0)x = x 
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(2.50) 
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Parameter 


Nominal Value 


Cart Mass, M 
Pendulum Mass, m 
Rod Length, L 
Linear Friction Coef, F ? 
Angular Friction Coef, Fp 


3.0 kg 
0.50 kg 
0 4 m 
6 0 kg/s 
0.005 kg m 2 


Value after 
“Failure” 

30.0 kg - 
5 0 kg 
0.8 m 

25 76 kg/s 
0.05 kg-m 2 


Initially, consider a SISO system with the force as the only input (x = 0). Then, the state 
space representation of the linearized system is 


Z 1 


'o 

0 

1 

z 2 


0 

0 

0 

z 3 


0 

-1.6345 

-2 

Aj 


0 

28.6037 

5 


0 

z . 


0 

1 

Z 2 

+ 

0 

0.0042 

Z 3 


0.3333 

-0.0729 

_ Z 4 _ 


-0.8333 
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A sliding mode controller is designed using the regular form approach. 

- Design the sliding surface 

- Change coordinates of given system to regular form (see Edwards & Spurgeon 91 , p.66 
for m-file) 

- Perform QR decomposition on input distribution matrix to get T r 

- Obtain A rcg and B reg using T r 

- Obtain matrix sub-blocks in the regular form equations 

0 10 O' 

-0 3714 0 0.8621 0.3448 

T r = 

1 0 9285 0 0.3448 0.1379 

0 0 -0.3714 0.9285 


%\ 


0 

0.3448 

0.1379 

0.9285 " 

X, 


0 



8.4543 

-0.3276 

-0.1310 

0.1179 

x 2 

+ 

0 

X 3 


3.3817 

0.7974 

0.3190 

-0.3528 

x 3 

0 

X4 _ 


27 1650 

46185 

18474 1 

-2.0643 

_X 4 


-0.8975 


- Use linear quadratic cost function to design the switching function matrix coefficients 
(see Edwards & Spurgeon, 91 p.74 for m-file) 

- Transform weighting matrix to regular form coordinates 

The weighing matrices in the original coordinates are the standard weights in an LQG 
problem. For this example, the values chosen are. 


10 0 0 0 
0 10 0 0 
0 0 10 
0 0 0 1 


R = 0.01 


- Compatibly partition weighting matrix with regular form description 
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- Form reduced order system description and associated weighting matrix 

- Solve the LQR problem 

- Transform solution back into original coordinates 

S = [3.1632 10.1219 2.7436 2.1745] (in original 

coordinates) 

Calculate U e(i : U eq = -(S B)" 1 SAx 

U eq =[o 64.3035 9.5233 11.1136] * 

- The final control parameter to design is the gain p. In the following examples, 
the value of p varies depending on the particular case being run Its value will be called 
out in each example. The block diagram of the sliding mode controller just designed is 
given in Figure2-13. 



Chapter 2: Sliding Mode Control 2.2 Application Example: Inverted Pendulum on a Translating Cart 





68 


For the sake of comparison, an LQR controller is designed using the same weights 
as above. The resulting feedback control gains for the LQR controller are. 

K L qr = [-31.6 -179.1 -38.7 -35.8] 

The next series of figures are ones which demonstrate the basic operation of the 
SMC just designed. A simulation with the linear model and no noise is used. These 
demonstrate the effects of equivalent control and the boundary layer. 


Figure 2-14 
Figure 2-15: 
Figure 2-16. 
Figure 2-17. 
Figure 2-18: 


LQR Regulator (Linear Model) 

SMC Regulator (Linear Model) with u 
SMC Regulator (Linear Model) with u 
SMC Regulator (Linear Model) with u 
SMC Regulator (Linear Model) with u 


U eq 

20 sat(a/0.005) 

U eq + 5 sat(o/0.0001) 
U eq + 5 sat(a/0. 1 ) 
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For the sake of comparison, this plot shows the results of a simulation using a standard 
LQR regulator using the same weights given above. As expected, an LQR controller 
does an acceptable job. This same controller is used on the full non-linear system 
simulation and achieves exactly the same results. That plot is not given because it looks 
essentially the same as this one 



Time, (sec) 


Figure 2-14: LQR Regulator (Linear Model) 
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This SMC controller uses only the equivalent control — it has no discontinuous control 
element Note, Ueq will not reach the sliding mode by itself A typical discontinuous 
element is used during the reaching phase (0 < t < 0. 15 sec) and turned off as soon as the 
sliding mode is reached. As expected, U eq maintains the sliding mode for this case 
because the simulation model is exactly the same as the design model, and there is no 
noise. When this same control is applied to the non-linear simulation (or if noise is 
present), stable control is lost. 



Time, (sec) 

Figure 2-15: SMC Regulator (Linear Model) with u = U eq 
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This SMC controller uses no equivalent control. Notice the high activity of the control 
(near infinite frequency) This is because this case uses a fairly small boundary layer. 
As the boundary layer approaches zero (ideal sliding mode), the frequency of the control 
activity becomes infinite. This same controller is used on the full non-linear system 
simulation and achieves exactly the same results. 




01 23456789 10 

Time, (sec) 


Figure 2-16: SMC Regulator (Linear Model) with u = 20 sat(a/0.005) 
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This SMC controller uses both the equivalent control and the discontinuous element. 
Note that the magnitude of the discontinuous part is considerably less than that required 
for the previous example. That is because the continuous equivalent control is doing 
much of the work. The boundary layer is very small in this case to show the effect of 
approaching the ideal sliding mode. This same controller is used on the full non-linear 
system simulation and achieves exactly the same results. 





Time, (sec) 


Figure 2-17: SMC Regulator (Linear Model) with u = U eq + 5 sat(o/0.0001) 
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This SMC controller uses both the equivalent control and the discontinuous element 
This case differs from the previous case in the size of the boundary layer. In this case 
the boundary layer is increased to the point that control appears completely continuous. 
If the boundary layer is increased too much, the performance becomes noticeably 
degraded. This same controller is used on the full non-linear system simulation and 
achieves exactly the same results. 





Figure 2-18: SMC Regulator (Linear Model) with u = U eq + 5 sat(a/0.1) 
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Next the invariance property of the sliding mode controller is tested. If the 
physical properties of the system are varied, the controller proves to be remarkably 
robust. However, it is not much more robust than the LQR controller — which also 
performs very well. This result is to be expected. SMC systems are not invariant to 
unmatched uncertainty. If the system in regular form is examined, it is clear that the 
physical properties all appear in the unmatched terms of A eq . This means, that while the 
SMC will be invariant to variations in the A 21 and A 2 2 terms, it will only be as robust to 
changes in An and A 12 as the system in equation ( 2.44 ) can be made to be — which was 
designed using an LQR approach. The plots for these cases are not presented, but both 
the LQR controller and the SMC controller perform very well with state noise, a non- 
linear system simulation, large state initial conditions, and varying system parameters. 
In general, the SMC controller is more robust, but not by much 

In order to properly test the invariance property, matched disturbances need to be 
input into the system The simplest way to implement this is to inject disturbances 
directly at the control input. A white noise is added directly to the control input 
Additionally, at / = 5 sec, a single step bias with an amplitude of 10 is injected at the 
control input. This might simulate a control “hard over.” 

The results of this simulation for the LQR system and the SMC system are given 
in the next two figures. 

Figure 2-19 LQR Regulator (Non-linear System) with Control Input Disturbances 
Figure 2-20: SMC Regulator (Non-linear System) with Control Input Disturbances 
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This is the LQR controller using the same weights as previously. White Gaussian noise 
is injected at the control input. A disturbance bias of 10 is added to the control input at 
/ = 5 sec. Pendulum angle tracking is noisy. Cart position has a noticeable bias after 
5 sec 



Time, (sec) 


Figure 2-19: LQR Regulator (Non-linear System) with Control Input Disturbances 
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Next, the control torque is included in the system to make the system square. 
Using the idea of feedback linearization outlined above, a sliding mode controller which 
decouples the outputs is designed The resulting design is 

a x = (x c - x) + 10 (x c - x) a 0 = (0 C -0) + lO (0 c -0) Switching surfaces 


u = 50 sat 



x = 1 0 sat 



Control laws 


As before, a LQR controller is also designed for comparison purposes. 


For the following plots, the command signal for x is a sine wave and the initial 
state of x(0) = 0. The command signal for 0 is a cosine and the initial state of 0(0) = 1 . 
All cases are run with the non-linear system 

The first plot is the SMC baseline case with a healthy plant. The next four plots 
show the results when the plant model experiences a “failure” as defined above. 

Figure 2-2 1 . Decoupled SMC Tracking, Baseline Case 

Figure 2-22: LQR Tracking with System Failure 

Figure 2-23. SMC Tracking with System Failure 

Figure 2-24 SMC Tracking with System Failure (Insufficient px & p0) 

Figure 2-25 SMC Tracking with System Failure (Insufficient p0) 
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(a ) 


f G 0 ^ 

This plot shows the SMC controller with: u = 50 sat 


and x = 1 0 sat 


L 0-4 J 


l04j 


non-linear system simulation is run with no system failure. Very nice decoupled 
tracking is demonstrated. The baseline LQR simulation performance looks the same 






Figure 2-21: Decoupled SMC Tracking, Baseline Case 
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This plot shows the results of the LQR controller with failures The disturbance bias at 
5 sec has magnitude of 50 A non-linear system simulation is am with system failure at 
3 sec. Tracking performance is severely degraded. 






Figure 2-22: LQR Tracking with System Failure 
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and x = 75sat 


K0.4) 


1,0.4 J 


This plot shows the SMC controller with: u = 175 sat 

The disturbance bias at 5 sec has magnitude of 50. A non-linear system simulation is 
run with system failure at 3 sec. Very nice decoupled tracking is demonstrated. Notice 
the gains on the discontinuous control elements are higher than the baseline case. This is 
because a larger value is needed to provide the reaching condition in the face of the 
increased uncertainty. 






Figure 2-23: SMC Tracking with System Failure 
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This plot shows the SMC controller with: u = 50 sat 



( a ^ 

and x = 1 0 sat 


u = 50 sat 

— 


10.4J 


1.0 4 J 


The 


disturbance bias at 5 sec has magnitude of 50. A non-linear system simulation is run 
with system failure at 3 sec. Note the controller is still tracking x after the failure at 3 
sec (although some degradation is apparent), but the disturbance bias at 5 sec causes x to 
diverge. The significance of this plot is that it demonstrates that if the disturbance level 
is higher than designed for, the SMC will not be able to stabilize the system — and, in 
fact, may drive a statically stable state unstable 






Figure 2-24: SMC Tracking with System Failure (Insufficient p x & p 0 ) 
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This plot shows the SMC controller with: u = 175 

The disturbance bias at 5 sec has magnitude of 50. A non-linear system simulation is 
run with system failure at 3 sec. Note the controller is tracking x throughout, but 9 
breaks tracking shortly after the system failure The significance of this plot is that it 
demonstrates the ability of the controller to track one state even though another 
(coupled) state is not on its sliding mode. 


r o A 1 

f g b ^ 

— — and x = 1 0 sat 


lo.4j 1 

{OAJ 






Figure 2-25: SMC Tracking with System Failure (Insufficient pe) 
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Parasitic Dynamics and SMC 

3.1 Introduction 

As introduced in Chapter 2, unmodeled parasitic dynamics and physical limits are 
a significant implementation problem for sliding mode controllers. In order for the 
system to remain on the ideal sliding mode, an infinite frequency switching must occur. 
Real actuators and sensors prevent this infinite frequency switching, and the sliding 
mode is lost (along with all its desirable characteristics). In order to demonstrate the 
effects of actuator dynamics on a SMC system, the problem of the inverted pendulum on 
a translating cart is revisited. The single control input case with the regular form design 
is used. In all cases, the linear system model is used in the simulation There is no noise 
and no bias disturbance. In each case, the initial state of the cart position is x(0) = -1.0, 

and the initial state of the pendulum angle is 0(0) = -10 deg. The control law is: 

\ 

u = 50 sat — with the same sliding surface as the first example in Section 2.2 Three 

UU 

cases are presented. The only difference between them is the bandwidth of a simulated 

actuator: Figure 3-1: SMC Control, No actuator; Figure 3-2: SMC Control, With 

5 20 5 13 

Actuator: — (s) = , Figure 3-3 . SMC Control, With Actuator: — (s) = . 

u c s + 20 u c s + 13 
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Commanded Force (N) xande(m, rad) 









Commanded Force (N) x and 0 (m, rad) 






One approach to solve the actuator problem is to include the actuator in the initial 
design process. In this approach, the actuator is lumped together with the plant, and the 
sliding manifold is defined with the increased order associated with the addition of the 
actuator dynamics For example, consider the following system. 


u 

p 

s + p 

6 

- - ^ 

k 

y 

w* 

W 

s + a 

W 


Figure 3-4: Simple Plant with Actuator 


Including the actuator dynamics, the input-output differential equation for the system is 

y + (p + a) y + (p a) y = p k u (3 1) 

Define the sliding manifold as 

o= my + y (3 2) 


The order of this manifold is consistent with a traditional SMC — one degree less than 
the relative order of the state variable to be controlled. If the actuator had been 
neglected, the order of the manifold would have been zero instead of one. Next, find the 
equivalent continuous control, Ueq, needed to maintain ideal sliding motion by setting 
6 = 0 and solving for u (see Section 2. 1 .4) 


u 




fm(p + a)-lV fa'] 

/ ■- Y+ r y 

mpk J \kj 


(3.3) 


Substitute u,*, from (3.3 ) into the original differential equation (3.1 ) to obtain the 
closed loop, unforced system dynamics. 
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y = -~ y (3.4) 

m 

Note that none of the parameters (p, k, or a) appear, implying that, when on the sliding 
surface, the system output is invariant to changes in both plant and actuator parameters. 
Indeed, simulations of this simple system confirm that this is the case. An interesting 
note is that, while the system output, y, is invariant to changes in both the actuator and 
plant parameters, the intermediate state, 5, is only invariant to changes in the actuator 
parameters. This can become an issue if there are physical limits for 5 and can result in 
system instability if those limits are reached 

The simple solution to the actuator problem would appear to be the inclusion of 
the actuator in the SMC design. The resulting system is robust to both actuator and plant 
variations. Unfortunately, this is difficult to implement in practice because of the 
increased order of the sliding manifold associated with the additional dynamics In 
general, the order of the manifold will increase by the same order as the modeled 
actuator dynamics. This means, for a second order actuator, at least two derivatives of 
the output signal are required. For a real system with measurement noise, these 
additional derivatives make this approach very unattractive. While it is true that the 
SMC is very effective at rejecting the input noise, it does so by high frequency control 
commands. A noisy signal passed through two derivatives and an SMC would have 
highly undesired characteristics at the actuator input. 
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Young, Utkin, and Ozgiiner’s classic survey paper, “A Control Engineer’s Guide 
to Sliding Mode Control ,” 97 briefly discusses the actuator bandwidth issue. 

In plants where control actuators have limited bandwidth, e.g., 
hydraulic actuators, there are two possibilities: First, the actuator is outside 
the required closed-loop bandwidth. Thus the actuator dynamics become 
unmodeled dynamics. . . . While it is possible to ignore the actuator dynamics 
in linear control design, doing so in VSC requires extreme care. By ignoring 
actuator dynamics in a classical SMC design, chattering is likely to occur 
since the switching frequency is limited by the actuator dynamics even in the 
absence of other parasitic dynamics. Strictly speaking, sliding mode cannot 
occur, since the control input to the plant is continuous. 

Second, the desired closed-loop bandwidth is beyond the actuator 
bandwidth In this case, regardless of whether SMC or other control designs 
are to be used, the actuator dynamics are lumped together with the plant, and 
the control design model encompasses the actuator-plant in series. With the 
actuator dynamics no longer negligible, often the matching conditions for 
disturbance rejection and insensitivity to parameter variations in sliding 
mode which are satisfied in the nominal plant model are violated. 


They then offer several approaches to deal with the problem of parasitic dynamics 

• Boundary Layer Control. By increasing the boundary layer thickness, the 
effective linear gain is reduced, and the undesirable oscillations about the 
sliding manifold can be eliminated in some cases. While this method is 
often proposed, Young et al are quick to point out that this is not a 
recommended approach. A worst case boundary layer control design is 
usually required, resulting in poor SMC disturbance rejection properties. 
For a tracking task, this means poor tracking. Further, boundary layer 
control can not always stabilize the chatter. 

• Observer-Based SMC. An asymptotic observer is placed in the feedback 
path for the SMC. This observer acts as a high frequency bypass loop and 
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can eliminate chatter due to the unmodeled actuator dynamics Section 3 3 
investigates this approach further. 

• Disturbance Compensation. An SMC disturbance estimator is used. The 
control law consists of a conventional linear feedback component and a 
discontinuous component based on the estimated disturbance. 

• SMC Design with Prefilter. Actuator dynamics are incorporated as a 
prefilter to the SMC. A form of this approach is utilized by Shtessel and is 
demonstrated in Section 3.2. It is actually more of a postfilter, but the 
concept is similar. He does not refer to his method as this, but the 
approach falls into this general category 

• Frequency Shaping. The sliding manifolds are defined as linear operators 
and are “introduced to suppress frequency components of the sliding mode 
response in a designated frequency band.” 

Of these different approaches, the two which appear most straightforward and applicable 
to the problem at hand are the observer-based and postfilter-based designs. A large body 
of work is available in the literature utilizing the postfilter approach, so this method is 
investigated first. It has some drawbacks which are discussed in the next section The 
observer-based design is developed in Section 3.3. This approach is the one used in the 
final design methodology of this work. 

3.2 Postfilter Design with Dynamic Boundary Layer 

Shtessel and his associates have a number of papers addressing the issue of finite 

bandwidth actuators with rate and position limits 63 ' 66,106 ’ 110 Actuator dynamics are 

handled by defining the sliding manifold in terms of the derivative of the output states, 
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creating a pseudo command which is the derivative of the desired actuator command. 
Note, this is similar to the approach of a 2-SMC. The actual actuator command is then 
obtained by integrating the pseudo command through a model of the actuator. It should 
be emphasized that this approach requires measurement of the actuator output. Again, 
although not referred as such, this method is essentially an SMC postfilter. As will be 
demonstrated in the examples to follow, this method works very well for a first order 
actuator with a known bandwidth. In order to handle rate and position limits, they use 
what they call “reconfigurable sliding modes.” Probably a more accurate description is 
“dynamic boundary layer.” The main idea is to dynamically adjust the boundary layer 
thickness to keep the controller operating in the linear region of a boundary layer 
saturation element. By increasing the boundary layer thickness, the tracking 
performance is degraded but the states remain within the boundary layer, stability it 
maintained, actuator limits are not violated, and integrator wind-up is avoided. For 
known limits, this method works very well, as will be shown in the following examples. 

3.2.1 Application Example: Pit 

In order to exercise the 
ideas of Shtessel’s dynamic 
boundary layer and validate the 
operation of an example SMC, the 
exact same model from Shtessel’s 
work 63 is used. 


ich Rate Tracking for an F-16 



Chapter 3: Parasitic Dynamics and SMC 


3.2 Postfilter Design with Dynamic Boundary Layer 


Consider the linearized longitudinal short period mode approximation of an F-16 at 
Mach = 0.7, and h = 10,000 ft. 
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= A 


+ B6 f 


5 e = -20 5 e + 20 u 


(3.5) 


where 


A = A n + AA, B = B n + AB 


A n - 


«n = 


) 0.9937 


0.0400 

AA = 


- 1 .2600 


1.8560 

- 0.1770 " 

AB = 

' 0.0885 

- 19.5000 


9.7500 


l(t~5) 


l(t - 5 ) 


The “failure” simulated at t = 5 sec represents a 50% loss of horizontal tail area 


The pilot command, u p , consists of 0.1 rad/s pitch rate pulses of 1.0 sec duration with 
polarities of +, and + at times of 0.0, 3.0, 6.0, and 9.0 seconds respectively. These 
commands are filtered by a reference model which is given in state space as: 


" -1.2693 

0.9531 ' 

B r = 

' -0. 1 770 ~ 

-9.4176 

-5.7307 


-19.5000 


The designed SMC control law is: 

° = (q r -q) + 5(q r -q) + 25 J(q r -q) dx 


6 „ =6 + 


v20, 


f 


-10 sat 
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The block diagrams for the overall model and the SMC controller are given below. 



Figure 3-5: F-16 SISO Longitudinal System Block Diagram 



CE> 




Figure 3-6: F-16 SISO SMC Controller Block Diagram 


The simulation results for the healthy aircraft are shown in Figure 3-7. The results for 
the failed aircraft (failure at t = 5 sec) are shown in Figure 3-8. In both cases, actuator 
dynamics are included, but actuator limits are not. The results are excellent tracking in 
both cases. 
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Figure 3-7: F-16 SISO Pitch Rate Tracking, SMC, No Failure, No Limits 
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Next, actuator positions and rate limits 5 max =0.37 rad, 5 max =1 Orad/s are 

modeled. All the following plots use these limits. At this point, the boundary layer is 
static — Shtessel’s dynamic boundary layer is not in use yet. 

First, the case with same input profile as the baseline cases is given. The only 
difference between this case and the previous case is the inclusion of actuator limits in 
the simulation actuator models. The results are shown in Figure 3-9. Notice that some 
minor rate limiting occurs, but the position limit is never reached. The result is slightly 
degraded tracking performance. This tends to be typical — some rate limiting is usually 
tolerated. 

Next, in order to cause the system to hard limit, the input profile is changed to 
0.2 rad/s pitch rate pulses (instead of 0 1 rad/s). The results are shown in Figure 3-10. 
Once the actuator hits the position limits, the system no longer tracks the input and it 
becomes unstable. This also tends to be typical — with a static boundary layer and 
nonredundant control effectors, position limits are almost always fatal to SMC when 
encountered. 
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Shtessel el a / 63 - 66110 propose a method to account for actuator limits in the SMC 
design. They use a dynamic boundary layer. The boundary layer thickness, e, is 
calculated to satisfy the following inequalities which guarantees compliance with 
actuator displacement and rate limits and avoidance of integrator windup : 63 


Integrator Windup: a, <e, Vi = 1, 3 


(3.6) 


Actuator 
Deflection Limit: 



r fa Y1 

8,+0/a.) 

r- +p yj 


Vi = 1, 3 


(3.7) 


Actuator Rate Limit: 




/V 

G l 

V +p 

v eq, K , 

— 




Vi = I, 3 


(3.8) 


The next case is the same as the previous case except it includes this dynamic 
boundary layer The results are shown in Figure 3-11 The dynamic boundary layer 
keeps the system stable and achieves remarkably good tracking performance in the face 
of the actuator limiting. Tracking performance is degraded somewhat, but this is to be 
expected when limits are reached. 

Next, the actuator bandwidth is decreased Both the actual actuator model 
bandwidth and the design actuator model term in Eqn (3.7) are halved (as = 10). 
Results are shown in Figure 3-12. Results are comparable to the last case. 
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Figure 3-11: F-16 SMC, u p = 0.2 rad/s, With Failure, With Limits, Dynamic BL 
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Unfortunately, the actuator model term in Eqn ( 3.7 ) is a design parameter and 
has a static value. If the actual actuator dynamics do not match this, there could be 
problems. Figure 3-13 shows the results when the actuator model bandwidth is 
decreased but the design parameter in Eqn ( 3.7 ) remains the same. The system fails. 

This presents a problem when considering an application for a reconfigurable 
design After an actuator failure, the actuator dynamics are not known. Therefore, the 
very situation demonstrated here is likely to happen Also, the rate and position limits 
are not known after a failure. Although not shown here, it turns out that off-design 
limits also result in an unstable SMC. 

Another major drawback to this method is the assumption of a first order actuator 
in the control law. If a second order (or higher) actuator is used, the current control law 
fails. It seems plausible to re-derive the control law postfilter assuming a second order 
actuator. However, this would require the second derivative of the output states for the 
sliding manifold — which is unreasonable from a practical standpoint. Plus, this solution 
would only be valid for a second order actuator 

The conclusion is this, for actuators with known order, bandwidth, and limits, this 
method works exceptionally well. The results are the most impressive of any in the 
current literature. In fact, these results initially defined the direction of this research. 
However, the issues relating to changes in the true actuator dynamics have proven to be 
difficult to overcome. 
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3.2.2 Application Example: Bank Angle Tracking for a DC-8 



The problem of sensitivity to changes in actuator dynamics can be partially 
addressed using a multiple-loop scheme also proposed by Shtessel. 66 In this approach, 
an SMC is designed to track outer loop flight attitude angle commands. The output is a 
pseudo command for an inner loop SMC which tracks angular rates. In addition, a “very 
inner” loop is wrapped around the actuator and is controlled by another SMC. 

For simplicity in an initial demonstration, a linear lateral-directional model of a 
DC-8 is taken from McRuer, Ashkenas and Graham. 133 The goal is to design a bank 
angle tracker. In this case, the longitudinal variables are not controlled, and the only 
actuators are the ailerons and rudder. This results in a square system, as required. 
Sideslip angle is also not controlled, but the yaw rate is regulated to zero. 

Using the stability and control derivatives for flight condition 8003, the state 
space representation of the linear lateral-directional perturbation equations for the 
healthy DC-8 are: 133 
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The healthy actuator dynamics are assumed be 

5 5 r 20 

(s) = — -(s) = 

5 5 s + 20 


Outer Loop Design. Although an SMC can be used for the outer loop, for 
simplicity a basic loop shaping technique is used here. It is assumed the inner loop SMC 
has provided decoupled control of the roll rate, so the effective system seen by the outer 

p , % _ 4 

loop is simply the roll rate model reference transfer function: — 00 - g 2 + 8) 2 + 2 2 

The result of the loop shaping design gives the outer loop compensator: 


G c (s) = 1778 


(s + 2) 2 
(s + 20) 3 


Inner Loop Design: The sliding mode control is designed using the derivative of 
the actuator position, and then the actual actuator command is calculated using the 

model of the actuator dynamics: 8 = -A 6 (5-5 c ) . This means the system being 

controlled has relative order r = (2, 2} T and the sliding surfaces are given by 

a = // + Ctf + C3 J17 dx, cr,//e R 2 
where /? = {(p r -p), (r r -r)} T . 


(39) 


The control law is. 


5„=8 i+ — 
3 8 


a i 

P, sat — 
e 


Vi = 1,2 


(3.10) 


Notice, again, this is a postfilter approach assuming a first order actuator. The chosen 
design parameters are: 
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C 2 = diag{5}, C 3 = diag{25}, e=l, ps a = -50, ps r = 25 

Actuator Loop Design: The sliding surfaces are given by: 

o«.=(8--«.) + 10j(S=-S.)dx 

a„. =(8„ -6,)+loJ(5„ -8, ) dt 

and the SMC control law is 



fa, ) 


f° s 1 

u Sa = 1 0 sat 

M 

u 5r = 10 sat 

l 1 J 


A schematic of the overall system is shown in Figure 3-14. 



Figure 3-14: Schematic of DC-8 Lat/Dir Control System 

A schematic of the inner loop SMC is shown in Figure 3-15. 



delta a 



Figure 3-15: Inner Loop Sliding Mode Controller 
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A schematic of the actuator loop SMC is shown in Figure 3-16 (note both actuator 
controllers are the same; only one is shown). 
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Figure 3-16: Actuator Loop Sliding Mode Controllers 


Results for a case with no failure and no noise is given in Figure 3-17. The 
controller achieves excellent tracking. 

Next, a failure is simulated to test the system’s ability to handle a failure The 
“failure” is defined as 
System A-matrix 
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Results for a case with a simulated failure at / = 4 0 sec and no noise is given in 
Figure 3-18. Again, the controller achieves excellent tracking. In fact, there is no 
difference in the tracking performance between this case and the un-failed case. The 
multi-loop system is invariant to the changes in the A-matrix, the B-matrix, and even to 
the reduced bandwidth of the actuator. 

Unfortunately, the required control action is completely unrealistic. If actuator 
rate and position limits are considered, this system fails. Of course, if Shtessel’s 
dynamic boundary layer is employed in the inner loop, stability can be maintained in the 
face of actuator limits (with the expected degradation in performance) However, in 
order for the dynamic boundary layer to work, the limits must be known While this 
approach handles changes to actuator dynamics nicely, it still does not solve the problem 
of unknown limits. It also does not solve the problem.^' unknown actuator order. 
Again, a first order actuator is assumed. If a higher order actuator is present, the control 
law fails. Another drawback to this approach is the assumption that the states of the 
actuator (position and rate) are available. This could present actual implementation 
issues — especially in a reconfigurable setting. 
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3.3 Observer-Based SMC 

Recognizing the essential triggering mechanism for chattering is due 
to the interactions of the switching action with the parasitic dynamics, an 
approach which utilizes asymptotic observers to construct a high-frequency 
by pass loop has been proposed. This design exploits a localization of the 
high-frequency phenomenon in the feedback loop by introducing a 
discontinuous feedback control loop which is closed through an asymptotic 
observer of the plant. Since the model imperfections of the observer are 
supposedly smaller than those in the plant, and the control in discontinuous 
only with respect to the observer variables, chattering is localized inside a 
high-frequency loop which bypasses the plant However, this approach 
assumes that an asymptotic observer can indeed be designed such that the 
observation error converges to zero asymptotically 97 

This is a good qualitative description of how an observer-based SMC works. 
However, what is not clear is exactly why the observer helps and how to choose the 
appropriate observer gains. Initial simulation runs confirmed the utility of the observer- 
based design. The observer does, indeed, enable an SMC to run in the presence of 
unmodeled finite bandwidth actuators. However, the performance is highly dependent 
upon the speed of the observer (the observer gains). In general, if the observer gains are 
too high, chatter and instability result. Sensor noise also becomes a problem. If the 
observer gains are too low, robustness to system parameter variations is lost. 

In order to investigate this further, consider the following system. 


Model Control 

Reference Controller Distribution Actuators Plant 



Figure 3-19: Observer-Based SMC Block Diagram 


Chapter 3: Parasitic Dynamics and SMC 


3.3 Observer-Based SMC 







(3 


Define the following state space representations for the given system: 

xe9t n , x c e9t llc , x r efR nf , i 0 e9t n ° 

ye9? n \ y a e<R m , y c e9i my , _y r e9t"\ yetR™’', ze9t m * 

Actuators 

x a — A a x a + B a K b y c 
ya = x a + D a K b y c 

Observer 

x = (A 0 - G C z )x+ B 0 K b ^ c + G z 

y = c 0 x 

For the sake of analysis, consider an equivalent plant in which the control distribution, 
actuators, nominal plant, and observer are lumped together as shown below. 


Linear Plant 
x = A x + B y a 

j/ = c y x+ D y j> a 

z = C z x + D z y a 

Reference Model 

X r = A r *r + B r 'p 

y r = C r x r + D r r p 



Figure 3-20: Equivalent Plant with Observer 

It is easy to show that the state space representation of this system is given by 



A 

BC a 

^nxn 0 


BD a K b 

x e = 

®n a *n 

Aa 

^n a *n 0 

x e + 

BaK b 


GC Z GD z C a 

(a 0 -G c z )j 


B 0 Kjj + G D z D a 


yc 


y = 


c v o , . 

y m y x ( n a +n o> 


0 0 c ft 

mxn mxn O 

y y a 


D y K b 


y c 


X 

X 


a 

x 


( 3 . 12 ) 


It is instructive to now examine the transfer function X-(s) and compare it to both the 


original nominal plant transfer function with no actuator and the plant with the actuator 
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For the MIMO case, the transfer functions of interest are the ones from the control 
outputs to their corresponding output variable. Since a square system is assumed, this 
direct correspondence will always exist. The cross-coupling transfer functions are not 
needed here 

3.3.1 Observer-Based Design, SISO Case 

Consider the F-16 SISO system given in Section 3.2.1 with a second order 

8 20 2 

actuator j~( s ) - r . Recall the controlled feedback variable is pitch 

rate, q. The Bode plots of the nominal system and the nominal system with the actuator 
are shown in the Figure 3-21 While this plot does not show anything surprising, there 
are two things to point out when considering a sliding mode controller for this system. 
First, the relative order of the nominal plant is 1 . Therefore, if the actuator dynamics are 
neglected, the order of a traditional 1-SMC sliding manifold has to be 0 (1 degree less 
than the relative degree of the controlled state). Unfortunately, the real system has the 
actuator, so the actual relative degree of the system is 3. Therefore, the SMC can not be 
guaranteed to stabilize the system with the actuator. Second, if the actuator is included 
in the design, the sliding manifold would need to be 2 nd order. While an SMC could be 
designed to stabilize this system, the second derivative of the pitch rate would be 
required. Now, compare the estimated output from the equivalent system shown above 
with the nominal plant output. The observer poles are arbitrary selected to be very fast, 
(A. = -500, -501). The resulting Bode plots are shown in Figure 3-22. 
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Figure 3-21: F-16 Bode Plots, — Nominal System and Actuator 




Figure 3-22: F-16 Bode Plots, — Observer System Output 

y c 
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Since the controller is operating on the estimated state, the behavior of this transfer 
function is of particular interest. At frequencies below the bandwidth of the actuator, the 
observer state will follow that of the actual plant+actuator response. At higher 
frequencies, the magnitude of the observer state will approach that of the nominal state, 
but with a large increase in phase lag This system, including the observer, cannot be 
stabilized by the SMC. However, an interesting observation can be made. The 
estimated state signal once again has the same relative order as the original system As 
long as the nominal plant used in the observer has the same relative order as the nominal 
system plant, this will always be the case. Since it is critical to know the relative order 
of the plant for which the SMC is to be designed, using an observer provides a large 
benefit. This is also encouraging because the knowledge of the actual order of the 
actuator dynamics is no longer needed. In addition, no measurement of actuator output 
is necessary. However, just using an observer does not guarantee system stability. The 
observer must be tuned such that the effects of the actuator are attenuated before being 
fed back to the SMC. Now, consider slowing down the observer. 

The Bode plots in Figure 3-23 show the results for observer frequencies of 500, 
100, 50, 20, 10, 5, and 1 rad/s. On the magnitude plots, the effect of decreasing the 
observer speed is to cause the plots to approach that of the nominal system with no 
actuators. In general, the differences between the two cannot be completely eliminated, 
but they can be minimized to a degree depending on the relative dominant frequencies of 
the actuators and the plant. On the phase plots, note there is an observer speed below 
which the phase does not exhibit the large lag. In general, the observer poles must be 
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Frequency (rad/s) 


* y 

Figure 3-23: F-16 Bode Plots, - — Various Observer Speeds 

Yc 

below this frequency in order for the SMC to stabilize the system The phase lag/lead 
shape near the bandwidth of the actuator is common and results in a characteristic 
lead/lag shape in the closed loop. These Bode plots are very useful for determining the 
appropriate observer gains during the design process. First, design the SMC for the 
nominal plant neglecting the actuator dynamics. Next, include the actuator dynamics 
and examine the Bode plots of the equivalent system with the observer. Start with 
observer speeds near the bandwidth of the actuator and decrease the observer speeds 
until further reduction does not produce any further benefit. Double-check the results by 
examining a step response of the controlled variable of interest. 

As an example, consider again the F-16 SISO pitch rate tracking problem with the 
actuator defined above. The SMC is already designed. Now tune the observer 
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Observer speeds of 20, 10, 5, 1 and 0.1 rad/s are tried. The Bode plots and resulting 
closed-loop step responses are shown in Figure 3-24 and Figure 3-25. There is almost 
no difference on the Bode plots between observer poles at 1 and 0 1 Based on the Bode 
plots alone, observer poles at 1.0 appear to be the best choice. This design tends to give 
too much overshoot, however, and the final choice of 5 is made based on the actual step 
responses. 




10 -' 10 ° 10 ' 10 ’ 
Frequency (rad/s) 


Figure 3-24: F-16 Bode Plots, — Choosing Observer Speed 

Yc 
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Figure 3-25: F-16 Pitch Rate Step Responses, Various Observer Speeds 


3.3.2 Observer-Based Design, MIMO Case 

Assuming the system is square and feedback linearizable, the controlled outputs 

can be decoupled and independently managed by the controller. In light of the desire to 

shape the observer feedback loop appropriately for the SMC, it logically follows that an 

independent observer may be needed for each feedback channel. Indeed, simulation has 

indicated that this is the case, especially when the feedback variables have sufficiently 

different time-scales. The MIMO case then becomes simply a set of independent SISO 

cases, and the approach outlined above is used for each feedback channel Parameter 

coupling can still be somewhat of a problem, so some care must be given to ensure the 

observer for one channel does not cause adverse interaction with another channel. An 

example of this is seen later in the ICE aircraft application. 
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Consider, now, a system with multiple observers (in this case, three), 



Figure 3-26: Multiple Observer-Based SMC Block Diagram 


with the following state space definitions 


xe9t n , x a e9i na , x^eW"', x r e<R n ' , x o/ e9? n ‘> 1 , 

>/ c eW n \ y r e W 1 "* , y\ e 9t m > r , y 2 e 9t niy , 
z, eSR m *', z 2 z‘R m *\ z 3 e9? m * J , 


y 3 € 9t m ’ f 


Linear Plant 

x = A x + B y a 
_y = Cy X + Dy y a 
zy = C zl X + D zl y a 
z 2 = ^ z2 1 ®z2 Va 

= C z3 * + D z3 


Actuators 

x a = A a x a 1 ®a 

3V/ — c a jc^ + D a K|, y c 


( 3 . 13 ) 


Reference Model 

x r = A r x, + B r r p 

y, = C r x r + D r r p 

Coi e 9? my><n with zeros in the rows 
corresponding to states not output by 
Observer 1 . Gi = Observer gains 


Observer 1 

x 1 = (a o i - G j C z i )xy + B oI K b ^ c + G i z i 
h = c oi *i 
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Observer 2 

*2 - ( A o 2 “ G 2 C z2 ) x 2 + B o 2 K bXc + G 2 z 2 
h = C o2 *2 

Observer 3 

x 3 ~ ( A o3 _ G 3 ) x 3 + B o 3 KbXc + G 3 z 3 

h = C o3 *3 

Output Feedback 

y = yi+y2 + h 


C02 e 9? my * n with zeros in the rows 
corresponding to states not output by 
Observer 2 G2 = Observer gains 


C«3 6 9t my * n with zeros in the rows 
corresponding to states not output by 
Observer 3 G3 = Observer gains 


Again, an equivalent plant in which the control distribution, actuators, nominal plant, 
and observers lumped together as shown in Figure 3-20 can be expressed as 
x e - A e x e + B e y c 


G V x(n -tn.+n^+n,) 

J y ' a ol o2 o3 ' 


y= 
y=c e x e 

where 


x e + DyK^ y c 




x 

x a 

X 1 

x 2 

x 3 


(3.14) 
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Again, this system gives the designer the ability to look at the — (s) transfer 

y c 

functions and tune the observers as needed 

3.3.3 Observers with Failed Plants 

In general, it is seen that slowing down the observer helps mitigate the effects of 
parasitic dynamics. Consider now the behavior of the observer when the actual plant is 
significantly different than the nominal plant. In this case, it is assumed that the 
observer continues to use the nominal model in its calculations while the actual plant has 

y 

suffered some failure. The same ( s ) transfer functions used for tuning the observer 

Yc 

provide useful insight here. The Bode plots in Figure 3-27 show the SISO F-16 
example The nominal system and a failed system (a 50% loss of horizontal tail area) 
are shown. The actuator dynamics are neglected. Also shown is the equivalent system 
assuming a very fast observer (X = -500, -501). Note that at low frequencies, the 
observed state closely matches that of the actual failed plant. However, the observed 
state then transitions to the nominal system and matches the nominal system at 
frequencies above the observer frequency. This implies that as long as the observed 
state is close to the actual state around the closed-loop bandwidth, the controller will 
exhibit acceptable performance. 
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Frequency (rad/s) 


Figure 3-27: F-16 Bode Plots, Nominal Plant — , Failed Plant Fast Observer — 

y c y c y c 


Next, consider slowing the observer down to the speed determined earlier as 
optimum for reducing the effects of the parasitic dynamics (X = -5, -6). Again, the 
observer follows the nominal system at frequencies above the observer frequency. 
Below that frequency, the observed state follows neither. The result in the closed loop is 
very poor tracking and a steady state error. 
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Figure 3-28: F-16 Bode Plots, Nominal Plant — - , Failed Plant ^ fai1 , Slow Observer — — 

Yc . Ye Yc 


These results lead to the conclusion that a very fast observer is desired for 
robustness to system parameter changes. Unfortunately, fast observers do not help 
eliminate the effects of the actuator. Therefore, these competing design objectives must 
either be somehow balanced, or another method of dealing with parasitic dynamics must 
be sought 

3.4 Model Reference Hedging 

Observer-based SMC helps deal with the problem of unmodeled parasitic 
dynamics. However, the addition of the observer does not directly address the issue of 
control saturation The observer does help with rate saturation to some degree because 
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of phase lag reductions at medium to high frequencies. However, position saturation is 
still an issue — especially when the system does not have redundant control effectors 
Initially, in an attempt to deal with control saturation, a method called model reference 
hedging is pursued The concept of hedging has been successfully demonstrated in a 
dynamic inversion design approach . 134 The concept in words is this. “The reference 
model is moved backwards (hedged) by an estimate of the amount the plant did not 
move due to system characteristics the control designer does not want the adaptive 
control element to ‘know’ about .” 134 The actual accelerations are subtracted from the 
expected accelerations (assuming no actuators). This difference represents the amount 
of desired acceleration which was not achieved due to the actuators and should capture 
nonlinear saturation of the actual actuators. This difference is then subtracted from the 
reference model acceleration. Since all this is done in a dynamic inversion setting, these 
accelerations are pseudo-commands for the dynamic inversion controller, and Johnson el 
al. call this “Pseudo-Control Hedging.” The concept has great merit and is very 
successful in their work. 

In order to employ this method in an SMC design, some modification to the actual 
implementation is required. Rather than subtracting the acceleration difference from the 
reference model pseudo command, the following approach is used. The actual 
controlled variable output is subtracted from the expected output (through a nominal 
system with no actuators). This signal (yh) represents the amount of unachieved 
performance due to the actuator This passes through a hedge gain (Kh) and is subtracted 
directly from the model reference states as shown below. 
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If a 3 rd order reference model is used, the hedge signal enters the reference model as 
shown below. 



In order to investigate how hedging affects the system, an analytical expression 
for the system with hedging included is needed. Consider, initially, a model architecture 
as shown below. 
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Hedge Gain Hedge Plant 



Figure 3-31: Hedged System Architecture 


In this system, the hedge plant, Gh, is simply the nominal plant model (assuming it is 
stable). The subsystem block labeled “Hedged Reference Model” is the subsystem 
shown in either Figure 3-29 or Figure 3-30 Note that hedging occurs on individual 
control variable channels, so it is possible to write the transfer function for the hedged 
reference model. It can be shown that this transfer function (for a 2 nd order reference 
model) is given by 


x r h = 


m 0 

v 

( 2 

2s +mjs 

2 

^s +mjs + mo J 

X C 

2 

(s +mis + mo J 


K h Yh 


(3.15) 


sG r (s)x c -G f (s)K h y h 
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This indicates that the hedged reference signal consists of two parts, a reference model 
part and a hedge signal passing through a “hedge filter,” Gf. With this definition, the 
block diagram can be redrawn as: 



Figure 3-32: Equivalent Hedged System Architecture 


Next, define the following state space representations for the given system: 

jre 9t n , * a e9? n \ x c e 9t"', x r e 9t n - , x h eW nh , xye9t nf , i o/ e 9 t n <”, i o2 G9t n ^, i w e9t n ° J 

yes™', y a eW m , y c e9l m *, y [ e91 m > r , y h etR m >', y f e9? m >\ 22 6 3i m ', y 3 e <K m >' 

z,e9t m z\ z 2 e9t m ^ j z 3 e9t m '\ ;p h e9t m >' 


Actual Linear Plant 

x = A x + B y a 
y =C y x + D y y a 

z i — C z j x + D z j y a 
z 2 = x + D z 2 >a 
z 3 ~~ C z3 x + 


Reference Model 

*r = A r *r +B r 

y r =C r x r + D r x c 


Actuators 
x a — A a x a + B a 

y a — ^a a ^a^b .Vc 


Hedge Plant 
*h = A h^+ B h K b3'c 
yh =C h ^ + D h K b y c 


Linear Compensator 
x c — A c x c + B c c 
y c ~ ^ c * c + D c e 


Hedge Filter 
x f = A r x f + B r K h yh 
yf=C ( xf +D f K h >- h 


Chapter 3' Parasitic Dynamics and SMC 


3.4 Model Reference Hedging 









128 


Observer 1 

X 1 =( A oI - C,C tl )*, +B ol K b _y 0 +G| Z; 
h = c «i x i 

Observer 2 

x 2 = ( A o2 _ G 2 C z2 ) x 2 + B o2 +G 1 z 2 

h = c 0l i 2 

Observer 3 

x 3 ~ ( A o3 — G 3 G z 3 ) J 'i + B o 3 K b)'c + G 3 Z 3 

h = c o3 

Output Feedback 

h = yi+h+h+y/ 


C«i e 9 t myxn with zeros in the rows 
corresponding to states not output by 
Observer 1 . G| = Observer gains 

C 02 g 9t myxn with zeros in the rows 
corresponding to states not output by 
Observer 2. G 2 = Observer gains 

C 0 3 e 9t myXn with zeros in the rows 
corresponding to states not output by 
Observer 3. G 3 = Observer gains 


Again, this can be cast into an equivalent system in which the control distribution, 
actuators, nominal plant, observer, and equivalent hedge are lumped together as shown 
below. 


y 



Figure 3-33: Equivalent Plant with Observer and Hedge 
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The equivalent system is defined in state space as 


X 


x e = A e X e + B e y c 


y=\ 


+ »!,) 


Xp + 


[ D ! K b 


y c 


y^j - G e x e + yc 


x 


a 


X I 


x 2 

x 3 

X f 

x h 


(3 16 ) 


where 


A 

BC a 

°n*n 0l 

°nxn o2 

O nxn o3 

®nxn f 

°nxn h 

^n a xn 

A a 

° n a xn ol 

°n a xn o2 

°n a xn o1 

°n a xn f 

°n a xn h 

G 1 C zl 

Gj D z j C a 

(Aoi-g, c zi ) 

^ n ol Xn o2 

° n oi xn o3 

0 n o , x H f 

° n oi x »h 

G 2 C z2 

g 2 c a 

° n o2 Xn ot 

(A o 2-G 2 C z2 ) 

° n o2 xn o3 

°n o2 xn f 

° n o2 xn h 

G 3 C z3 

G 3 C a 

° n o3 Xn ol 

0 n 0 3 Xn o2 

(Ao3 “G 3 C^) 

°a o1 xn f 

° n o3 xn i, 

xn 

°n f xn a 

°n f xn ol 

°n f xn o2 

°n r xn o3 

A f 

Bf K h C h 

°n h xn 

° n h xn a 

°«h xn oi 

0 n h xn O 2 

° n h xn o3 

° n h x «f 

A i, 


B 


e ~ 


BD, K b 

B, «< b 

B 0 l K b +G( D zl D, K b 
B 0 2 K b +G 2 Djj D, K b 
B„ 3 K b +G 3 Dj 3 D, K b 
B f K b D„ K b 
Bh K b 


C e = 


°m y xn °m y xn a c 0 l C o2 C q3 C f D f C h 


D c = D f K h D h K b 
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Now the effects of hedging can be examined on the Bode plots of the transfer 

y h 

functions — ( s ) as is done with the observer. For the sake of illustration, consider a 

y c 

SISO model of the F-18/HARV longitudinal system for a pitch rate tracking task. 
Details for the system are given later in Section 4.5. The important aspects of the system 
for this discussion are that the nominal plant is stable and has a relative order of 1 for 
pitch rate; it has second order actuator dynamics; and it has a single observer which is 
very fast (X = -1000, -1001). The observer ensures the SMC sees a relative order of 1, 
however, it is running too fast to eliminate the instability caused by the parasitic 
dynamics of the actuators. Note, the F-16 model used in previous illustrations cannot be 
used here because it has an unstable plant — which causes yields an unstable hedge 
signal. 

The Bode plots below show the nominal system and the — ( s ) feedback signal. 

yc 

Also shown is the transfer function K|, Gh(s) Gf(s) which is the hedge signal. In this 
example the hedge gain is very low to show the system behavior with no hedging. Note 
the large phase lag. The SMC designed for this system cannot stabilize it In this 
example, the hedge plant, Gh, is simply the nominal system open loop transfer function 
for pitch rate, q. The hedge filter, Gf, is the transfer function defined earlier in 
Eqn. (3 15). 
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Figure 3-34: System with Hedging, Bode Plots — , K|,=0.01 

Yc 


Next, increase the hedge gain to see its effect on the system. This is shown in the 
Bode plots below. Of course, increasing hedge gain does not change the loop shape of 
the hedge signal — it simply shifts its magnitude curve vertically As Figure 3-32 
indicates, the hedge signal is additive to the observer signal. Therefore, as it approaches 
the observer magnitude curve from below, it increases the magnitude of this latter curve. 
The result is that the large drop-off in the magnitude curve (due to the actuators) is 
reduced. Note, also, that the phase no longer has the characteristic increase in phase lag 
This system is stable and has acceptable tracking performance. 
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Figure 3-35: System with Hedging, Bode Plots — , Ki,=0.8 

y c 

The hedge plant and the observer plant have the same relative order. This means 
that when the hedge gain is large enough to eliminate the drop-off in the magnitude 
curve, the additive effect of the hedge signal is too large at high frequencies. Note how 

the ~( s ) magnitude curve is above the nominal system at high frequencies. This can 

y c 

result in overshoot in the closed loop response. If the hedge signal had additional roll- 
off at high frequencies, this effect could be eliminated. 

After examining several different systems, it is noted that the hedge signal has 
basically the same loop shape for each system. It resembles a derivative at low 
frequencies; it peaks; and it rolls off at the relative order of the plant at high frequencies. 
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This is not surprising considering the hedge filter has the same form in all cases — it 
looks like a high-pass filter. All real plants have a high frequency roll-off When these 
two are placed in series, the result is this characteristic “hump” shape. The question then 
is this: can the same beneficial effects of hedging be achieved with a simplified 

“equivalent hedge” transfer function? This is highly desired because, in its current form, 
hedging can only be used with a stable plant 

Consider a hedge filter of the form (a high-pass filter). 

Gf (s) = — - — (3.17) 

s + af v ' 

And a hedge plant of the form (a low-pass filter): 

G h( s ) = ^ — — (3.18) 

s H-a^s + bh 

The hedge plant, Gh, in Eqn (318) has relative order 2 and is intended to be used with a 
system with a relative order of 1 The extra pole is added in order to have the additional 
high frequency roll-off, which is noted as being desirable in the previous discussion. 
This appears to be a good rule. The hedge plant should be one relative degree higher 
than the nominal system plant. Figure 3-36 shows a comparison of this simplified hedge 
model (Kh-Gh Gf) with a hedge signal using the nominal plant. In general, the simplified 
model can be designed to be fairly close to the original hedge system. In addition, 
desirable high frequency roll-off is obtained and lightly damped zeros are eliminated 
Note that since the simplified models operate on individual pseudo-command channels, 
the input to the hedge plant is simply y c (instead of B^). In other words, the take-off 
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point for the input to the hedge model in the block diagram in Figure 3-32 moves to a 
position before the control distribution block. The appropriate changes to the state space 
matrices (Eqn. (3.13)) are easily made. 
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Figure 3-36: Simplified Hedge System, — - 

Consider now, using this simplified hedge system with the F-18/HARV model 
demonstrated above. The results are shown in Figure 3-37. Compare this plot with the 
original hedge scheme results in Figure 3-35. The results here are very good In fact 
these results are better than the original hedge scheme. The observed state being fed 
back is very close to the original nominal plant with no actuators. 
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Figure 3-37: System with Simplified Hedging, Bode Plots K|,=8.0 

Ye 


3.4.1 Creating the Hedge Model 

The following technique is proposed for creating the hedge model. Begin by 

y h 

plotting the Bode plots of ( s ) and the nominal system as before with zero hedge gain. 

y c 

The basic form of the desired hedge model loop shape is as follows: 

• +20 dB/dec slope at low frequencies 

• -20*r dB/dec slope at frequencies where the actuators distort the magnitude 
curve (r = relative degree of nominal system with no actuators) 

• -20*r-20 dB/dec slope at high frequencies 
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This shape can be accomplished (for a system with relative degree 1) using the hedge 
plant and hedge filter given in Eqns (3.17 )and (3.18 ) with the pole in the hedge fdter 
at the high frequency end of the magnitude distortion and the two poles of the hedge 
plant at the low frequency end of the distortion. For example, in the F-18/HARV 

s 

example given above, the hedge filter and hedge plant are G f (s) = and 

s + 100 

G h (s) — — — ^22 . This is illustrated below in Figure 3-38. 

s + 20s + 1 00 



Frequency (rad/sec) 


Figure 3-38: Creating the Hedge Model, — 

y c 
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Obviously, the selection procedure of the pole locations for the hedge system is not 
exact. Some tuning is usually required. However, the basic rules of thumb given above 
are a good starting place. Some key factors which affect the optimum pole placement of 
the hedge system include: system plant parameters, actuator dynamics, and observer 
speed. Once the hedge system loop shape is set, adjust the hedge gain as required. 

3.4.2 Selecting Hedge Gain 

Note that the hedge gain in the simplified hedge model shown above is much 
higher than that in the original hedge approach (Kh = 8.0 instead of 0.8). The magnitude 
of the hedge gain required is a function of the bandwidth of the hedge plant and filter 
Higher bandwidth hedge systems require lower hedge gains; lower bandwidth hedge 
systems require higher hedge gains. In general, it is better to have a lower bandwidth 
hedge system because the larger range of acceptable hedge gains makes the overall 
system less sensitive to hedge gain— -thus easier to design and potentially easier to 
implement an adaptive hedge scheme. 

If the observer is very fast, or if no observer is used, there is a minimum hedge 
gain which will stabilize the system. Using this method of examining the Bode plots of 

y h 

^(S), it is easy to determine the minimum hedge gain. Examining the Bode plots will 

Yc 

also give an indication of the upper limit for acceptable hedge gain In general, if the 

y h 

hedge signal moves the magnitude plot of ( s ) above the nominal plant magnitude 

y c 

plot, the closed-loop system will have unacceptable overshoot and phase lag This can 
lead to instability if outer control loops are closed around the SMC system The 
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y h 

use of Bode plots of — ( s ) , combined with time domain step responses of the SMC 
inner loop control variables provides a useful tool for choosing hedge gains. 



Figure 3-40: System with Simplified Hedging, Step Responses 


3.4.3 Hedging in M1MO Systems 

Since the control variables are decoupled by the controller, it is a simple matter 
to treat each channel independently and design a hedge model for each separate channel 
However, recall that with the initial concept of hedging, the control output fed into a 
nominal model of the plant. This means that the original hedge signal included not only 
the primary variable transfer function, but it also included the effects of all the cross- 
coupled transfer functions as well. For example, in a lateral directional model (with roll 
rate, p; and sideslip angle, p), the primary transfer function for roll rate would be roll 

rate to roll-rate-command, — . The cross-coupled transfer function would be roll rate to 
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P-command, — When using the simplified hedge plant as introduced in Section 3.4 1, 
U P 

the effects of the cross-coupling are lost Unfortunately, if the effects of cross-coupling 
are strong, the interaction of commands of one channel with the unmodeled parasitic 
dynamics can lead to instability in another channel Therefore, in MIMO applications, 
the cross-coupled transfer functions need to be examined. If necessary, a hedge model 
for the cross-coupling term can be designed and the feedback loop properly shaped. The 
same guidelines introduced above are also used for the cross-coupled hedge models. 
This cross-term hedge signal is then added to the primary hedge signal For example, 
see the hedge model shown in Figure 3-41 This shows the hedge model for a lateral- 
directional system in which stability axis roll rate, P s , and sideslip, P, are the control 
variables. The cross-coupled terms are included. 



Figure 3-41: Hedge Model Example with Cross-Coupled Terms 
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Examples of hedging with cross-coupling are given later in an F-18/HARV and a 
tailless fighter aircraft application The potential effects of neglecting these cross- 
coupled terms are also shown. 
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3.4.4 Observer Gains and Hedging with Failed Plants 

Recall from Section 3.3.3 that robustness to failures of the system plant and 

tolerance of parasitic dynamics present competing objectives when choosing observer 
gains. Fast observers are more robust to system parameter variations, but they do not 
eliminate the adverse interaction of the SMC and parasitic dynamics Slow observers 
enable the SMC to operate in the face of parasitic dynamics, but they provide poor 
robustness to system parameter changes. Some trade-offs must be made when designing 
the observer. 

The use of hedging can sometimes help meet both objectives. Since the hedge 
system contains no nominal system model, it is less sensitive to system parameter 
variations. The basic loop shape of the hedge system is always the same. Therefore, as 
long as the parameter variations do not significantly change the location of the distortion 


y h 

in the — ( s ) Bode magnitude curve, the design hedge system poles will provide the 
y c 


correct hedge loop shape for the failed system With the correct hedge system loop 
shape, the effects of the parasitic dynamics can be mitigated with the correct hedge gain. 
After a system failure, the hedge gain will probably need to be adjusted using some kind 
of adaptive mechanism. 

This is encouraging and means that it is possible to run the observer with fast 
eigenvalues. In fact, one may be tempted to eliminate the observer altogether and use 
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only hedging. In some cases, this is true. However, in general, hedging alone is not the 
best answer To illustrate the effects of observer speed and hedging, consider several 
different types of failures, plant parameter variations (changes to system state space A 
and B matrices), actuator bandwidth variations; and pure time delay The following set 
of figures use the F-18/HARV SISO longitudinal pitch rate tracking model introduced 
above and given in detail in Section 4.5. All use the same SMC and include actuator 
dynamics. 



0 0.5 1 1.5 2 2.5 3 3.5 4 4 5 5 

Time (sec) 


Figure 3-42: F-18/HARV Step Response, Nominal System, Nominal Actuators, 

No Time Delay 

Figure 3-42 shows the nominal system step response assuming a “fast” observer with 
hedging and a “slow” observer with no hedging. In all the following cases, “Fast 
Observer” means observer poles at X = -100,-101; and “Slow Observer” means observer 
poles at X = -4,-5. While the SMC alone cannot stabilize the system when actuators are 
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included, the observer-based SMC and the hedged SMC model both stabilize the system 
as illustrated in Figure 3-42. Note that they provide equally good performance. The 
question, then, is. which will provide the most robustness to system failures? 
Figure 3-43 shows the same SMC, observer, and hedge systems with a “failed” system 
plant. 



Figure 3-43: F-18/HARV Step Response, FAILED System, Nominal Actuators, 

No Time Delay 

The failure is a 50% loss of horizontal tail area (which affects the system state space A 
and B matrices). As can be seen in Figure 3-43, the hedged system is much more robust 
than the system with the slow observer. While this tends to be true in the case of system 
variations, it is not true if the actuator dynamics change. For example, Figure 3-44 
shows the results for a nominal plant and failed actuators. The actuator failures are a 
50% reduction in the actuator bandwidths. 
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Figure 3-44: F-18/HARV Step Response, Nominal System, FAILED Actuators, 

No Time Delay 

As can be seen, in this case the slow observer has better performance than the hedged 

y h 

system. This is because the location of the distortion in the — ( s ) Bode magnitude 

y c 

curve occurs at a lower frequency and the hedge signal loop shape is incorrect for the 
failed system. Varying the hedge gain does little to help. The hedge system poles 
require tuning for the new condition. 

Next, consider the case of a nominal plant and actuators and a pure time delay A 
50 ms time delay is inserted between the controller and the actuators For the 
simulation, a second order Pade approximation is used. 
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Figure 3-45: F-18/HARV Step Response, Nominal System, Nominal Actuators, 

50 mS Time Delay 


In this case, neither the slow observer nor the hedged system with a fast observer 
provide responses comparable to that of the nominal system. The slow observer has 
20% overshoot, and the hedged model is too oscillatory. As illustrated, increasing the 
hedge gain eliminates some of the oscillations, but changing the hedge gain alone can 
not restore the desired performance. The addition of a time delay has the effect of 

y h 

moving the location of the distortion in the ( s ) Bode magnitude curve to a lower 

y C 

frequency. Again, this type of failure requires tuning of the hedge model poles. 

Finally, consider a case that combines all three types of failures. Figure 3-46 
shows the case of a failed plant, 50% reduction in actuator bandwidths, and a 50 ms time 
delay. Results are clearly unacceptable. 
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Figure 3-46: F-18/HARV Step Response, FAILED System, FAILED Actuators, 

50 mS Time Delay 

As illustrated, the slow observer system is unstable, the hedged system is highly 
oscillatory. In order for the slow observer system to be effective, the plant model used 
in the observer needs to be adjusted. In order for the hedged system to be effective, the 
hedge model poles need to be tuned Hedge gain alone will not help 

The conclusion here is that design trade-offs need to be made based on the types 
of failures expected In general, hedging provides more robustness to plant variations 
Slow observers provide more robustness to actuator bandwidth variations and time 
delays. If a system is to be robust to all three types of failures, the observer should be 
moderately fast and hedging should be included. In addition, the hedge model can be 
tuned such that there is a wider separation between the low and high frequency poles in 
the hedge model (brand ah in Eqns (3.18) and ( 3.17 ) ) than is required for the nominal 
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system This results in a hedge model which provides degraded performance in the 
nominal case but is more robust to actuator variations and time delays. 
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3.4.5 A Hedge By Any Other Name 

It is worth noting that the hedging method proposed no longer resembles the 
concept of hedging given in the initial reference work. 134 Note in Figure 3-32 that all the 
hedge signal dynamics occur parallel to the observer loop and no longer enter the 
reference model The hedging method here is really a form of observer loop shaping. In 
fact, it is recognized that this approach is very similar to Loop Transfer Recovery (LTR) 
as used in Linear Quadratic Gaussian (LQG) control An LQG controller combines an 
optimal regulator (LQR) with an optimal observer (Kalman filter) and can result in 
arbitrarily poor stability margins. LTR is a method by which optimality of the 
observer is traded for increased stability margins. The hedge signal used here also has 
the effect of “tuning down” the observer at certain frequencies in order to recover 
desirable stability margins. 

It can also be shown that the hedge system is attempting to invert the actuator 
dynamics. This is why it is so dependent on the actuator bandwidths. When viewed in 
this light, one could argue that this is simply another form of an SMC prefilter-type 
design. In order to prevent confusion, the method presented here should probably be 
given a different name. On the other hand, recall that if the hedge signal enters the 
reference model as originally illustrated in Figure 3-29 (which is supposed to implement 
the intent of hedging 134 ), the final method being proposed is exactly equivalent. 
Therefore, this analysis has illustrated in the frequency domain why model reference 
hedging works and how to proceed with the design. There are two advantages of this 
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approach First, the actuator states are not required as they are in pseudo-command 
hedging . 134 Second, if the simplified hedge plant is used, the model of the nominal 
system is not required. 

3.5 Robustness Issues 

So far, the effects of asymptotic observers and model reference hedging on an 
SMC system have been investigated, and design procedures for synthesizing these 
models have been offered. General observations about performance and robustness are 
made via application examples and inferences from frequency domain Bode plots. A 
key issue which remains unanswered is a formal quantitative statement about overall 
system robustness There are several difficulties which hinder the derivation of such a 
proof. These are introduced below. 

3.5.1 Difficulties with SMC Robustness Proofs 

3.5.1. 1 Practical design issues with standard Lyapunov proofs 

Almost all robustness proofs for SMC involve the use of Lyapunov stability 

criteria to guarantee global attractiveness of the sliding manifold (i.e. the reaching 
condition) A classic simple example of how such a proof proceeds is given in The 
Control Handbook 92 Consider the following nonlinear plant: 


x = A(x) x + B u(x) 



The sliding manifold is given by 

o(x) = [ Si s 2 ] x 
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A variable structure control law is used 


where 


u(x) = k|(x) xi + k 2 (x) x 2 


fa,(x), if a(x) x, > 0 
jp, 00, ifa(x)x,<0 


A candidate Lyapunov function is selected 



If the control law constants are selected such that ( 3 19 ) is satisfied, the sliding mode is 
guaranteed 

Note that the criteria in (3.19) only accounts for maximum excursions of the 
state variables and the (known) nonlinearities in the system. It does not account for any 
other parameter variations. A common method to include parameter variations is to 
write the state space system as 
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'/ 


x — (A 4 - AA) x + (3 + AB) u 

where AA and AB are unknown, bounded, smooth perturbations There are several 
examples of works which approach the problem in this manner 63 66 The analysis 
proceeds as above; however, it becomes extremely involved. In the end, bounds for the 
control law parameters which guarantee the reaching condition are found in terms of 
maximum state excursions and norms of the perturbation parameterizations. 

These types of proofs are necessary to complete a full analytical treatment of 
SMC robustness, however, they are not entirely practical. The example given above 
from The Control Handbook is a simple pendulum for which all the possible states of the 
system are known (and bounded). In general, for an aircraft application, the maximum 
excursions of all the states (and potentially their derivatives) are unknown and could be 
practically unbounded. Further, parameterizing all potential failure events is a difficult 
task Norms of the perturbation matrices for the state parameters is common. However, 
parameterizing actuator failures is especially difficult due to the nonlinear nature of the 
failure modes. From a practical standpoint, these issues make even this traditional 
approach to proving robustness very difficult 

3.5.1. 2 Boundary layers 

The Lyapunov analyses discussed above are generally applied to the pure sliding 
mode problem. The SMC control law is derived to guarantee the sliding motion, and 
once on the sliding mode, the system is invariant to matched uncertainty. When a 
boundary layer is included, pseudo-sliding takes place, and the invariance property is 
lost. It is possible to perform the Lyapunov analysis using a boundary layer in the 
control law definition. In this case, stability to the boundary layer surface (rather than 
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the sliding manifold) is guaranteed However, a question which remains is “How robust 
is the system inside the boundary layer?” The system can be made arbitrarily close to 
invariant by decreasing the boundary layer width, but it is not clear how much 
robustness is retained when practically useful boundary layer widths are used Further, 
boundary layer thickness appears to have competing robustness objectives. On one 
hand, very small boundary layers make the system more robust to plant parameter 
variations On the other hand, large boundary layers make the system more robust to 
unmodeled parasitic dynamics So, even if a complete formal Lyapunov stability proof 
is performed, the measure of robustness is not exactly known when a boundary layer is 
employed. 

3.5.1.3 Matching conditions 

Another difficulty with SMC robustness proofs is the issue of matching 
conditions The system is shown to be invariant to matched uncertainties when on the 
sliding surface. However, an issue which is not usually addressed is robustness to 
variations in the unmatched parameters. There are two recent works which provide 
design approaches for systems with unmatched uncertainties , 136 ’ 137 however, a formal 
treatment of robustness to unmatched uncertainties in a traditional SMC setting has not 
been found. There are two interesting observations concerning matching conditions 
which have been noted in the process of this research. First, matching conditions change 
depending on the state space selected. Second, if parasitic dynamics are neglected in the 
design, matching conditions vary with the bandwidth of the neglected dynamics Each 
of these is illustrated in turn. 
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First, the state variables chosen affect which parameters are matched and which 
are unmatched Consider again the simple system introduced in Section 3.1, Figure 3-4, 
pg. 87. By lumping the actuator with the plant, the resulting state space variables appear 
in phase variable form: 


*•' 


0 1 

x i 

+ 

"o" 



-pa -(p + a) 

_ X 2 _ 


k 


* 




(3.20) 


By inspection, it is expected that the parameters p, k, and a are matched because they lie 
in the row of the control input. The sliding manifold is defined in terms of the state 
variables (xi=plant output; X 2 =derivative of plant output): 

cr = m x 2 + x, (3.21 ) 


The equivalent control is found 



' m (p + a) - P 

Y 4- 



l mpk } 

A ^ • 



(3.22) 


Using this equivalent control, the closed loop, unforced response is determined 


X 1 
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with eigenvalues at X = 0, . 


m 

changes in parameters p, k, and a. 


Indeed, as noted before, this system is invariant to 
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Next, consider the same system; but instead of writing the state variables in phase 

variable form, write the states with the actuator and plant in series 

xl [-a k If x 1 fo 
= + u 

L° -PjL X aJ Lp_ 

(3.24) 

y =[* of 

L x . 


By inspection, it is expected that the parameter p is matched, and the parameters k and 
a are unmatched. The sliding manifold is defined in terms of the state variables 
(x=plant output, Xa^actuator output): 

a = mx a +x (3 25 ) 



Using this equivalent control, the closed loop, unforced response is determined 



( 3 27 ) 


with eigenvalues at X = 0, - 



The closed loop response is a function of both a 


and k. This system is not invariant to these unmatched parameters. This shows that 
when the sliding manifold is defined in terms of different state variables, invariance can 
be lost and system robustness is affected. The work of Chan, et al., 136 takes advantage 
of this by defining a subspace where unmatched parameters become matched— resulting 
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in a new sliding manifold on which invariance is achieved. Another work defines the 
existence of sliding surfaces in terms of Linear Matrix Inequalities. 137 The conclusion of 
this observation is that care must be given to the selection of state variables In this 
work, by specifying the system architecture to be square, all the key system parameters 
are matched However, there can still be uncontrolled states which remain unmatched, 
and robustness to changes in these parameters is not clear. 

A second observation is that matching conditions vary with the bandwidth of 
parasitic dynamics when these dynamics are neglected in the design Consider a second 
order plant model (with no actuator) 



( 3.28 ) 

y = D o][ Xl 
L X 2 J 

The sliding manifold is defined in terms of the state variables (xi=plant output; 
X2=derivative of plant output): 

o = m x 2 + x, ( 3.29 ) 

The equivalent control is found 

\ 

x 2 (3.30) 

) 


u eq =-T a l x l + a 2 


Using this equivalent control, the closed loop, unforced response is determined 



(3.31) 
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with eigenvalues at X- 0, . As expected, the system is invariant to changes in a t , 

m 

a 2 , and k. Now, a first order actuator model, — , is placed in series with the plant, 

s + p 

resulting in the following state space system. 

0 1 ° ITx,] [o' 

-a, -a 2 k x 2 + 0 u 

0 0 Pj |_ x a J |_P. 

( 3 32 ) 

X >1 


The actuator dynamics were not included in the design, but consider driving the “actual” 
system in ( 3.32 ) with the control derived for the system in ( 3 28 ). The equivalent 
control in ( 3 30 ) is not capable of maintaining the sliding mode because of the 
unmodeled disturbance, so a traditional discontinuous control element is included in the 
feedback control law. 


y = [i o o] 





+ 




-psgn(o) 


( 3.33 ) 


Further, in order to analyze this system using linear methods, replace the signum 
function with a boundary layer in which the limits of the saturation function are never 
reached. The approximated linear control takes the form 



(3.34) 
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The closed loop, unforced response is then given by 
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( 3 . 35 ) 


The characteristic equation for this system is 


s 3 + (a 2 + P) s 2 + 


p ppmk 

— + — + a| 

m £ 


S + B£i = 0 


( 3 . 36 ) 


First of all, notice that by neglecting the actuator and using the control law derived for 
the nominal plant, the resulting system is not invariant to either the plant parameters or 
the actuator parameters The relative values of the plant parameters and actuator 
parameters can result in oscillatory behavior and instability. Unfortunately, it is invalid 
to consider the system behavior as the bandwidth of the actuator, p, approaches zero 
because the assumption of a linear control becomes invalid The system, with a 
discontinuous control element in the control law, quickly becomes unstable as p is 
reduced; while the continuous control law results in closed loop eigenvalues at the 
locations of the nominal plant. However, as p is increased, the linear assumption holds 

and it can be shown that the closed loop eigenvalues approach: X = - — ; - — ^ m ; - oo . 

ms 


In fact, if the system is exactly on the sliding mode, the eigenvalues as p->oo are 


X = ;0;-oo, which is exactly the same result as obtained when there was no 

m 

actuator. In conclusion, actuator bandwidth affects the matching conditions On the low 
frequency side, matching conditions are lost and instability occurs. On the high 
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limit. 
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3.5.2 Difficulties with Hybrid Nature of the Observer-based SMC 

As discussed, proving robustness is difficult for the sliding mode controller. The 

nonlinear control necessitates the use of Lyapunov analyses which are problematic in a 
practical sense; boundary layers introduce an unknown amount of robustness 
degradation; and there are questions surrounding robustness to unmatched uncertainties. 
Now add an asymptotic observer with hedging. The classical linear tools for analyzing 
the robustness of the observer cannot be used due to the nonlinear controller. The 
already difficult methods of analyzing the SMC become intractable due to the addition 
of an entire closed-loop feedback structure in the SMC feedback path. The controller 
never “sees” the actual plant outputs, so it is no longer clear which parameters are 
matched, if any. 

3.6 Robustness Summary 

It is clearly unsatisfying to offer a design approach for which quantitative 
measures of robustness are not known — especially when the main goal of the method is 
to handle large scale failures. All efforts to quantify the robustness have turned up 
empty due to the issues described above. However, experience has shown that the SMC 
approach provides superior robustness and deserves continued development. While the 
formal robustness proofs are left for another, attention is now given to the development 
of a practical design procedure and application examples which demonstrate the SMC’s 
robustness. 
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Design Procedure 

4.1 Interpretation of Linear SMC Design in the Frequency Domain 

Considering a linear plant or vehicle model, the SMC system including boundary 

layers can be interpreted as a linear system and analyzed in the frequency domain That 
is, classical loop-shaping principles can be applied to the design, both in terms of 
choosing appropriate sliding manifolds and in evaluating the characteristics of the final 
SMC configuration. Note that, if a boundary layer is employed, and the boundary layer 
thickness is selected such that the limits of the saturation element are never reached, the 
control law is essentially a high gain linear controller. For example, consider a sliding 
function of the form 

a = Cj e + Cq e + c_| j edt (4.1) 


and an SMC control law (with a boundary layer) 

u = psat 






(4.2) 


If the saturation element always operates in the linear region, the feedback control law in 
the Laplace domain looks like a PID controller: 


M P P 
u(s) = — o = - 


C| S H-Cq + 


s 


e(s) 


(4.3) 
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This fact can be used to help choose both the sliding function constants (c.i, Co, and Ci) 
and the control law gain (p/s). This interpretation of the sliding function allows the use 
of familiar loop shaping techniques and provides insight into the behavior of the 
controller 

The full closed-loop system including the reference model, linear compensator, 
actuators, plant, observer(s), and hedge is given in state space by 


I CL = A CL I a. + B CL , y 


y- 


C 0 

y m * (n + n -»« +- n f aj +n+n+n) 

y a ol o2 o3 f h c r 


X CL ~ 


CL 


(4.4) 


where 


‘CL 


A, B.D.SC, B t C ( -B,D SD C, B t D SC r 
B t SC e A ( -B c SD C c B c SC r 

0 0 A. 


B f 


'B. D c SD r 
B c SD r 

B 


S s (l in + D c D c ) 1 ; and all e-subscripted variables are defined in Eqn (3.16). 


4.2 Instantaneous Adaptation 

One can rightly argue that, in treating the SMC as a linear high gain controller as 
introduced in Section 4.1, the system should no longer be called a sliding mode system. 
No true sliding mode is ever achieved, and all the properties associated with the sliding 
mode are lost If the system is nothing more than a linear one, why not use classical 
linear design techniques for which design methods and robustness bounds are well 
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known? There are several answers to this question. First of all, it is common practice to 
still treat this as an SMC. While this is a very unsatisfying response, it is true. When 
applying SMC to real systems with physical actuators, a continuous control signal is 
almost always sought. Since an infinite frequency switching is required to maintain the 
sliding mode, none of the researchers who employ a boundary layer tally have a sliding 
mode system. They really only have a high gain linear controller. Yet, these controllers 
are still referred to in the literature as sliding mode controllers. Another reason is that 
the sliding manifold provides known error dynamics for the system While a pseudo- 
sliding only provides asymptotic convergence to the sliding manifold, it nevertheless is 
converging to a known dynamic response. A third and more important reason is that this 
observer-based system with a pseudo-sliding mode is actually adapting to parameter 
changes. It instantaneously adjusts the gain and crossover frequency of the equivalent 
linear controller A conventional linear controller (designed by loop shaping, for 
example) does not do this. 

In order to illustrate this phenomenon, consider again the F-16 SISO longitudinal 
system from Section 3 2.1. Figure 4-1 shows the Bode plots of the pitch-rate to 
pitch-rate-command transfer function for the nominal system and a failed system. The 
system failure is a 50% loss of horizontal tail area (resulting in changes to the plant A 
and B state space matrices) Notice that the failure results in a gain reduction in the 
open-loop system transfer function. 
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An observer-based SMC with hedging is designed for the nominal model. The SMC 

M c 

control law is u c = -30 sat - , o = 1 (q r -q) + 20 j (q r -q) dx . The observer poles are 

1 J J 


( 64 V s > 

at X = -500,- 501 . The hedge system is —50 — . Using 

yS + 16 s + 64 JyS+300y 

Eqn ( 4.7 ), the closed loop transfer function, Gcl(s), is determined for both the nominal 
and failed cases. Next, determine the equivalent unity feedback loop transmission as 
illustrated below 
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The unity feedback loop transmission, L, is given simply as 


L( s ) = 


Gq.fr) 


(45) 


The resulting Bodes plots of the equivalent unity feedback loop transmissions for the 
nominal and failed cases are shown in Figure 4-2. 




Frequency (rad/s) 


Figure 4-2: Bode Plots, F-16 Equivalent Unity Feedback Loop Transmission 
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Note, first of all, the desirable — characteristic at the crossover frequency. Note also, 

s 

that both the nominal and failed systems look almost identical at crossover. Indeed, both 
systems have very good (almost identical) tracking performance. 

Consider, next, the determination of an equivalent serial compensator as illustrated 
below. 


+ 


►o — ► 


= 7 > 



where G a (s) is the actuator and P(s) is the system plant (nominal or failed). The element 
Gc e (s) then represents a serial compensator which is equivalent to the combined effects 
of the (linear) SMC, observer, and hedge. It can be calculated as 


G„(s) = 


L(s) 


G a (s) P(s) 


(46) 


This equivalent serial compensator is calculated for the nominal and failed systems The 
results are shown in the Bode plots in Figure 4-3. First of all, note that the two 
compensators are different. This means that the observer-based SMC has changed. 
Note also that the equivalent serial compensator for the failed case has an increased gain 
to compensate for the lower gain of the failed system. So, the controller is varying in a 
manner that tends to compensate for the changes in the vehicle dynamics. This behavior 
is attributable to the presence of the observer(s) in the feedback loop Of course, the 
changes are instantaneous. 


Chapter 4: Design Procedure 


4.2 Instantaneous Adaptation 







uotjvjdvp y stiOD uvjuvjsuj £ p 


djnpsoou j uSisdQ :p uajdtnjj 


ui uMoqs are °r)) pu e ( TO c i* 0 * ° 0 ) J° suoipury jqjsu^j) doo| uado aqx 

sjspiuBJBd uiajsAs SuiAjba aqi qjiM 3 §uBqo iou saop qoiqM joiBsuadiuoo padcqs dooj 
b SuiuSissp 01 snoSo|BUB si siqx !UB[d pspBj aqi qiiM JoiBsuaduioo |BU 3 S iU 3 [BAinb 3 
^psidBpB,, aqi Suisn qijA\ sqnsaj sqi ajBdiuoo puB iUB|d pajiBj sqi qi|A\ joiBSusdiuoo 
[BU 9 S 1 U 9 |BAinb 3 JBUIUIOU 3 qi guisn J 9 piSU 03 1 S 3 J 91 UI JO SI UO| 1 BAJ 3 SqO 1 SB| 3 U 0 

poyirpT, a 

3 r> pirn °£> ‘jojt;siiadui(r) |BU3S |U3|BAinbg[ 9|-j ‘sioij apog :£-£ 34n§;j 


(s/pej) Aouanbajj 




01 

(Q 

ZJ 

c 

CL 

CD 


CL 

CD 


*91 




Figure 4-4. The “adapted” equivalent serial compensator is superior. It has a 
considerably higher crossover frequency with little degradation in stability margins. 
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Figure 4-4: Bode Plots, Nominal Equiv Serial Compensator with Failed Plant 


4.3 Design Procedure Summary 

A proposed SMC design methodology can now be offered A step-by-step 
procedure is given, followed by three application examples of increasing complexity 
Recalling that the objectives of SMC design are finding appropriate sliding surfaces and 
switching logic, the procedure to be described approaches this problem from a 
frequency-domain based, model-reference approach. 
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(1) The vehicle model is obtained, along with an estimate of the frequency beyond 
which parasitic dynamics (or unstructured uncertainties) are likely to come into 
play. This frequency is referred to as the limit frequency in this discussion 
Actuator dynamics are not included in the nominal plant model. 

(2) A reference model is chosen for each control variable channel. Since piloted 
flight control is of interest in the present application, this reference model should 
be selected with an eye towards Level 1 handling qualities with no pilot-induced 
oscillation (PIO) tendencies. This can be accomplished based upon a pilot 
model-based handling qualities and PIO prediction technique. 138 An example is 
available. 43 

(3) The desired feedback structure of the control system is determined with a square 
system architecture. For example, if a pitch-rate command flight control system 
is desired, then pitch rate (q c ) becomes the output of the reference model, and 
estimated pitch rate (q) is fed back to the SMC system from the observer 

System error is then defined as e(t) = q c (t) — q(t) . Special care must be given 
to ensure the dynamics of the uncontrolled variables remain stable This 
assumption is analogous to the minimum phase assumption used in feedback 
linearization 66 Also, if there are redundant control effectors for the desired 
moments/pseudo-commands, a control distribution matrix must be defined This 
research considers only static control allocation. Simple ganging is used in the 
first two examples. A pseudo-inverse method is used in the third application 
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(4) The sliding manifold, a, is chosen based upon the following principles 

a) a is derived from a tracking error expression as 

a = e(t) P_l +K p _ 2 e(t) P “ 2 +K 0 e(0 + K_,fe(l) dt (4.7) 

where p is the relative order of the system, i.e., the number of 
times the vehicle output must be differentiated for the input to 
appear. Note that the (p-l) st derivative of the error signal is 
included in the definition of a. An integral term also appears in 
Eqn ( 4.7 ) to counter the steady-state bias often created with the 
use of a boundary layer. 

b) Recognizing that a boundary layer is to be implemented, the 
control law is expressed as a linear transfer function as discussed 
earlier. 

u(s) = — a 

e 

= k p( sP "' +k p 

The parameters K, are chosen to provide desirable properties in 
the frequency domain. This means creating a loop transmission 
with broad K/s-like characteristics around crossover. 139 This will 
always be possible since enough derivatives are included in 
Eqn.( 4.3 ) to create exact K/s characteristics beyond a certain 



frequency (at least as high as the limit frequency). Parasitic 
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dynamics are deliberately excluded in this formulation. This step 
will involve obtaining an estimate of K p , as this value will 
determine the crossover frequency of the loop transmission This 
crossover frequency is selected to provide acceptable stability 
margins as obtained from a Bode plot of the loop transmission but 
using a value of K p at least as large as the largest amplitude limit 
of any of the control effectors. The latter criterion is included to 
accommodate maximum trim positions of the control effectors. 
As opposed to typical designs involving loop shaping, very high 
crossover frequencies may result from this step. Indeed these 
frequencies may be well beyond the limit frequency. This result 
is of no immediate concern If a 1VQMO system is being designed, 
a classic sequential loop closure technique is used, thereby 
sequentially and independently determining the coefficients for 
each sliding manifold 

(5) Using the Kfs just determined in the definition of the sliding function, the 
existence of a sliding mode is verified in the inner loops using a true SMC. This 
step is completed without the observer, actuators, reference model or pilot 
model, i.e., assuming that no outer-loop is being utilized. If necessary, p is 
increased until sliding behavior is created. The initial value of p = K p obtained 
in step (4b) should be considered a lower limit in this process. While an 
analytical approach to determine p is certainly possible here, a more expedient 
route of establishing the sliding mode using a computer simulation of the system 
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is also possible. Near perfect tracking in the face of large parameter variations 
should be observed. The control signal, however, will exhibit very high 
frequency switching 

(6) A boundary layer is included in the controller by replacing the signum function 


sign(a) with the saturation function satl — J . While maintaining an approximate 


constant — = K p , increase the boundary layer thickness, £, until no high 

s 

frequency switching is evident. Again, a simulation of the SMC system is a 
convenient way of finding this s. Near-perfect tracking (with a continuous 
control signal) in the face of large parameter variations should be observed. 

(7) Parasitic dynamics are included in the model. The SMC controller will very 
likely be unstable at this juncture 

(8) An asymptotic observer is created for each control variable channel as discussed 
in Section 3.3. 

(9) If desired, a hedge model is designed as described in Section 3 .4 

(10) The frequency domain characteristics of the open and closed-loop SMC system 
with observer, boundary layer and reference model are examined to ensure that 
stability of the linear system is in evidence 


4.4 Task-Dependent Analytical Handling Qualities Assessment 

The goal of the research described herein involves the stability and performance 

robustness of flight control systems associated with damaged aircraft/actuators With the 
exception of uninhabited vehicles, all these aircraft will be under piloted control. Thus, 
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an analysis of the predicted handling qualities of damaged aircraft/actuators with the 
SMC systems is of considerable importance in determining the utility the proposed 
design methodology. Space does not allow a detailed description of the handling 
qualities assessment technique to be utilized here. Several key works provide a 
description of the pilot-model based assessment tool that can be used. 43 ' 138,140 Using the 
handling qualities level and PIO levels discussed therein, task-oriented handling qualities 
and PIO levels can be predicted using the Structural model of the human pilot 

In this pilot modeling approach, a Handling Qualities Sensitivity Function 
(HQSF)— defined as the magnitude of a transfer {unction obtained from the pilot model— 
is obtained. The magnitude is plotted on linear amplitude vs. frequency axes and 
compared with bounds delineating handling qualities Levels I -3. The predicted 
handling qualities level is determined by the minimum level bound violated by the 
HQSF. 

A similar plot involving the power spectral density of a proprioceptive signal in 
the pilot model is used to delineate predicted PIO levels PIO "levels" are defined 
relative to a six-point PIO rating (PIOR) scale 138 with PIOR = 1 meaning no PIO 
tendencies and PIOR = 6 meaning very PIO prone 

It should be emphasized that once a suitable pilot model has been chosen for the 
nominal or undamaged vehicle, it is not altered when pilot/vehicle performance is 
predicted for the damaged vehicle. The design examples to be discussed next include 
these handling qualities and PIO assessments. 
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4.5 Application Example: F-18/HARV Longitudinal SISO Model 



The F-18 High Angle-of- Attack Research Vehicle (HARV) is a modified pre- 
production-model F-18 Hornet used by NASA Dryden Flight Research Center for high 
angle of attack studies The aircraft completed 385 research flights during a three-phase 
program from 1978 to 1996. Key modifications included a multi-axis thrust vectoring 
system, advanced flight control laws, and movable strakes. 141 The only special feature 
utilized in the models which follow is the thrust vectoring system. 

The linearized longitudinal short period approximation for the nominal F-18/HARV 
at the flight condition, M = 0.6, h = 30k ft, straight and level, is given in state space as 142 


a 

= 

'-0.5088 

0.9940 " 

a 

+ 

'-0.9277 

-0.01787' 

5 stab 

q_ 


-1.1310 

-0.2804 

_q_ 


-6.575 

-1.525 

5 tv _ 


where: a = angle of attack (deg), q = pitch rate (deg/s), 5 sta b = stabilator deflection (deg), 
and 5 tv = pitch thrust vector deflection (deg). 

The failure for the longitudinal model of the F-18/HARV is defined as a 50% loss 
of stabilator area. This will be referred to as the plant failure. Later, actuator failures 
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will also be defined. An approximate state space representation of this failed system at 
the same flight condition can be given by 
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a 


'-0.48 

0.997 ' 

a 


-0.4639 

-0.01787 

_q_ 


-0.39 

-0.182 

_q_ 

+ 

-3 2875 

-1.525 



1 


U stab 


5 


L U TV J 


(4.10) 


The healthy actuator dynamics and limits are given in Table 4-1 



Dynamics 

Amplitude 

Limit 

Rate Limit 

Stabilator 

(30) 2 

(s 2 + 2 0.707-30s + 30 2 ) 

+30 deg 

60 deg/s 

Pitch Thrust Vector 

(20) 2 

(s 2 + 2 • 0.6 • 205 + 20 2 ) 

±30 deg 

60 deg/s 


Table 4-1: F-18/HARV Longitudinal Actuator Dynamics and Limits 


A pitch rate tracking control system is desired. The next step in the design 
process is to define a reference model. For this application, the following reference 
model is chosen. 40 


G r (s) = — (s) = 


100 


(s + 5) (s + 20) 


(4.11) 


Next, the desired square system feedback architecture is defined. This is given in 
Figure 4-5. The inner loop is a pitch rate feedback loop which is controlled by the SMC. 
A pilot closes an outer loop to control the Euler pitch angle. 
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Figure 4-5: F-18/HARV Longitudinal System Block Diagram 


The controller to be designed commands pitch rate This demand is allocated to the two 
control effectors using an arbitrary ganging method The control distribution matrix is 
defined as 40 


B 


k 



(4.12) 


The pitch-attitude command profile is a sum-of-sines profile which gives a random-like 
tracking task with frequencies representative of those encountered by a pilot 143 

ec(0 = ^Z A - sin ( (D - t ) (4.13) 

1 = 1 


The overall amplitude gain is £ = 5. The sum-of-sines parameters are given in 
Table 4-2. In addition, the pitch-attitude signal is introduced in exponential fashion by 

0 8 

multiplying the sum-of-sines signal by a unit step passed through a filter of . 

s + 0.8 

This results in a continuous, zero-mean, random-like signal with a zero initial condition 
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i 

A, 

©i 

i 

-1.00 

0.19947 

2 

1.00 

0.49867 

3 

1 00 

0 89760 

4 

0.50 

1.39626 

5 " 

-0 20 

2.39359 

6 " 

0.20 

4.19970 

7 

-0.08 

8.97598 

Table 4-2: 

Sum-of-Sines Parameters 


The measured states are angle of attack, a, and pitch rate, q The state measurement 
noise for each channel is assumed to be band-limited white noise fdtered by 


100 

s 2 + 20s + 100 


The RMS value of the noise is approximately 0.5 deg/s. 


Assuming the inner loop SMC provides the desired model reference tracking, the 
apparent system model that the pilot “sees” is the model reference. Therefore, the pilot 
model, 138 ’ 144 ’ 145 based on the reference model given, is shown below. 



Figure 4-6: Pilot Model for F-18/HARV Pitch Angle Tracking Task 
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compensator and plant with these values is shown in Figure 4-7. Note the desired — 

s 

shape at crossover. The crossover frequency is about (0 = 250 rad/s. The initial design 
of the SMC is complete and is shown in Figure 4-8. 



Figure 4-8: SMC Pitch Rate Controller for F-18/HARV 


The Fifth step in the design process is to verify that the SMC is working properly. 
Either replace the saturation element in Figure 4-8 with a signum function or make the 
boundary layer very small (s wO.OOOl). Figure 4-9 shows the time response for the inner 
loop pitch rate tracking of the SMC with a boundary layer of e = 0 0001 . The Simulink® 
simulation is am with an ODE2 solver using a fixed time step of At = 0.0005 sec. A 
system plant failure occurs at t = 10 sec. 
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Figure 4-9: F-18/HARV Inner Loop SMC Pitch Rate Tracking 


As expected, the SMC performs very well and is invariant to the system parameter 
changes. The control output shows the classic high frequency switching Note, if a 
signum function was actually used, this plot would exhibit the expected near infinite 
frequency switching. 

Next, the boundary layer is increased until a continuous control signal is achieved 
For this model, e = 0 1 is the approximate minimum boundary layer thickness that will 
result in a continuous control signal. The resulting time history of the inner loop 
tracking (with a failure at t = 10 sec) is shown in Figure 4-10. 
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Figure 4-10: F-18/HARV Inner Loop SMC Pitch Rate Tracking with Boundary 

Layer 


Although not apparent, the tracking performance is degraded slightly by increasing the 
boundary layer thickness. The performance is still excellent, even in the face of a large 
system failure. More importantly, the control signal is now continuous. During the 
design process, the boundary layer thickness, e, is varied without changing the SMC 
gain, p, which appears to be contrary to the guidance given in Section 4.3, Step 6 (of 

trying to maintain a constant K, = — ) This is because the amplitude limits on the 

1 e 
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actuators are ±30. The magnitude of p should always be equal to or greater than the 
maximum expected actuator amplitude limit. Therefore, in this case, |p| can not be 
reduced below 30. In addition, as long as e is less than 1, the transmission loop 
crossover will be at least as high as that designed for. Therefore, decreasing e below 1 
without varying p is usually acceptable. If e is increased above 1, p should be increased 

accordingly to try to maintain an approximately constant K =— . In practice, the 

£ 

easiest way to identify the minimum boundary layer thickness is to begin with p = K () 
and e = l. Then decrease s (leaving p alone) until the undesired switching begins. It is 
rare that e needs to be greater than 1 for the nominal design. 

The SMC design is now complete, the resulting tracking performance is 
remarkably good, and the system is (nearly) invariant to system parameter changes. The 
next step is the inclusion of the actuators. This is the point where sliding mode methods 
falter. Figure 4- 11 shows the inner loop tracking results with the actuators now 
included The system (without plant failure) goes unstable in less than a second. 
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Figure 4-1 1: F-18/HARV Inner Loop SMC Pitch Rate Tracking, Actuators 

Included 


The next steps are to design the observer and hedge models. As introduced in 
Section 3.3, the Bode plots of the equivalent plant/observer with unity feedback ( — ) 

u c 

are examined for varying observer speeds (with no hedging). Observer speeds of 100, 
50, 20, 10, 5, and 1 rad/s are tried The Bode plots are shown in Figure 4-12. 
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Figure 4-12: Bode Plots, F-18/HARV Pitch Rate — , Various Observer Speeds 

U c 


As shown on the Bode plots, observer speeds of 50 and 100 rad/s result in an unstable 
system. The unit step response for the inner loop with observer speeds of 20, 10, 5, and 
1 rad/s are given in Figure 4-13. Based on the Bode plots and the unit step response 
plots alone, a good choice for the observer speed is between 5 and 10 rad/s. 


Chapter 4 * Design Procedure 


4.5 Application Example: F- 1 H/H A RV Longitudinal SISO Model 





182 



Figure 4-13: Unit Step Response, F-18/HARV Pitch Rate, Various Observer 

Speeds 

Observer poles of X = -7,-8 are selected, and the design (with no hedging) is complete. 
Now exercise the full system with measurement noise and the pilot in the loop. The 
outer loop pitch angle tracking performance for the nominal vehicle is shown in 
Figure 4-14 Performance is very good. The large time delay between the commanded 
pitch angle and the achieved pitch angle is almost entirely due to the pilot. Note in 
Figure 4-15 that the inner loop tracking is also very good. Of course, with the actuators 
and observer now included, the SMC does not give the almost perfect tracking achieved 
before. Figure 4-16 and Figure 4-17 show the actuator deflections and rates, 
respectively No limits are hit As shown in Figure 4-18 and Figure 4-19, Level I 
handling qualities are expected with no adverse PIO tendencies. 
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Pitch Rate, q (deg/s) Pitch Angle, 0 (deg) 
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Figure 4-14: Outer Loop Tracking, Nominal F-18/HARV, with Actuators, 

Observer, and Noise 



Figure 4-15: Inner Loop Tracking, Nominal F-18/HARV, with Actuators, 

Observer, and Noise 
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Figure 4-16: Actuator Deflections, Nominal F-18/HARV, with Observer and Noise 



(/) 



Figure 4-17: Actuator Rates, Nominal F-18/HARV, with Observer and Noise 
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Figure 4-18: Nominal F-18/IJLARV Handling Qualities Assessment, 0 Tracking 

Task 



Figure 4-19: Nominal F-18/HARV PIO Assessment, 0 Tracking Task 






Next, exercise this control system with a vehicle failure. The failure used in the 
following simulations is defined as a plant failure (50% loss of horizontal tail area) plus 
a thrust vector hard-over of 5tv = -3 deg. Initially, the simulation starts with a healthy 
vehicle. Then, at t = 10 sec, the system experiences a sudden “unknown” failure. No 
changes are made to the control law or observer after the failure. 

Figure 4-20 shows the outer loop pitch angle tracking performance. Figure 4-21 
shows the inner loop pitch rate tracking performance After the failure, the inner loop 
tracking is noticeably degraded, although it is still relatively good. The degradation of 
the inner loop performance results in degraded outer loop performance. Figure 4-22 and 
Figure 4-23 show the actuator deflections and rates. Note the stabilator is working much 
harder after the failure. Still no limits are reached Even though the tracking 
performance is degraded, Figure 4-24 and Figure 4-25 predict Level II handling qualities 
(almost Level I) and no adverse PIO tendencies. This is actually remarkably good 
considering the severity of the failure. 
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Figure 4-22: Actuator Deflections, F-18/HARV, with SMC, Observer, and Noise, 

Failure at t = 10 sec 



1 1 ) 



Figure 4-23: Actuator Rates, F-18/HARV, with SMC, Actuators, Observer, and 

Noise, Failure at t = 10 sec 
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As discussed earlier, the slow observer can result in noticeably reduced robustness 
to plant parameter changes. Therefore, the observer gains are increased to help deal with 
this failure This, however, results in adverse interaction between the SMC and parasitic 
dynamics. With a fast observer, hedging is required. The methods described in. 
Section 3.4 are used to redesign the observer speed and build the hedge model. The final 
values chosen are 

Observer poles: A, = -50, -51; Hedge model: -5 — - — 

P V s + 1 00 

and the resulting design Bode plots are shown in Figure 4-26 




Figure 4-27 and Figure 4-28 show the tracking performance with this new design. Again 

the plant failure plus thrust vector hard-over occurs at t = 10 sec. Figure 4-29 and 

Figure 4-30 show the handling qualities and PIO bounds for the failed system. 
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Figure 4-28: Inner Loop Tracking, F-18/HARV, with SMC, Actuators, Observer, 
Hedging, and Noise, Failure at t = 10 sec 


n 
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This design, with the faster observer and hedge model included is clearly superior to the 
slow observer design for this particular failure With this design and failure, tracking is 
almost as good as the nominal case, and Level 1 handling qualities are predicted. 

As seen in Figure 4-26, the hedge model for this system is designed to work well 
with the known linear actuators. The portion of the hedge signal on the Bode magnitude 
plot which needs to have the -20 dB/dec slope is fairly narrow — just enough to remove 
the effects of the actuators. Unfortunately, if the actuator dynamics slow down (by 
changing bandwidth or by reaching rate limits), or if a pure time delay is introduced, the 
hedge model shown above will not be able to stabilize the system. A slower observer 
will help with robustness to changes in actuator dynamics but will result in lower 
robustness to system parameter changes. Therefore, a compromise is sought. Consider 
slowing down the observer to reduce the amount of hedging needed and also changing 
the hedge signal loop shape. Create a hedge signal loop shape with a broader region of - 
20 dB/dec slope. Since slower actuator dynamics are of concern (not faster), extend the 
region of -20 dB/dec slope toward lower frequencies. For example, consider the 
following observer and hedge systems. 


s 1 

f 4 ^ 

v s + 100 J 

v s 2 +4s + 4. 


Observer poles. X - -40, -4 1 , Hedge model: - 100 

The resulting design Bode plots are shown in Figure 4-3 1 As seen in 
Figure 4-31, the effects of the actuator dynamics are not removed as well as they are in 
the previous design. The hedge gain could be increased to help remove the distortion in 
the Bode magnitude plot; however, since the hedge loop shape has such a broad region 
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Figure 4-31: Bode Plots — , F-18/HARV Modified Observer & Hedge Design 

u c 

of -20 dB/dec slope, any more increase in hedge gain will cause the low frequency 
portion of the — (s) Bode magnitude plot to move above the desired nominal plot This 

u c 

results in unacceptable overshoot in the nominal case. Overall, this design appears to be 
a reasonable compromise. The nominal case with this new design has very slightly 
degraded performance and still has Level I handling qualities. 

Consider a new failure condition (Failure #2) as follows failed plant; stabilator 
rate limit = 30 deg/s; thrust vector rate limit = 6 deg/s, thrust vector 
bandwidth = 10 rad/s; and a 50 ms time delay in both actuators. The simulation results 
with a failure at t — 10 sec follow 
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Figure 4-32: Outer Loop Tracking, Failed F-18/HARV, with SMC, Actuators, 
Observer, Hedging, and Noise, (Failure #2) 



Figure 4-33: Inner Loop Tracking, Failed F-18/HARV, with SMC, Actuators, 
Observer, Hedging, and Noise, (Failure #2) 
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Figure 4-34: Actuator Deflections, Failed F-18/HARV, with SMC, Observer, 

Hedging, and Noise, (Failure #2) 



(0 



Figure 4-35: Actuator Rates, Failed F-18/HARV, with SMC, Observer, Hedging, 

and Noise, (Failure #2) 


Chapter 4: Design Procedure 


4.5 Application Example: F-18/HARV Longitudinal SISO Model 






197 



Figure 4-36: Failed F-18/BLARV Handling Qualities Assessment, 0 Tracking Task, 

With Hedging (Failure #2) 



Figure 4-37: Failed F-18/HARV PIO Assessment, 0 Tracking Task, With Hedging 

(Failure #2) 
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Several points are worth noting. The tracking performance is degraded after the 
failure, but the system remains stable In fact. Level II handling qualities (almost 
Level l) are predicted. This is very encouraging considering the severity of the failure 
and the fact that the thrust vector is in almost constant rate saturation and the stabilator 
hits its reduced rate limits. 

For comparison, consider a loop shaped design. Using classic loop shaping 


techniques, a controller is designed: G c =-2.5 


V +0.8s + 1.2 A 
s(s + 0.3) 


The resulting forward 


loop transmission, G c (s)P(s), with no actuators is given in Figure 4-38. 




Frequency (rad/s) 


Figure 4-38: Bode Plots, F-18/HARV Pitch Rate Loop Transmission for Loop 

Shaped Design 

As expected, this results in excellent tracking and Level I handling qualities for the 
nominal case. These plots are not given here, but they are essentially the same as the 
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SMC design results. This design is remarkably robust. In fact, for Failure #1, this loop 
shaped design does almost as well as the SMC design. However, this design can not 
handle Failure #2 as illustrated in Figure 4-39. Failure occurs at t = 10 sec. 



Figure 4-39: Outer Loop Tracking ( 9 ) for Loop-Shaped Controller, Failed 
F-18/HARV, with Actuators and Noise, (Failure #2) 


In conclusion, an SMC controller is designed for a SISO system with limited 
bandwidth actuators and is found to be superior to a classically designed system for a 
class of failures. Design tradeoffs are made in the observer and hedge models which 
help make the system more robust to different types of failures. Next, a basic MIMO 
system design is demonstrated 
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4.6 Application Example: F-18/HARV Lat-Dir MIMO Model 

The linearized lateral 
directional (Lat-Dir) state space 
model for the F-18/HARV 142 at the 
flight condition, M = 0.2, 
h = 1 Ok ft, straight and level, 
otirim = 29.7 deg, is given in Eqn. 

( 4.16 ). Note this low speed, high angle of attack, flight condition is a very demanding 
flight condition with extreme yaw/roll coupling. 




-0.05904 

-5.513 

0.06838 


0.4959 

-0.9391 

0.02632 


-0.8703 

0.6655 

-0.1038 


P 


P 

r 


+ 


0.005629 

1.879 

-0.1092 


0.005764 

1.328 

-0.09645 


0.003685 

0.02922 

-0.08404 


0 

0.6754 

0.006811 


0.0904 

0.217 

-2.974 


'DT 


'A1 


'RU 


RTV 


' Yrv , 


(4.16) 


where 

P — sideslip angle (deg), p = roll rate (deg/s), r = yaw rate (deg Is), Sqt 55 differential tail 
deflection (deg), 5 A i = aileron deflection (deg), 8 ru = rudder deflection (deg), 8rtv = roll 
thrust vector deflection (deg), and §ytv 2 yaw thrust vector deflection (deg). 
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The failure for the lateral-directional model of the F-18/HARV is defined as 


A fail = 2 0 A 
B fail = 0.75 B 


(4 17) 


where A is the state plant matrix and B is the control power matrix in Eqn ( 4. 16 ). 

This will be referred to as the plant failure. Later, actuator failures will also be defined. 
This failure is not intended to represent an actual physical failure. It simply gives a 
100% change in the system A-matrix and a 25% reduction in all actuator control power 
derivatives. 


The healthy actuator dynamics and limits are given below. H2 



Dynamics 

Amplitude 

Limit 

Rate Limit 

Differential Tail 

(30) 2 

i (s 2 + 2 0.707- 30s + 30 2 ) 

+ 17.5 deg 

60 deg/s 

Aileron 

(75) 2 

(s 2 +2 0 6-75s + 75 2 ) 

±27.5 deg 

100 deg/s 

Rudder 

(7 2) 2 

(s 2 +2-0.69-72s + 72 2 ) 

±30 deg 

100 deg/s 

Roll Thrust Vector 

(20) 2 

(s 2 + 2-0.6-20s + 20 2 ) 

±30 deg 

60 deg/s 

Yaw Thrust Vector 

(20) 2 

(s 2 + 2 -0.6 -20s + 20 2 ) 

±30 deg 

60 deg/s 


Table 4-3: F-18/HARV Lat/Dir Actuator Dynamics and Limits 
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A roll rate tracking control system with zero sideslip is desired. The next step in 
the design process is to define a reference model For this application, the following 
reference model is chosen: 43 


G p (s) = £s -(s) = 
Pc 


1000 

(s + 4.1)(s 2 + 2-0. 829-15. 6s + 15.6 2 ) 


(4.18) 


Next, the desired square system feedback architecture is defined. This is given in 
Figure 4-40 The inner loop has a roll rate and a sideslip angle feedback loop which are 
controlled by the SMC. A zero sideslip is commanded, and a pilot closes an outer loop 
to control the Euler roll angle. 



Figure 4-40: F-18/HARV Lateral-Directional System Block Diagram 


The controller to be designed outputs commands for sideslip angle and stability axis roll 
rate. These demands are allocated to the five control effectors using an arbitrary ganging 
method. The control distribution matrix is defined as 43 
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0 0.6 
0 1 


B 


k 


1 0 

0 0.6 

0.6 0 


(4.19) 


The input command profile for roll angle is a sum-of-sines function defined in 
Eqn. (4 13 ). The sum-of-sines parameters are those given in Table 4-2, and the overall 
amplitude gain in Eqn. (4.13 ) is £, = 6. 

The measured states are sideslip, P, roll rate, p; and yaw rate, r The state measurement 
noise for each channel is assumed to be band-limited white noise filtered by 

— and multiplied by a gain. Roll and yaw rate noise have a gain of 0 1414, 

s 2 + 20 s + 100 

resulting in an RMS noise value of approximately 0 07 deg/s. The sideslip noise has a 
gain of 0.0707, resulting in an RMS noise value of approximately 0 035 deg. 

Assuming the inner loop SMC provides the desired model reference tracking, the 
apparent system model that the pilot “sees” is the model reference. Therefore, the pilot 
model, 138,144,145 based on the reference model given, is shown below. 



Figure 4-41: Pilot Model for F-18/HARV Roll Angle Tracking Task 
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The next step is to define the desired sliding manifolds. The system (without 
actuators) has a relative order of 1 for both roll rate and sideslip. Therefore, the form of 
the sliding manifolds are 

a = K 0 e + K_ 1 J c dr; e = {p r -(3, p, - p } (4.20) 

The control laws, assuming the use of a boundary layer, then can be expressed in lineai 
form as 



This gives one zero to place and a gain to set for each loop during the loop-shaping 
design. Since the roll rate is a higher bandwidth loop, begin with this loop first 
Assuming the beta loop is left open, the values chosen for the p-loop are K pp = 500, 
Ko P = 1, and K. ip = 50. The resulting forward p-loop transmission of the compensator 

and plant (no actuators) with these values is shown in Figure 4-42. Note the desired — 

s 

shape at crossover. The crossover frequency is about 0) = 1000 rad/s. 

Using a standard sequential loop closure technique, the p-loop is closed and the 
P-loop is designed. The values chosen for the P-loop are K p p = 3000, Kop = 1, and 
K.ip = 5 The resulting forward P-loop transmission of the compensator and effective 

K 

plant with the p-loop closed is shown in Figure 4-43. Note the desired — shape at 

s 


crossover. The crossover frequency is about co = 200 rad/s. 
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Figure 4-42: p-Loop Transmission of Compensator* Plant, F-18/HARV, 

p-Loop Open 




Figure 4-43: p-Loop Transmission of Gompensator*EfTective Plant, F-18/HARV, 

p-Loop Closed 
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The initial design of the SMC is complete and is shown in Figure 4-44. 



Figure 4-44: SMC Roll Rate and Sideslip Angle Controller for F-18/HARV 

The fifth step in the design process is to verify that the SMC is working properly. 
Either replace the saturation elements in Figure 4-44 with a signum function or make the 
boundary layer very small. Figure 4-45 and Figure 4-46 show the time responses for the 
inner loop p s and P tracking of the SMC with a boundary layer of s = 0.000001. The 
Simulink® simulation is run with an ODE2 solver using a fixed time step of 
At = 0.0005 sec. A system plant failure occurs at t = 10 sec. Although the P input will 
be a constant zero for the final system, a simple sine wave is used here to verify 
decoupled tracking of p s and p. Notice that, in this case, the SMC gain, p, is reduced for 
the very low value of s. Sometimes, when the SMC gain is very high, the numerics of 
the solution become unstable. One solution is to reduce the simulation time step 
Another solution is to reduce p It may not be possible to reduce p in a manner such that 

K p =— = constant because, as discussed in Section 2.1.6, there is a minimum p which 
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As expected, the SMC performs very well and is invariant to the system parameter 
changes. The control output shows the classic high frequency switching. Note, also that 
the controller achieves decoupled tracking of p s and P, as desired. Considering the high 
level of yaw/roll coupling at this flight condition, the decoupled tracking performance of 
this controller is remarkable 

Next, the boundary layer is increased until a continuous control signal is achieved. 
For this model, e = 1 is chosen for both channels. The resulting time history of the inner 
loop tracking for p s and P is shown in Figure 4-47 and Figure 4-48. 




Figure 4-47: F-18/HARV Inner Loop SMC Roll Rate Tracking with Boundary Layer 
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Figure 4-48: F-18/ELARV Inner Loop SMC ft Tracking with Boundary Layer 


Again, the performance is excellent, even in the face of a large system failure. Also, the 
control signal is now continuous. Note the large spike in the control signal at the time of 
failure and the increased overall magnitude of the required control signal after the 
failure. 

The SMC design is now complete, the resulting tracking performance is 
remarkably good, and the system is (nearly) invariant to system parameter changes The 
next step is the inclusion of the actuators. Figure 4-49 and Figure 4-50 show the inner 
loop tracking results with the actuators now included. As expected, the nominal system 
goes unstable in less than half a second. 
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SMC Control Output, 
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Figure 4-49: Inner Loop p-Tracking, Nominal F-18/HARV, SMC with Actuators 




Figure 4-50: Inner Loop P-Tracking, Nominal F-18/HARV, SMC with Actuators 
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The next steps are to design the observers and hedge models. As introduced in 
Section 3.3, the individual Bode plots of the equivalent plant/observer transfer functions 
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^ a 

with unity feedback (-2- and ) are examined for varying observer speeds (with no 


cp 


cP 


hedging). Observer speeds of 100, 50, 20, 10, 5, and 1 rad/s are tried. The Bode plots 
for p s and (3 are shown in Figure 4-5 1 and Figure 4-52 respectively. 




Frequency (rad/s) 


Figure 4-51: Bode Plots, , F-18/HARV, Various Observer Speeds 
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Frequency (rad/s) 


Figure 4-52: Bode Plots, , F-18/HARV, Various Observer Speeds 

U cp 

As shown on the Bode plots, observer speeds of 50 and 100 rad/s in the p-loop result in 
an unstable system. Also, it is clear that the observer for the p-loop can not be slowed 
down enough to stabilize the system. Hedging is required for this system, at least in the 
P-loop. 

Before designing the hedge models, the transfer functions of the cross-coupled 
variables are examined The Bode plots of the open loop transfer functions for the 
nominal plant are given below. Figure 4-53 shows the stability axis roll rate transfer 
functions (roll-rate to roll-rate-command, and roll-rate to P-command). Figure 4-54 
shows the sideslip angle transfer functions ( P to P-command, and p to roll-rate- 
command) 
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Figure 4-53: Bode Plots, and — , F-18/HARV Nominal Vehicle 

Ucp U c/J 




Frequency (rad/s) 


Figure 4-54: Bode Plots, and F-I8/HARV Nominal Vehicle 
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Very strong coupling is seen The Bode magnitude plots of the cross-coupled transfer 
functions are of the same order of magnitude as the primary transfer functions. It is 
expected, therefore, that cross-coupled hedging models may be required. Following the 
design procedures outlined in Section 3.4, the following parameters are chosen for the 
observer and hedging models. 



Table 4-4: Observer and Hedge Models for Lateral Directional F-18/HARV 

The resulting design Bode plots are shown below. 
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Plant + Acuators + Observer + Hedge 

Hedge Signal 




Frequency (rad/s) 


Figure 4-55: Bode Plots, , F-18/HARV Hedge Design 
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Figure 4-56: Bode Plots, — , F-18/HARV Hedge Design 
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Note from Figure 4-55 that the p-loop is stable without hedging if the observer is run at 
10 rad/s. Hedging helps some but is not required for the nominal model. Figure 4-56 
shows that the P-loop is unstable without hedging. The hedging used on this channel 
stabilizes the system and provides a desirable feedback loop shape. 

Next, consider the cross-coupling hedge models. Figure 4-57 shows the p to 
P-command hedge design. As can be seen, hedging helps some but is probably not 
required. 



Figure 4-57: Bode Plots, 



(cross term), F-18/HARV Hedge Design 


Figure 4-58 shows the P to p-command hedge design. In this case, the effect of the 


p-command on P is strong enough to require hedging. The hedge model chosen provides 


a desirable feedback loop shape for this cross-coupled channel 
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Figure 4-58: Bode Plots, 



(cross term), F-18/HARV Hedge Design 


With the observer and hedge models designed, the system is now tested First the 
step responses for the inner loop are examined. The step responses for the inner loop 
with the given observer and hedge models are shown below. The roll rate command is a 
unit step, while the P command is a step ofO. 1. The two steps are commanded 
simultaneously in order to test the effects of cross-coupling. Figure 4-59 shows the 
results if the cross-coupled hedging models are neglected. Figure 4-60 includes the 
cross-coupled hedging models. 
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For the nominal model, the step responses in Figure 4-59 and Figure 4-60 indicate that 
neglecting the cross-coupled hedge models yields slightly better results It is seen later, 
however, that including the cross-coupled hedge models increases the system robustness 
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to variations in the actuator dynamics and reduces the activity in the P channel. Unless 
specified, all the following results use the coupled hedge models. 

Now exercise the full system with measurement noise and the pilot in the loop. 
The outer loop bank angle tracking performance for the nominal vehicle is shown in 
Figure 4-61. Performance is very good The large time delay between the commanded 
roll angle and the achieved roll angle is almost entirely due to the pilot. Note in 
Figure 4-62 that the inner loop tracking is very good. Of course, with the actuators and 
observer now included, the SMC does not give the almost perfect tracking achieved with 
the plant only There is still essentially decoupled control of p s and P, although P is not 
exactly zero. Considering the aggressive roll rates and the addition of state noise, P is 
acceptably small. Figure 4-63 and Figure 4-64 show the actuator deflections and rates, 
respectively. No position limits are reached, although there is some minor rate limiting 
in the rudder and yaw thrust vector. As shown in Figure 4-65 and Figure 4-66, Level I 
handling qualities are expected for the bank angle tracking task with no adverse PIO 
tendencies 
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Sideslip Angle, p (deg) Roll R at e, (deg/s) 


Time (sec) 


Figure 4-61: Outer Loop Tracking ( <j) ), Nominal F-18/HARV, with SMC, 
Actuators, Observers, Hedging, and Noise 




Figure 4-62: Inner Loop Tracking (p s and P), Nominal F-18/HARV, with SMC, 
Actuators, Observers, Hedging, and Noise 







Figure 4-63: Actuator Deflections, Nominal F-L8/HARV, with SMC, Observers, 

Hedging and Noise 
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Figure 4-64: Actuator Rates, Nominal F-18/HARV, with SMC, Observers, Hedging 

and Noise 
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Figure 4-65: Nominal F-18/HARV Handling Qualities Assessment, <j» Tracking Task 



Frequency (rad/s) 


Figure 4-66: Nominal F-18/HARV PIO Assessment, <|) Tracking Task 




224 


In order to compare performance of a failed vehicle with other types of designs, 
two different failures are defined. Failure #1 is defined as follows: failed plant, failed 
yaw and roll thrust vectors (max deflection = 15 deg, max rate = 30 deg/s, 50 ms time 
delay). This is a dramatic failure: 100% change in plant parameters, 25% reduction in 
effectiveness of all controls, a 50% reduction in thrust vector limits, and a sizeable time 
delay. 

The Simulink® simulation results are shown below. The simulation starts with a 
nominal vehicle followed by a failure which occurs at t = 25 sec. Figure 4-67 shows the 
outer loop bank angle tracking performance. Tracking is degraded after the failure, but 
the system remains stable. Figure 4-68 shows the inner loop tracking. P excursions 
change from negligible (before the failure) to a maximum value of about 1 deg (after the 
failure) Figure 4-69 shows the actuator deflections. No position limits are encountered, 
although the rudder is very near its limit. Figure 4-70 gives the actuator rates. The yaw 
thrust vector, roll thrust vector, and rudder are in almost constant rate saturation after the 
failure The differential tail and ailerons also reach their rate limits occasionally after the 
failure. This would seem to indicate that the system is on the edge of its maximum 
achievable performance. Figure 4-71 and Figure 4-72 indicate that a Level II handling 
qualities rating is expected and PIO tendencies are low. Given the severity of the failure, 
the demands of the flight condition, and the fact that all the control effectors are 
operating at their limits, it is doubtful that any control system could achieve Level I 
handling qualities for this condition. 
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Figure 4-67: Outer Loop Tracking ( <J> ), F-18/TIARV, with SMC, Actuators, 
Observer, Hedging, and Noise, Failure #1 at t = 25 sec 




Figure 4-68: Inner Loop Tracking (p s and P), F-18/HARV, with SMC, Actuators, 
Observer, Hedging, and Noise, Failure #1 at t = 25 sec 
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Figure 4-69: Actuator Deflections, F-18/HARV, with SMC, Observer, Hedging, 

and Noise, Failure #1 at t = 25 sec 
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Figure 4-70: Actuator Rates, F-18/HARV, with SMC, Actuators, Observer, 
Hedging, and Noise, Failure #1 at t = 25 sec 
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Figure 4-71: F-18/UARV, Failure #1, Handling Qualities Assessment, 

<j) Tracking Task 
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Figure 4-72: F-18/HARV, Failure #1, P10 Assessment, 4> Tracking Task 
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In order to compare these results to a well-known design approach, a classical loop 
shaped controller is designed. The first loop to design is the p-loop (with the p-loop 


open). The roll rate control transfer function selected is G (s) = 5.0 


^s 2 +3s + 2 25 a 
s + 0. 5s 


The resulting loop transmission for the control and plant is shown in Figure 4-73. This 
gives a crossover frequency of about co = 10 rad/s. 




Frequency (rad/s) 


Figure 4-73: Bode Plots, F-18/HARV Roll Rate Loop Transmission for Loop 

Shaped Design, P-Loop Open 


With the p-loop closed, the P-loop is designed: G„(s) = 63.0 


V + 0.9s + 0 2^ 
s 2 + 20s 


. The 


resulting loop transmission of the P controller and effective plant is shown in 


Figure 4-74. The crossover frequency for p is about © = 5 rad/s. 


Chapter 4: Design Procedure 


4.6 Application Example: F-18/HARV Lat-Dir MIMO Model 





230 




Figure 4-74: Bode Plots, F-18/HARV (3 Loop Transmission for Loop Shaped 

Design, p-Loop Closed 

This design results in excellent tracking in the nominal case. Although not shown here, 
the nominal tracking performance is comparable to the SMC design and gives Level l 
handling qualities. Next, this design is exercised with a system failure (Failure #1) at 
time t = 25 sec. Figure 4-75 shows the outer loop bank angle tracking performance The 
system goes unstable immediately after the failure. Figure 4-76 shows that the inner 
loop variables are unstable after the failure. Figure 4-77 and Figure 4-78 show that all 
the actuators are driven to their limits. 
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Figure 4-75: Outer Loop Tracking ( (j) ), F-18/HARV, with Loop Shaped 
Controller, Actuators, and Noise, Failure #1 at t = 25 sec 




Figure 4-76: Inner Loop Tracking (p s and P), F-18/HARV, with Loop Shaped 
Controller, Actuators, and Noise, Failure #1 at t = 25 sec 
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Figure 4-77: Actuator Deflections, F-18/HARV, with Loop Shaped Controller and 

Noise, Failure #1 at t = 25 sec 
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Figure 4-78: Actuator Rates, F-18/HARV, with Loop Shaped Controller and Noise, 

Failure #1 at t = 25 sec 
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This flight condition and failure provide an excellent example to illustrate the 
effects of cross-coupled hedging. In all the previous examples, the SMC design included 
hedging of the cross-coupled terms. It was noted that, in the nominal case, excluding the 
cross-coupled terms results in slightly better inner loop performance (see Figure 4-59 
and Figure 4-60). Flowever, when the actuator dynamics change (as they do in this 
failure case), the effects of the cross-coupled transfer functions can become significant. 
Consider a direct comparison of including the cross-coupled hedge models and 
excluding them Figure 4-79 shows the roll angle tracking performance for these two 
cases. Failure #1 occurs at t = 20 sec. 



Figure 4-79: Outer Loop Tracking ( <J> ), SMC with and without Cross-Coupled 

Hedging 
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Figure 4-80 shows the (3 tracking performance. 



Figure 4-80: Inner Loop Tracking ( P ), SMC with and without Cross-Coupled 

Hedging 

' Several points are worthy of mention. First of all, while the roll rate tracking is 
essentially the same for both cases in the nominal system, the hedge model with no 
cross-coupling results in an unstable system in the failed case. Further, the cross- 
coupled hedge model provides much better P tracking even in the nominal case. Sideslip 
excursions are an order of magnitude less. It is concluded that if the nominal system has 
strong cross-coupling, hedging of these cross-coupled channels should be considered. 

Next, to compare this SMC controller with another reconfigurable controller, a 
new failure condition is used. Failure #2 is defined as follows nominal plant, 90% 
reduction in yaw thrust vector effectiveness (gain = 0.1, position limits = ±3.0 deg, rate 
limits = ±6.0 deg/sec). This is a failure used in a previous work. 146 In order to directly 
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compare results, the same input signal is also used. The input is a series of steps filtered 


by 


100 


V i 


s 2 + 20s + 1 00 


s 2 + 2s + 1 


In the reference work, the reconfigurable design 


approach utilizes dynamic inversion, QFT, and least-mean square adaptive filtering 146 
The control laws take about 10 to 20 seconds to converge to their reconfigured states. 
The results in Figure 4-81 show the fully reconfigured QFT model 146 and the SMC 
design developed here Figure 4-82 shows the corresponding tracking for P with the 
SMC controller. 



Figure 4-81: Outer Loop Tracking ( <J> ), F-18/HARV, SMC vs Reconfigured QFT, 

Failure #2 at t = 0 sec 
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Figure 4-82: P Tracking, F-18/HARV, SMC, Failure #2 at t = 0 sec 


The inner loop tracking, actuator position, and actuator rate plots do not add anything to 
this comparison, so they are not shown As seen in Figure 4-8 1, the observer-based 
SMC with hedging provides very comparable tracking results to the fully reconfigured 
QFT system. 146 More importantly, the SMC accomplishes this with no adaptation time 
Additionally, the control architecture of reconfigurable QFT system 146 is considerably 
more complex than the design developed here. 
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4.7 Application Example: ICE MIMO 6-DOF Model 


tlirust vectoring 
nozzle 



leading edge flap 


spoiler slot deflector 
symmetric pitch flap 
elevon 


/ — all-moving 
ft wingtip 


Figure 4-83: Innovative Control Effectors (ICE) Aircraft 


The Innovative Control Effectors (ICE) aircraft is developed by Lockheed under 
an Air Force Research Laboratories (AFRL) sponsored program and is the vehicle of 
choice for many controls applications in the current literature. 18 ' 19 ’ 21 ' 25 ’ 48 ’ 36 ' 66,79 ' 801 10,147 It 
is a single engine, multi-role, supersonic, tailless fighter aircraft in the 38,000 lb (gross 
takeoff weight) class. It has a 65 degree sweep delta wing and an internal weapons 
carriage bay. The configuration incorporates all-aspect low observable technologies and 
is sized for a 1,100 nautical mile air-to-ground mission. 25 As illustrated in Figure 4-83, 
the conventional control effectors include elevons, symmetric pitch flaps, and outboard 
leading edge flaps. The innovative control effectors include pitch and yaw thrust 
vectoring, all moving tips, and spoiler slot deflectors. The all moving tips and spoiler 
slot deflectors have zero lower deflections limits. 

The static aerodynamic force and moment data were collected by NASA Langley 
Research Center and AFRL using wind tunnel tests with a 1/1 8 th scale model. 
Additional wind tunnel tests during Phase II of the ICE program provide updated data 
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for simulation models. 25 There are strong multi-axis effects and highly nonlinear 
interactions between the close-coupled control surfaces. The full nonlinear simulation of 
the ICE vehicle is proprietary; however, linearized models at various flight conditions 
(which do not include nonlinear effector interactions) are made available to this research 
directly from the NASA Langley Research Center. 

The flight condition used for this application design is. lg, wings-level, M = 0.3, 
h = 15,000 ft. The linearized model for this flight condition in state space is given by 
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The state vector is defined as. 
x T =[u w qb theta v pb rb pi,i] 
where 

u = velocity along body x-axis (ft/s), w = velocity along body z-axis (ft/s), qbsbody 
axis pitch rate (deg/s), theta = Euler pitch angle (deg), betaw = sideslip angle (deg), 
v = velocity along body y-axis (ft/s), pb = body axis roll rate (deg/s), rb = body axis yaw 
rate (deg/s), phi = Euler bank angle (deg). 
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The output vector is defined as 

t t = [vel alphaw qb theta betaw ps rs phi axcg aycg azcg] 
where 

vel = true airspeed (ft/s), alphaw = angle of attack (deg), qb = body axis pitch rate 
(deg/s), theta = Euler pitch angle (deg), betaw = sideslip angle (deg), ps = stability axis 
roll rate (deg/s), rs = stability axis yaw rate (deg/s), phi = Euler bank angle (deg), 
axcg = longitudinal acceleration at e g. (g’s), aycg = lateral acceleration at eg. (g’s), 
azcg = normal acceleration at c.g. (g’s) 

The control vector is defined as 

m t =[de3 de 1 3 de4 de5 del5 de9 del9 de2 del2 delO de20] 
where (see Figure 4-84) 

de3/del 3 = left/right elevon deflection (deg), de4 = symmetric pitch flap deflection 
(deg), de5/del5 = left/right all moving tip deflection (deg), de9/del9 = left/right spoiler 
slot deflector angle (deg), de2/del2 = lefl/right outboard leading edge flap deflection 
(deg), delO = pitch nozzle deflection (deg), de20 = yaw nozzle deflection (deg). 



Figure 4-84: ICE Control Deflection Definitions 
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The healthy actuator dynamics and limits are given in Table 4-5. 



Dynamics 

Amplitude 

Limit 

Rate Limit 

Elevon 

(40) (100) 

(s + 40) (s + 100) 

±30 deg 

1 50 deg/s 

Symmetric Pitch Flap 

(40) (100) 

(s + 40) (s + 100) 

±30 deg 

50 deg/s 

All Moving Tip 

(40) (100) 

(s + 40) (s + 100) 

[0, 60] deg 

1 50 deg/s 

Spoiler-Slot Deflector 

(40) (100) 

(s + 40) (s + 100) 

[0, 60] deg 

150 deg/s 

Leading Edge Flap 

(17.828) (100) 
(s+ 17.828) (s + 100) 

±40 deg 

40 deg/s 

Pitch/Y aw Thrust 
Vectoring 

(37.1 68) (41.3) 

(s + 37.186) (s + 41.3) 

±15 deg 

60 deg/s 


Table 4-5: ICE Actuator Dynamics and Limits 


A control system which tracks roll rate, angle of attack and sideslip is desired. 
The next step in the design process is to define a reference model. For this application, 
the following reference models are chosen: 


G.00 



100 

(s 2 + 20s + 100) 



100 

(s 2 + 20s + 100) 


G (s) = — (s) = — — — 

P P c (s 2 + 25s + 100) 


(4.22) 
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Next, the desired square system feedback architecture is defined. This is given in 
Figure 4-85. The inner loop has an angle of attack, roll rate, and a sideslip angle 
feedback loop which are controlled by the SMC. A zero sideslip is commanded; and a 
pilot closes an outer loop to control the angle of attack and Euler roll angle. 



Figure 4-85: ICE System Block Diagram 


The controller to be designed outputs commands for angle of attack, stability axis roll 
rate, and side slip angle. These demands are allocated to the eleven control effectors 
using a pseudo-inverse approach (see Eqn ( 1 .23 )). If the preferred actuator positions 
are assumed to be zero, the solution to the control allocation problem using a pseudo- 
inverse is 

B k = W;'B r (BW;'B r ) 1 (4.23) 

With weighting matrix, W u = diag[2,2,l,l,l,l,l,l,l,l,l], and the spoilers artificially 
removed from the pitch channel, the control distribution used in the design is 
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The input command profiles for angle of attack and roll angle are sum-of-sines functions 
defined in Eqn. ( 4. 13 ). The sum-of-sines parameters are those given in Table 4-2 The 
overall amplitude gain in Eqn. ( 4.13 ) for the a-channel is £ = 5 The overall amplitude 
gain for the 4>-channel is £=10. Additionally, the /= 1,7 sine functions for the <j>- 


channel have phase offsets of 


i — 

C 


in order to prevent correlation with the a-channel. 


The measured variables include the full output vector given above The measurement 
noise for each channel is assumed to be band-limited white noise filtered by 
100 


s + 20s + 100 


and multiplied by a gain. All variables except a and (3 have a noise gain 


of 0. 1414, resulting in an RMS noise value of approximately 0.07 deg/s. a and P have a 
noise gain of 0.0707, resulting in an RMS noise value of approximately 0.035 deg 

Assuming the inner loop SMC provides the desired model reference tracking, the 
apparent system models that the pilot “sees” are the model references. Therefore, the 
pilot models, 138 ’ 144 ’ 145 based on the reference models given, are shown below 
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Figure 4-86: Pilot Models for ICE Angle of Attack and Roll Angle Tracking Tasks 


The next step is to define the desired sliding manifolds. The system (without 
actuators) has a relative order of 1 for all three control variables. Therefore, the form of 
the sliding manifolds are 

o = K 0 e + K_ 1 jedx, e = {(«, -d), (p r -p), (p r -fi)f (4.25) 

The control laws, assuming the use of a boundary layer, then can be expressed in linear 
form as 


U C a(s) = K pa 

u cp(s) = K pp 

u c /»(s) = K p/J 


K 0 «s + K_ la 


V s J 

f K 0 „s + 


(4.26) 


This gives one zero to place and a gain to set for each loop during the loop-shaping 
design. Using a traditional sequential loop shaping technique, the parameters of 
Eqn. ( 4.26 ) are determined— with roll rate being the first loop to close, followed by a; 
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then p. The designed manifolds are given in Eqn. ( 4.27 ), and the design Bode plots are 
shown below. A crossover frequency of © = 1000 rad/s is set for each channel. 


U ca( s ) = 5000 


u cp (s) = 1000 


u c /j ( s ) = 5000 


^ s + io A 
V s J 

"s + 20^ 

v s J 

f s+10^ 

l s ) 


(4.27) 




Figure 4-87: ICE a-Loop Transmission of Compensator*PIant, 
p-Loop Open, p-Loop Open 
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Figure 4-88: ICE p-Loop Transmission of Compensator*Effective Plant, 

a-Loop Closed, P-Loop Open 




Frequency (rad/s) 


Figure 4-89: ICE P-Loop Transmission of Compensator* Effective Plant, 

a-Loop Closed, p-Loop Closed 
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The initial design of the SMC is complete and is shown in Figure 4-90 



Figure 4-90: SMC Controller for ICE 


The fifth step in the design process is to verify that the SMC is working properly. 
Either replace the saturation elements in Figure 4-90 with a signum function or make the 
boundary layer very small. Figure 4-91 and Figure 4-92 show the SMC inner loop 
tracking of the control variables and the SMC control effort. This Simulink® simulation 
is run with signum functions in the SMC and a time step of At = 0.0005 sec. A system 
plant failure occurs at t = 10 sec. Although the P input will be a constant zero for the 
final simulation, a simple sine wave is used here to verify decoupled tracking of a, p s 
and P 
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As expected, the SMC performs very well and is invariant to the system parameter 
changes. The control output shows the classic high frequency switching. Note, also that 
the controller achieves decoupled tracking of a, p s and P, as desired. 

Next, the boundary layer is increased until a continuous control signal is achieved. 

For this model, e = 1 is chosen for all three channels. The resulting time history of the 
inner loop tracking of a, p s and P is shown in Figure 4-93. The corresponding SMC 
control efforts are shown in Figure 4-94. Again, the performance is excellent, even in 
the face of a large system failure. Also, the control signal is now continuous. Note the 
large spike in the control signal at the time of failure and the increased overall magnitude 
of the required control signal after the failure 



-20 I i i J i _j I 

0 5 10 15 20 25 30 
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Time (sec) 

Figure 4-93: ICE Inner Loop SMC Tracking with Boundary Layer 
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Time (sec) 

Figure 4-94: ICE SMC Control Outputs with Boundary Layer 


The SMC design is now complete, the resulting tracking performance is 
remarkably good, and the system is (nearly) invariant to system parameter changes The 
next step is the inclusion of the actuators. Figure 4-95 shows the inner loop tracking 
results with the actuators now included. As expected, the nominal system goes unstable. 
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Tim e (sec) 


Figure 4-95: Inner Loop Tracking, Nominal ICE, SMC with Actuators 


The next steps are to design the observers and hedge models. As introduced in 
Section 3.3, the individual Bode plots of the equivalent plant/observer transfer functions 


~ ~ O 

with unity feedback and ) are examined for varying observer speeds (with 

U U U o 

ca cp cp 


no hedging). Observer speeds of 100, 50, 20, 10, 5, and I rad/s are tried. The Bode 
plots for a, p s and P are shown in Figure 4-96 , Figure 4-97, and Figure 4-98 
respectively 
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Figure 4-96: Bode Plots, 


ICE, Various Observer Speeds 



Frequency (rad/s) 


Figure 4-97: Bode Plots, 


P 


ICE, Various Observer Speeds 
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These plots show that the a-channel observer must run at 10 rad/s or less; the p-channel 
observer can run as fast as 1 00 rad/s, and the P-channel must run at 5 rad/s or less It 
appears that hedging will help some in the a-channel, is not necessary in the p-channel, 
and is probably required in the P-channel. Select 5 rad/s, 10 rad/s, and 1 rad/s for the 
a, p, and P observers respectively and check the step responses without hedging Step 
inputs of u a = 1.5, Up = 1.0, and up = 0 5 are used. The results are shown in Figure 4-99. 
It is clear that hedging is required in the P-channel, and there appears to be undesirable 
cross-coupling between p s and p. 
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Figure 4-99: Inner Loop Step Responses, ICE, SMC with Observers 

Before designing the hedge models, the transfer functions of the cross-coupled 
variables are examined. The Bode plots of the open loop transfer functions for the 
nominal plant are given below Figure 4-100 shows the angle of attack transfer 
functions (a to a-command, a to p-command, and a to (3-command) Figure 4-101 
shows the stability axis roll rate transfer functions (p to a-command, p to p-coinmand, 
and p to P-command). Figure 4- 102 shows the sideslip angle transfer functions (P to 
a-command, p to p-command, and P to p-command). 
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Figure 4-102: Bode Plots, 





ICE Nominal Vehicle 


(X 

Figure 4-100 shows some coupling between ^-command and a. Therefore, a — — term 

u c p 

will be included in the hedge model. The influence of roll rate command on angle of 
attack is a couple orders of magnitude less than angle of attack command, so that cross- 
coupling term will be neglected. Figure 4-101 shows very strong coupling between 
(i-command and p s . In fact the cross-coupled term is an order of magnitude greater than 


the primary transfer function. The — — term must be included in the hedge model 


u 


Although Figure 4-102 shows only a weak coupling between P and roll rate command, a 
term will be included in the hedge model All other cross-coupling terms are 


P h 


cp 
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negligible and will be ignored. Following the design procedures outlined in Section 3.4, 
the following parameters are chosen for the observer and hedging models. 


Channel 


Observer Poles 


A = -10.0, -10 1 -10.2, -103, -10.4, -10 5, -10.6, -10.7 


A = -40.0, -40. 1 , -40.2, -40 3, -40.4, -40.5, -40.6, -40 7 


A, = -1.0, -1. 1, -1.2, -1.3, -1.4, -1.5, - 1 .6, -1.7 


Hedge Model 


s + lOAs + 4s + 4 


s + 50 As + 4s + 4 


s + 5 A s +3s + 2 


0.0001 — a cross term 

s + lOAs 2 +4s + 4 J 


Y 64 A 

0.2 — p s cross term 

ls + 50 As + 16s + 64 J 


0.02 — p cross term 

. s + 5 A s +2s + ly 


Hedge Gain 


Table 4-6: Observer and Hedge Models for ICE 


The design Bode plots for each channel are given below. 


Chapter 4: Design Procedure 


4. 7 Application Example: ICE MIMO 6-DO F Mode! 





Phase (deg) Magnitude (dB) Phase (deg) Magnitude (dB) 


Hi 



Nominal Plant 


\r.:- : : 


Plant + Actuators +• Observers 

i- i-i- 

::l:. . l 


Plant + Actuators + Observers + Hedge 

— - Hedge Signal 

r'rr 

rrH — 


■ 

-Ml 

; H iii; 

1 

L •- ■_ 

uii j. j 


l Ini' l ■ Vi ■, t I ■ T 


10 ° 10 ' 

Frequency (rad/s) 


Figure 4-103: Bode Plots, 


, ICE Hedge Design 
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Figure 4-104: Bode Plots, — — , ICE Hedge Design 
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With the observer and hedge models designed, the system is now tested. First the 
step responses for the inner loop are checked. The results with the given observer and 
hedge models are shown below. Again, step inputs of u a = 1.5, u p = 1.0, and up = 0.5 are 
used. The three steps are commanded simultaneously in order to test the effects of 
cross-coupling 



01 23456789 10 

Time (sec) 


Figure 4-109: Inner Loop Step Responses, ICE, SMC with Observers and Hedging 

It is possible to achieve much better inner loop tracking by shaping the hedge models 
and adjusting the hedge gains; however, some design trade-offs are made to enhance 
robustness to actuator damage. The hedge signals are designed with a fairly broad range 
of -20 dB slope and hedge gains that are a little higher than needed for the nominal case. 
This results in slightly slower tracking response and some overshoot in the unfailed 
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vehicle. Again, all the design decisions are made using the Bode plots in Figure 4-103 
through Figure 4-108 and the step responses in Figure 4-109. 
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Finally, exercise the full system with measurement noise, system failure, and the 
pilot in the loop The Simulink® simulation is run with an ODE2 solver using a fixed 
time step of At = 0.0005 sec. Vehicle failure occurs at t = 20 sec. Vehicle failure is 
defined as: 

• Plant failure (A-matrix is multiplied by 2 — with the exception of the 
elements describing kinematic relationships; B-matrix is multiplied by 
0 75) 

• All actuators experience a 50 ms time delay 

• Left elevon rate limits are reduced from 150 deg/s to 10 deg/s, position 
limits are reduced from +30 deg to ±15 deg 

• Symmetric pitch flap undamped natural frequency is reduced from 63 rad/s 
to 1 0 rad/s 

• Left leading edge flap is jammed at +5 deg 

• Pitch nozzle actuator undamped natural frequency is reduced from 39 rad/s 
to 1 0 rad/s 

• Yaw nozzle actuator undamped natural frequency is reduced from 39 rad/s 
to 10 rad/s 

The resulting outer loop tracking is shown in Figure 4-1 10 Figure 4-1 11 shows the 
inner loop SMC tracking 
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Figure 4-1 10: Outer Loop Tracking ( a, <j>, and P ), ICE, with SMC, Actuators, 
Observers, Hedging, and Noise, Failure at t = 20 sec 




Figure 4-111: Inner Loop Tracking ( a, and p s ), ICE, with SMC, Actuators, 
Observers, Hedging, and Noise, Failure at t = 20 sec 
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Before the failure, the controller does a very good job. Desirable, decoupled tracking of 
the controlled variables is seen. Roll rate tracking could be better, but as explained 
earlier, this is a design tradeoff that is made to enhance robustness As shown in 
Figure 4-1 12 and Figure 4-1 13, Level 1 handling qualities are predicted and PIO 
tendencies are mild for the nominal vehicle. 




Figure 4-112: a-Tracking Task HQ and PIO Predictions, ICE Nominal System 




Figure 4-113: ({(-Tracking Task HQ and PIO Predictions, ICE Nominal System 

After system failure, tracking is noticeably degraded, but the vehicle remains stable. 
The actuator deflections are not shown since there are eleven. Upon examination, 
however, it is seen that nearly all the actuators are in almost constant rate saturation after 
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the failure. As seen in Figure 4-114 and Figure 4-115, Level I handling qualities are still 
predicted for the a-tracking task , although PIO susceptibility is increased Level II 
handling qualities are predicted for the <j)-tracking task, and PIO tendencies are 
moderately strong. 




Figure 4-114: a-Tracking Task HQ and PIO Predictions, ICE Failed System 




Figure 4-1 15: ^-Tracking Task HQ and PIO Predictions, ICE Failed System 


In viewing these results, take into account a several factors. First of all, this is a 
significant system failure, and it is remarkable that the controller is even able to maintain 
stability. Second, the input maneuvers are very aggressive. After a failure of this 
magnitude, a pilot would not be commanding such a demanding profile. Third, the HQ 
and PIO level predictions assume a static pilot model In actuality, a real pilot would 
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compensate for the perceived change in vehicle dynamics. Therefore, the results shown 
above are conservative. 

In order to compare these results with a baseline controller, a classic loop-shaped 
controller is designed. Using a sequential loop-shaping technique, the following 
controller is selected 


G a (s) = 56 


V + 1.7s + 0.52^ 
s 2 + 32s , 


G 0 (s) = 10 


( s 2 + 0.4s + 0.04 1 


Gp(s) = 50 


^s 2 + 2s + I ^ 
s 2 + 50s 


The design Bode plots are given below for reference. 


(4.28) 



Figure 4-116: ICE Loop Shaped Controller, a-Loop Transmission of 
Compensator*Plant, p-Loop Open, P-Loop Open 
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igure 4-117: ICE Loop Shaped Controller, p-Loop Transmission of 
Compensator*Effcctive Plant, a-Loop Closed, P-Loop Open 
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Figure 4-118: ICE Loop Shaped Controller, P-Loop Transmission of 
Compensator*Eftective Plant, a-Loop Closed, p-Loop Closed 
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This design results in very good tracking of the nominal system. Level I handling 
qualities are predicted for both the a and <j> tracking tasks. Figure 4-119 shows the outer 
loop tracking using this controller. The same failure as given above occurs at t = 20 sec 
As seen, the system goes unstable almost immediately after failure 
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Figure 4-119: Outer Loop Tracking ( a, 4>, and P ), ICE, with Loop-Shaped 
Controller, Actuators, and Noise, Failure at t = 20 sec 

In conclusion, a MIMO observer-based SMC with hedging is demonstrated for a 
challenging, highly coupled, six degree-of-freedom system. While providing excellent 
decoupled tracking of the controlled variables, it is shown to be extremely robust to both 
system plant variations and actuator failures and is seen to be superior to a traditional 
design. 
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Chapt er 5 

Conclusions and Recommendations 

5.1 Summary 

Sliding mode controllers are extremely attractive for use in a reconfigurable flight 
control system. Their invariance to matched uncertainty enable them to adapt to huge 
unmodeled plant failures without any changes to the control law. In fact, if they could 
handle parasitic dynamics, the job would be done. Unfortunately, all real mechanical 
systems with limited bandwidth actuators make the application of SMC more difficult. 
Using asymptotic observers with SMC has been shown to help mitigate the adverse 
effects of parasitic dynamics, and a design method for choosing observer gains has been 
presented. A form of model reference hedging has been shown to be equivalent to an 
SMC pre-filter which helps shape the feedback loop— thus partially removing the effects 
of actuator dynamics. A complete design procedure has been presented which 
incorporates a frequency domain approach to select 1) the sliding manifold, 2) the 
observer eigenvalues, and 3) the hedging dynamics. Finally, three aircraft applications 
have been demonstrated, each showing exceptional robustness to both system parameter 
changes and actuator variations. As with most research, however, this work has 
introduced as many questions as answers. There are many directions that future research 
can take. Below are several key areas 
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5.2 Potential Future Research 

5.2.1 Robustness Proofs 

The control schemes used are not fail-proof, and a quantitative measure of their 
robustness has not been found. Robustness proofs for true sliding mode controllers can 
be obtained using Lyapunov analyses, however, even these are problematic for practical 
problems. It is not clear which parameters are “matched” and which are “unmatched ” 
It appears that most of the parameters are partially matched, and it also appears that 
matching conditions change depending on frequency. It is clear that the observer-based 
SMC retains considerable robustness, and it is clear that the observer-based SMC is 
performing some kind of internal adaptation, however, no formal proofs of stability and 
robustness are known. This is a major research area which requires serious 
consideration 

5.2.2 Full Flight Envelope Issues 

The examples given are for single point flight conditions, and no effort is made to 
extend the analysis to the full flight envelope. One advantage to using a sliding mode 
controller is that the control law does not need to change for a class of system plants. As 
long as the SMC gain, p, is chosen large enough, the SMC can easily handle varying 
flight conditions. There are certainly conditions for which the sliding manifold itself 
would require adjustment, but most conditions would probably require no changes to the 
SMC. Unfortunately, a pure SMC can not be directly implemented due to the actuators, 
and the methods developed here require a nominal model of the plant in the observer. It 
is shown that the observer-based SMC is robust to changes in the plant parameters; 
therefore, changes in flight condition can be handled to a certain extent. However, if the 
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aircraft then experiences damage in an off-design flight condition, stability is doubtful 
If, however, the nominal model in the observer changed with the flight condition, the 
observer-SMC system should perform very well. This could be accomplished by using a 
scheduled observer with various pre-determined nominal plant models as a function of 
flight condition (similar to controller gain scheduling). It could also be accomplished in 
a more sophisticated manner using some type of adaptive observer. This would 
probably require some form of parameter identification and would bring with it all the 
complexities of that problem However, since flight condition parameters change 
relatively slowly and the observer-SMC is robust to modeling inaccuracies, some of the 
difficulties with parameter identification would be less critical here. These are questions 
which deserve further investigation. 

5.2.3 Adaptive Hedging 

It is observed that varying the hedge gain in the hedging scheme introduced here 
can dramatically affect system stability and performance. In general, increasing the 
hedge gain reduces the control activity and introduces a perceived phase lag to the pilot 
Reducing the hedge gain increases performance to a point, but further reduction can lead 
to instability due to the parasitic dynamics It would be desirable to find an adaptive 
scheme by which the performance, robustness, and control effort tradeoffs could be 
made on-line Several attempts were made to implement an adaptive scheme, but none 
of them were entirely successful. This is because there are at least two fundamental 
problems with the implementation. First, the system is very sensitive to dynamic 
changes in hedge gain. Even if an optimum hedge gain could be calculated on-line, it 
does not appear to be possible to change the hedge gain quickly enough to counter a 
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system failure because a dynamic hedge gain drives the system unstable. Second, as 
seen in the application examples, varying hedge gain alone is not sufficient Changes in 
actuator dynamics may require changes to the hedge plant and filter — not just the hedge 
gain. The issue of an adaptive hedge gain is an open topic which may deserve further 
study. 

5.2.4 Dynamic Hedged Observer 

Recognizing that the purpose of hedging is to remove the effects of actuator 
dynamics, a completely different scheme which tries to capture the original spirit of 
hedging was tried. Here, a method (which still makes use of an observer) attempts to 
completely hide the effects of the actuators from the SMC. This scheme requires 
knowledge of the actuator positions. Although not specifically mentioned before, a soft 
constraint on all the work done here is that no knowledge of actuator positions should be 
assumed. If this constraint is relaxed, the following architecture for an observer can be 
realized. 



Figure 5-1: Hedged Observer 
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In this model, u is the vector output of the SMC controller (which is fed to the 
actuators), del is the vector of actual actuator positions; y is the vector of measured 
states, and xhat is the vector of observed states which are fed back to the controller. The 
input to the block labeled Nominal Plant is the difference between the commanded 
control and the achieved control. Its output, then, represents the deficit in the system 
state due to the unachieved control input. Usually, in an observer, the measured states 
are compared with the observed states in a feedback loop In this model, the observer 
feedback loop also contains the effects of the actuator dynamics. This means the output 
of the observer block represents the state vector of the system if there were no actuators. 
This signal is the one that is passed to the controller, and the effects of the actuator 
dynamics are completely hidden from the SMC. Since the SMC does not “see” the 
actuators, the observer can be run very fast to ensure excellent convergence of the 
observer error. This scheme works very well — the actuators are hidden from the SMC 
and the SMC gives its expected near-perfect tracking and invariance to changes in 
system parameters. For example, the F-18/HARV longitudinal model from Chapter 4 is 
run with this new scheme. Failure #2 occurs at time t = 10 sec. Results are shown in 
Figure 5-2 Although not detectable on this plot, q r and qt, a t are coincident, i.e , the 
observed state perfectly tracks the reference signal before and after failure. However, 
the actual state does not track the reference signal after the failure. Although unexpected 
at first, this result makes sense. The SMC is tracking the system it sees. Nowhere do the 
actual states get fed back to the controller. Even though this approach does not solve the 
ultimate problem, it does appear to be a very promising first step. A hedged observer 
has been used to reconstruct the states of the system as it would exist if there were no 
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Figure 5-2: inner Loop Pitch Rate Tracking, F-18/HARV, Dynamic Observer 

Hedging, Failure #2 at t = 10 sec 


actuators. The near-invariance of the SMC is recaptured. A potential solution to 
complete this scheme is to wrap an outer loop around this inner loop, feeding back actual 
measured states through an additional controller. The form that this controller should 
• take is not clear Introduction of another sliding mode controller is tempting, however, 
the apparent system that this controller must control has a relative order of plant plus the 
actuators. This presents the same problem as including the actuators in the model in the 
first place. For a system with second order actuators, two derivatives of the state 
variable would be required in the calculation of the sliding manifold. This is completely 
impractical. Traditional linear controllers could be used. Even these, however, present 
problems because some derivative action is going to be required due to the loop shape of 
the apparent inner loop system. Therefore, noise is a primary concern. Initial tests using 
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an inner loop SMC with a dynamic hedged observer combined with an outer loop 
controller show some promise However, there are implementation issues which require 
further work 
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5.2.5 Dynamic Control Allocation 

All the control allocation done in this work has been static. It seems clear that 
additional performance could be realized if dynamic control allocation is employed 
This, however, requires knowledge of several system elements that are not easily 
obtained. The system B matrix must be known. This would require some form of 
parameter identification. Also, the actuator positions must be known. Additionally, 
actuator rate and position limits must be known. For these reasons, dynamic control 
allocation has been avoided. However, other research efforts utilize dynamic control 
allocation, and incorporation of it in this work may warrant further investigation. 

5.2.6 Actuator Position Limiting Issues 

The issue of actuator position limiting has not been satisfactorily addressed. All 

the systems used in this work have redundant control effectors. Therefore, if one 
effector saturated in position, the system did not go unstable because the redundant 
effector was able to compensate for the loss of control effort. If there were no redundant 
effectors, or if all the effectors reached position limits, the SMC would fail. It still 
seems that the initial concept of hedging the reference model may be a solution to this 
problem. This limits the amount of the commanded signal by the amount that is 
unattainable due to the saturated actuators. No successful implementation of this method 
was obtained Whether by model reference hedging or some other scheme, the issue of 
control position saturation requires much work. 
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